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Abstract 



Despite the fact that pulsars have been observed for almost half a century, until now many 
questions have remained unanswered. One of the fundamental problems is describing the 
physics of pulsar radiation. By trying to find an answer to this fundamental question 
we use the analysis of X-ray observations in order to study the polar cap region of radio 
pulsars. The size of the hot spots implies that the magnetic field configuration just above 
the stellar surface differs significantly from a purely dipole one. By using the conservation 
of the magnetic flux we can estimate the surface magnetic field as of the order of 10^^ G. 
On the other hand, the temperature of the hot spots is about a few million Kelvins. Based 
on these two facts the Partially Screened Gap (PSG) model was proposed to describe the 
Inner Acceleration Region (lAR). The PSG model assumes that the temperature of the 
actual polar cap is equal to the so-called critical value, i.e. the temperature at which the 
outflow of thermal ions from the surface screens the gap completely. 

We have found that, depending on the conditions above the polar cap, the generation 
of high energetic photons in lAR can be caused either by Curvature Radiation (CR) or 
by Inverse Compton Scattering (ICS). Completely different properties of both processes 
result in two different scenarios of breaking the acceleration gap: the so-called PSG-off 
mode for the gap dominated by CR and the PSG-on mode for the gap dominated by ICS. 
The existence of two different mechanisms of gap breakdown naturally explains the mode- 
changing phenomenon. Different characteristics of plasma generated in the acceleration 
region for both processes also explain the pulse nulling phenomenon. Furthermore, the 
mode changes of the lAR may explain the anti-correlation of radio and X-ray emission in 



very recent observations of PSR B0943-I-10 (Hermsen et al. , 2013). 



Simultaneous analysis of X-ray and radio properties have allowed to develop a model 
which explains the drifting subpulse phenomenon. According to this model the drift takes 
place when the charge density in lAR differs from the Goldreich- Julian co-rotational 
density. The proposed model allows to verify both the radio drift parameters and X-ray 
efficiency of the observed pulsars. 



Streszczenie 



Pominio, ze pulsary sej; badane juz od prawie pol wieku, do dzisiaj nie udalo si? znalezc 
odpowiedzi na wiele pytan. Jednym z fundament alnych problemow jest opis fizyki promi- 
eniowania pulsarow. Probuj^c znalezc odpowiedz na to fundament alne pytanie, wykorzys- 
tujemy analiz? obserwacji rentgenowskich w celu badania obszaru czapy polarnej pulsarow. 
Rozmiar obserwowanych gor^cych plam wskazuje, ze konfiguracja pola magnetycznego na 
powierzchni gwiazdy rozni si? znacznie od pola czysto dipolowego. Wykorzystujq^c prawo 
zachowania strumienia magnetycznego mozemy oszacowac sil? pola magnetycznego w ob- 
szarze czapy polarnej, ktore dla obserwowanych pulsarow jest rz?du 10^^ G. Z drugiej 
strony obserwowana temperatura gorq^cej plamy jest rz?du kilku milionow kelwinow. Opi- 
erajq^c si? na tycli dwocli faktach wykorzystujemy model cz?sciowo-ekranowanej przerwy 
akceleracyjnej (z ang. Partially Screened Gap - PSG), aby opisac wewn^trzuq^ przerw? 
akceleracyjnq, (z ang. Inner Acceleration Region - lAR). Model PSG zaklada, ze temper- 
atura czapy polarnej jest bliska do tak zwanej wartosci krytycznej tzn. takiej przy, ktorej 
termiczny odplyw jonow z powierzchni w pelni ekranuje przerw? akceleracyjnq,. 

W zaleznosci od warunkow jakie panujci, w obszarze czapy polarnej, mechanizmem 
odpowiedzialnym za generowanie wysokoenergetycznych fotonow w lAR moze bye pro- 
mieniowanie krzywiznowe (z ang. Curvature Radiation - CR) lub odwrotne rozpraszanie 
Comptona (z ang. Inverse Compoton Scattering - ICS). Calkowicie rozne wlasciwosci 
obu tych procesow prowadza^ do sytuacji, w ktorej mozemy wyroznic dwa scenariusze 
zamkni?cia przerwy akceleracyjnej: tzw. PSG-off dla przerwy zdominowanej przez pro- 
mieniowanie CR, oraz tzw. PSG-on dla przerwy zdominowanej przez ICS. Istnienie dwoch 
roznych mechanizmow zamkni?cia przerwy w naturalny sposob tlumaczy zjawisko zmi- 
any trybu promieniowania pulsarow (z ang. mode-changing). Rozna charakterystyka 
plazmy generowanej w obszarze akceleracyjnym dla obu tych trybow tlumaczy zjawisko 
sporadycznego braku pojedynczych pulsow (z ang. pulse nulling) w obserwacjach radiow- 
ych. Co wi?cej zmiana trybu w jakim pracuje przerwa akceleracyjna moze zostac powiq^z- 
ana z antykorelacjq promieniowania radiowego i rentgenowskiego wykazanq, w ostatnich 



obserwacjach PSR B0943+10 (Hermsen et al. 2013) 



Jednoczesna analiza wlasciwosci promieniowania rentgenowskiego i radiowego poz- 
wolila na opracowanie modelu dryfuj^cych skladowych pulsu pojedynczego (z ang. sub- 
pulses). Model ten zaklada, ze dryf jest wynikiem roznicy g?stosci ladunku w lAR w 
stosunku do g?stosci korotacji. Proponowany model pozwala zarowno na weryfikacj? 
wyznaczonych parametrow dryfu oraz na weryfikacj? np. efektywnosci promieniowania 
rentgenowskiego. 
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Introduction 



The history of neutron stars began in the early 1930s when Subrahmanyan Chandrasekhar 
calculated the critical mass for a white dwarf. As soon as the mass of a white dwarf exceeds 
the critical value (e.g. due to accretion of matter from a companion star) it collapses and a 
neutron star is formed. Chandrasekhar estimated that the critical mass was approximately 
1.4 solar masses (M©). Even before James Chadwick's discovery of neutrons (1932), Lev 



Landau anticipated the existence of neutron stars by writing about stars in which "atomic 



nuclei come in close contact, forming one gigantic nucleus". In 1934 Baade and Zwicky 
proposed that the "supernova process represents the transition of an ordinary star into a 



neutron star". Five years later Oppenheimer and Volkoff ( 1939 ), using the work of Tolman 



(1939), computed an upper bound on the mass of a star composed of neutron-degenerate 



matter. They assumed that the neutrons in a neutron star form a cold degenerate Fermi 
gas which leads to an upper bound of approximately O.7M0. Modern estimates of the 



critical mass for neutron stars range from approximately 1.5 Mq to 3 Mq (Bombaci, 1996). 
This uncertainty reflects the fact that the equation of state for extremely dense matter is 
not well known. Let us note that the radius of a neutron star should he R ^ 10 km. On 
the other hand nobody expected to detect any emission from neutron stars due to their 
small size and the lack of theoretical predictions about any radiation processes, except 
for thermal radiation. Thus, it took almost forty years to detect emission from a neutron 
star. 

The breakthrough came on 28 November 1967 with the radio observations that were 
performed by Jocelyn Bell-Burnell and Anthony Hewish. They observed radio pulses 
separated by 1.33 seconds. The world "pulsar" was adopted to reflect the specific property 
of these celestial objects. The suggestion that pulsars were rotating neutron stars was put 



forth independently by Gold (1968) and Pacini (1968), and was soon proved beyond a 



reasonable doubt by the discovery of a pulsar with a very short (33-millisecond) pulse 
period in the Crab nebula. It was suggested that this pulsar powers the activity of the 



nebula (Pacini, 1968). Nearly 2000 pulsars have been found so far. Observations of pulsars 



provide valuable information about neutron star physics, general relativity, the interstellar 

medium, celestial mechanics, planetary physics, the Galactic gravitational potential, the 

magnetic field and even cosmology. Studying neutron stars is therefore a very broad issue 

and it is beyond the scope of this thesis to describe the current status of the theory of 

neutron stars or pulsar population studies in detail. We rather refer the reader to the 
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literature (Michel 1991; Meszaros 1992 Glendenning 1996 Weber, 1999) and provide 



only a basic theoretical background that is relevant to the subject of this thesis. 



Following the ideas of Pacini (1968) and Gold (1968) radio pulsars can be interpreted 



as rapidly spinning, strongly magnetised neutron stars radiating at the expense of their 
rotational energy. Neutron stars consist of compressed matter with density in its core 
exceeding nuclear density pnuc = 2.8 x lO^^gcm"^. Direct and accurate mass meas- 
urements come from timing observations of binary pulsars and are consistent with a 
typically assumed neutron star mass M ^ 1.4 Mq. Most models predict a radius of 
i? ~ 10 km, which is consistent with the theoretical upper and lower limits. However, the 
measurements of neutron star radii are much less reliable than the mass measurements. 
Therefore, the moment of inertia for these canonical values (M = 1.4 Mq, R = 10 km) 
/ ~ (2/5) MR"^ ^ lO^^gcm^ may be uncertain by ~ 70%. The increase rate of a pulsar 
period, P = dP/dt, is related to the rate of rotational kinetic energy loss (spin-down 
luminosity) E = L^b = Att^IPP^^. In most cases only a tiny fraction of E can be 
converted into radio emission. The efficiency, Xradio = -^radio/-E, in the radio bands is 
typically in the range of ~ 10^^ — 10^^. It is assumed that the bulk of the rotational 
energy is converted into magnetic dipole radiation. The expected evolution of the angular 
velocity {Q = 27r/P) of a rotating magnetic dipole can be described as ^2 ~ fi", and the 
breaking index is n = 3 for the pure dipole radiation. Indeed, the observed values of 



the breaking index (e.g. Becker, 2009) confirm the above statement, e.g.: for the Crab 
n = 2.515 ± 0.005, for PSR B1509-58 n = 2.8±0.2, for PSR B0540-69 n = 2.28±0.02, for 
PSR J1911-6127 n = 2.91 ± 0.05, for PSR J1846-0258 n = 2.65±0.01, and for the Vela 
pulsar n = 1.4±0.2. On the other hand the observations of pulsar wind nebulae suggest 
that a significant fraction of the pulsar rotational energy is carried away by a pulsar wind. 
Furthermore, recent observations of high energy radiation from pulsars show that signific- 
antly more energy is radiated in the form of X-rays and 7-rays than in the form of radio 



emission (e.g. Abdo et al. 2010). Thus, pure magnetic breaking does not provide full 



information about the physical processes that take place in the pulsar magnetosphere. 

Despite the fact that pulsars have been observed for almost half a century, many 
questions still remain unanswered. One of the fundamental problems concerns the physics 
of pulsar radiation. Radio observations alone cannot point to the model (e.g. vacuum gap, 
slot gap, outer gap, free outfiow, etc.) that correctly describes the source of pulsar activity. 
Observations carried out by relatively new high-energy instruments, e.g. Chandra, and 
XMM-Newton, significantly extended the spectra over which we can study pulsars and 
their environments. There is no consensus about the origin of pulsar X-ray emission 



(Michel, 1991 ). We can distinguish two main types of models: the polar gap and the outer 
gap. The polar gap models suggest that the emission region is located in the vicinity of 
the neutron star polar caps, while the outer gap models assume that particle acceleration 



Vlll 



and X-ray emission take place close to the pulsar light cylindeiQ In both types of models 
high energy radiation is generated by relativistic particles accelerated in charge-depleted 
regions, while the high energy photons are emitted by means of Curvature Radiation 
(CR), Synchrotron Radiation (SR) and Inverse Compton Scattering (ICS). Both models 
are able to interpret existing observational data. 



In this thesis we will use the Partially Screened Gap (PSG) model (Gil et al. , 2007a). 
The PSG model assumes the existence of the Inner Acceleration Region (lAR) above 
the polar cap (a region penetrated by the open field lines) where the electric field has 
a component along the magnetic field. In this region particles (electrons and positrons) 
are accelerated in both directions: outward and toward the stellar surface. Consequently, 
outflowing particles are responsible for generation of magnetospheric emission (radio and 
high-frequency) while the backflowing particles heat the surface and provide the required 
energy for thermal emission. The PSG model is an extension of the Standard Model 



developed by Ruderman and Sutherland (1975) and takes into account the thermionic 



ion flow from the stellar surface heated up to a high temperature (a few million Kelvins) 
by the backstreaming particles. In such a scenario an analysis of X-ray radiation is an 
excellent method of obtaining insight into the most intriguing region of the neutron star. 



-'^The light cylinder with radius i?LC = cP/2ti is defined as a place where the azimuthal velocity of the 
co-rotating magnetic field lines is equal to the speed of light (c) 
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Chapter 1 

X-ray emission from Radio Pulsars 



1.1 Brief historical overview 

X-ray photons can only be detected by telescopes operating at high altitudes or above the 
Earth's atmosphere, thus detectors should be mounted on high-flying balloons, rockets or 
satellites. The first (i.e. carried out from space) X-ray observations were performed by 
a team led by Herbert Friedman in 1948. The team estimated the luminosity of X-ray 
radiation from the solar corona. They found that X-ray luminosity is weaker by a factor 
of 10^ than luminosity in the optical wave range. Up until the early 1960s it was widely 
believed that all other stars should be so faint in the X-rays that their observations 
would be hopeless. The situation changed in 1962 when a team led by Bruno Rossi 
and Riccardo Giacconi, when trying to find fluorescent X-ray photons from the moon, 
accidentally detected X-rays from Sco X-1. Subsequent flights launched to confirm these 
first results detected Tau X-1, a source in the constellation Taurus which coincided with 



the Crab supernova remnant (Bowyer et al. , 1964). The search for similar sources became 



a source of strong motivation for the further development of X-ray astronomy. 



Before the first direct detection of a neutron star by Hewish et al. (1968), it was 
predicted that neutron stars could be powerful sources of thermal X-ray emission due 
to a high surface temperature (Tg). The expected value of the surface temperature was 



estimated as Tg ~ 1 MK (Chiu and Salpeter, 1964 Tsuruta, 1964). The first X-ray 



observations of isolated neutron stars [j were initiated by the Einstein Observatory, which 
was launched by NASA in 1978. Using a high-resolution imaging camera sensitive in the 
0.2 — 3.5 keV energy range provided unprecedented levels of sensitivity (hundreds of times 
better than had previously been achieved). The Einstein detected X-ray emission from 
a number of neutron stars (mainly as compact sources in supernova remnants) such as 
the middle-aged radio pulsars B0656+14, B1055-52 and the old pulsar B0950+08. The 
Einstein observatory re-entered the Earth's atmosphere and burned up on 25 March 1982. 

"'^The term "isolated" is omitted hereafter in the text however aU X-ray observations presented in this 
thesis concern isolated neutron stars 



2 Chapter 1. X-ray emission from Radio Pulsars 

The next "decade of space science" was opened in the 1990s with the launch of the ROSAT 
mission that was sensitive in the 0.1 — 2.4 keV energy range. One of the major results 
achieved with the ROSAT was the identification of the 7-ray source Geminga as a pulsar, 



hence a neutron star (Halpern and Holt 1992). 



The current era of X-ray observations of neutron stars was begun with the launch 
of two satellites: the XMM-Newton owned by the European Space Agency and the 
Chandra owned by the National Aeronautics and Space Administration. These two 
grazing-incidence X-ray telescopes were placed in orbit in 1999. They were equipped 
with cameras and high-resolution spectrometers sensitive to low-energy X-rays: from 0.08 
to lOkeV for the Chandra and from 0.1 to 15keV for the XMM-Newton. While the two 
observatories have similar designs, they are not identical. The XMM-Newton observat- 
ory has three X-ray telescopes that provide six times the collecting area and a broader 
spectral range in images than the Chandra, while the Chandra has a much finer spatial 
resolution and a broader spectral range in its high-resolution spectroscopy than does the 
XMM-Newton. Both observatories are in a highly-elliptical orbit that permits continuous 
observations of up to 40 hours. The Chandra and XMM-Newton have greatly increased 
the quality and availability of observations of X-ray thermal radiation from neutron star 
surfaces. The total number of isolated neutron stars of different types detected in X-rays 
is hard to find since not all data have been published. Some authors estimate that about 



one hundred rotation-powered pulsars were detected in the X-rays (Zavlin, 2007a Becker 



2009). 



1.2 X-ray emission from isolated neutron stars 

X-ray emission is a common feature of all kinds of neutron stars. Furthermore, X-ray 
observations have led to the discovery of other types of neutron stars that for various 
reasons were missed in the standard searches for radio pulsars. These new classes, such 
as X-ray Dim Isolated Neutron Stars, Central Compact Objects in supernovae remnants. 
Anomalous X-ray Pulsars, and Soft Gamma-ray Repeaters, are only a small fraction of 
the whole number of observed pulsars but provide valuable information on the diversity 
of the neutron star population. 

X-ray radiation from an isolated neutron star can in general consist of two distinguish- 
able components: thermal and nonthermal emissions. The thermal emission can originate 
either from the entire surface of a cooling neutron star or from spots around the magnetic 
poles on the stellar surface (polar caps and adjacent areas). The temperature of a neut- 
ron star at the moment of its formation is extremely high - its value is even as high as 
10^° — 10^^ K. Such a high initial temperature leads to very fast cooling, and after several 
minutes the temperature of the star interior falls to 10^ — 10^° K. After 10 — 100 yr the 
neutron star will cool down to a few times 10^ K. At this point, depending on the still 
poorly known properties of super-dense matter, the temperature evolution can follow two 
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different scenarios. The standard cooling scenario predicts that the temperature decreases 
gradually, down to ~ (0.3 — 1) x 10^ K by the end of the neutrino cooling era and then falls 
exponentially to temperatures lower than ~ 10^ K in ~ 10^ yr. In the accelerated cooling 
scenario, which implies higher central densities (up to lO^^gcm"^) and/or exotic interior 
composition (e.g. quark plasma), at the age of ~ 10 — 100 yr the temperature decreases 
rapidly down to ~ (0.3 — 0.5) x 10® K and is followed by a more gradual decrease down to 



the same ~ 10 K in ~ 10 yr (Becker, 2009). The thermal evolution of neutron stars is 



very sensitive to the composition (and structure) of their interiors, therefore, measuring 
surface temperatures is an important tool in studying super-dense matter. In addition 
to a thermal component emitted from the entire surface, other thermal components can 
also be seen. One of these additional components could be related to the reheating of the 
polar cap region by relativistic backflowing particles (electron and/or positrons) created 
and accelerated in the so-called polar gaps (see Chapter [3]). The temperature of these hot 
spots does not obey the same age dependence as the thermal evolution of neutron stars. 
Thus, depending on the pulsar age the thermal radiation may be dominated by either the 
entire surface (for younger neutron stars) or the hot spot components (for older neutron 
stars). The nonthermal component is usually attributed to the emission produced by Syn- 
chrotron Radiation (SR) and/or Inverse Compton Scattering (ICS) of charged relativistic 
particles accelerated in the pulsar magnetosphere. As the energy of these particles follows 
a power-law distribution, nonthermal emission is also characterised by power-law spectra. 
The X-ray spectrum of a neutron star (thermal and nonthermal) depends on many 
factors, e.g. the age of the star (r), inclination angle, strength and geometry of the 
magnetic field, etc. In most of the very young pulsars (r ~ 1 kyr) the nonthermal com- 
ponent dominates, thus making it impossible to accurately measure the thermal flux - 
only the upper limits for the surface temperature can be derived. As a pulsar becomes 
older, its activity (nonthermal luminosity) decreases roughly proportionally to its spin- 
down luminosity Lsd- A spin-down luminosity generally decreases with the increasing 
star age, as Lsd oc r""^, where m ^ 2 — 4 depends on the pulsar dipole breaking index 



(Zavlin, 2007b). With the increase of the pulsar age the luminosity of the surface thermal 
radiation decreases more slowly than the luminosity of the nonthermal one. Thus, the 
thermal radiation from an entire stellar surface can dominate a soft X-ray spectrum of 
middle-aged (r ~ 100 kyr) and some younger (r ~ 10 kyr) pulsars. For the old neutron 
stars (r > IMyr), a surface temperature Tg < 0.1 MK makes it impossible to detect the 
thermal radiation from the entire surface by available observatories. However, most of the 
pulsar models predict the heating up of polar caps to very high temperatures (Tg > 1 MK) 
by relativistic particles which are created in the pulsar acceleration zones. Conventionally, 
it is assumed that the polar cap radius is R^p = y2TiE?/cP. 

Since the spin-down luminosity Lsd is the source for both nonthermal (magneto- 
spheric) and thermal (polar cap) components, it is hard to predict which one would 
prevail in the X-ray flux of old neutron stars. Figure |1.1^ shows the ratio of a thermal 
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luminosity to a nonthermal one as a function of the pulsar age. Since calculating this 
ratio is possible only for pulsars with blackbody plus power-law fit, only these pulsars 
are included in the Figure. There is also a significant number of pulsars (16) with the 
spectra dominated by nonthermal components. Let us note that it is impossible to de- 
termine the thermal components for these pulsars. Most of them are young neutron stars 
~ 10^ — 10^ yr, but there are also much older ones (~ 10^ yr). In addition, there is a 
group of 4 pulsars with the spectra dominated by thermal components (without a visible 
nonthermal component). Their age also varies in quite a wide range 10'^ — 10^ yr. 
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Figure 1.1: Ratio of X-ray luminosities (thermal and nonthermal components) as a function of 
r (panel a) and B^ (panel b). The plots contain only those pulsars for which the BB+PL (Black- 
Body plus Power-Law) spectral fit exists. The number labels at the points correspond to the pulsar 
numbers in Table fmi 



As it follows from the left panel of Figure |1.1[ there is no obvious relation between 
pulsar age and the ratio of luminosities. The spectra of pulsars with a similar age may be 
dominated either by nonthermal (e.g. PSR B1951-I-32, PSR B1046-58) or thermal (e.g. 
PSR B0656-I-14, PSR J0538-I-2817) components. It is difficult to provide a more detailed 
analysis because, on the one hand, the observational errors are large and, on the other 
hand, a separation of thermal and nonthermal components is often not possible. The 
ratio of luminosities also does not show any correlation with the strength of the dipolar 



magnetic field (see the right panel of Figure 1.1 ). Let us note that the value of the dipolar 



magnetic field is conventionally calculated by adopting that the spin-down luminosity is 
equal to the power of magneto-dipole radiation (neglecting the influence of a pulsar wind). 
Then, assuming a dipolar structure of the neutron star magnetic field down to the stellar 
surface, we estimate its strength (measured in Gauss) at the pole as 



B. 



2.02 X 10'^ PP^ 



-15 



0.5 



:i-i) 



Here P is a period in seconds and P-15 = P x 10 . The actual strength of the surface 
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magnetic field can greatly exceed the above value (see Chapter |2 



Table 1.1 presents the basic parameters of the 48 pulsars that we use in this thesis, 



while the results of the X-ray observations of these pulsars are listed in Tables 1.2 and 
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Table 1.1: Parameters of rotation powered normal pulsars with detected X-ray radiation. The 
individual columns are as follows: (1) Pulsar name, (2) Barycentric period P of the pulsar, (3) Time 
derivative of barycentric period P, (4) Canonical value of the dipolar magnetic field B^ at the poles, 
(5) Spin-down energy loss rate Lsd (spin-down luminosity) , (6) Dispersion measure DM, (7) Best 
estimate of pulsar distance D (used in all calculations), (8) Best estimate of pulsar age or spin-down 
age T = Pj (2P), (9) Pulsar number (used in the Figures). Parameters of the radio pulsar have 
been taken from the ATNF catalogue. 



Name 


P 


P 


Bd 


log LsB 


DM 


D 


T 


No. 




(s) 


{w-^^') 


(1012 G) 


(crgs-l) 


(cm^3 pc) 


(kpc) 






J0108-1431 


0.808 


0.077 


0.504 


30.76 


2.38 


0.18 


166 Myr 


1 


J0205+6449 


0.066 


193.9 


7.210 


37.43 


141 


3.20 


5.37 kyr 


2 


B0355+54 


0.156 


4.397 


1.675 


34.65 


57.1 


1.04 


564 kyr 


3 


B0531+21 


0.033 


422.8 


7.555 


38.66 


56.8 


2.00 


1.24 kyr 


4 


J0537-6910 


0.016 


51.78 


1.846 


38.69 


- 


47.0 


4.93 kyr 


5 


J0538+2817 


0.143 


3.669 


1.464 


34.69 


39.6 


1.20 


30.0 kyr 


6 


B0540-69 


0.050 


478.9 


9.934 


38.18 


146 


55.0 


1.67 kyr 


i 


B0628-28 


1.244 


7.123 


6.014 


32.18 


34.5 


1.45 


2.77 Myr 


8 


J0633+1746 


0.237 


10.97 


3.258 


34.51 


- 


0.16 


342 kyr 


9 


B0656+14 


0.385 


55.00 


9.294 


34.58 


14.0 


0.29 


111 kyr 


10 


J0821-4300 


0.113 


1.200 


0.743 


34.52 


— 


2.20 


3.7 kyr 


11 


B0823+26 


0.531 


1.709 


1.924 


32.65 


19.5 


0.34 


4.92 Myr 


12 


B0833-45 


0.089 


125.0 


6.750 


36.84 


68.0 


0.21 


11.3 kyr 


13 


B0834+06 


1.274 


6.799 


5.945 


32.11 


12.9 


0.64 


2.97 Myr 


14 


B0943+10 


1.098 


3.493 


3.956 


32.00 


15.4 


0.63 


4.98 Myr 


15 


B0950+08 


0.253 


0.230 


0.487 


32.75 


2.96 


0.26 


17.5 Myr 


16 


B1046-58 


0.124 


96.32 


6.972 


36.30 


129 


2.70 


20.3 kyr 


17 


B1055-52 


0.197 


5.833 


2.166 


34.48 


30.1 


0.75 


535 kyr 


18 


J1105-6107 


0.063 


15.83 


2.020 


36.40 


271 


7.00 


63.3 kyr 


19 


J1119-6127 


0.408 


4022 


81.80 


36.36 


707 


8.40 


1.61 kyr 


20 
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Table 1.1 
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Name 


P 


P 


5d 


log LsD 


DM 


D 


T 


No. 




(s) 


(10-15) 


(1012 G) 


(ergs-i) 


(cm-3 pc) 


(kpc) 






J1124-5916 


0.135 


747.1 


20.31 


37.08 


330 


6.00 


2.87 kyr 


21 


B1133+16 


1.188 


3.734 


4.254 


31.94 


4.86 


0.36 


5.04 Myr 


22 


J1210-5226 


0.424 


0.066 


0.338 


31.53 


- 


2.45 


102 Myr 


23 


B1259-63 


0.048 


2.276 


0.666 


35.91 


147 


2.00 


332 kyr 


24 


J1357-6429 


0.166 


360.2 


15.62 


36.49 


128 


2.50 


7.31 kyr 


25 


J 1420-6048 


0.068 


83.17 


4.810 


37.00 


360 


8.00 


13.0 kyr 


26 


B1451-68 


0.263 


0.098 


0.325 


32.32 


8.60 


0.48 


42.5 Myr 


27 


J1509-5850 


0.089 


9.170 


1.824 


35.71 


138 


2.56 


154 kyr 


28 


B1509-58 


0.151 


1537 


30.73 


37.26 


252 


4.18 


1.55 kyr 


29 


J1617-5055 


0.069 


135.1 


6.183 


37.20 


467 


6.50 


8.13 kyr 


30 


B1706-44 


0.102 


92.98 


6.235 


36.53 


75.7 


2.50 


17.5 kyr 


31 


B1719-37 


0.236 


10.85 


3.234 


34.52 


99.5 


1.84 


345 kyr 


32 


J1747-2958 


0.099 


61.32 


4.972 


36.40 


102 


5.00 


25.5 kyr 


33 


B1757-24 


0.125 


127.9 


8.075 


36.41 


289 


5.00 


15.5 kyr 


34 


B1800-21 


0.134 


134.1 


8.551 


36.34 


234 


4.00 


15.8 kyr 


35 


J1809-1917 


0.083 


25.54 


2.936 


36.26 


197 


3.50 


51.3 kyr 


36 


J1811-1925 


0.065 


44.00 


3.407 


36.81 


- 


5.00 


23.3 kyr 


37 


B1823-13 


0.101 


75.06 


5.575 


36.45 


231 


4.00 


21.4 kyr 


38 


J1846-0258 


0.326 


7083 


97.02 


36.91 


- 


6.00 


0.73 kyr 


39 


B1853+01 


0.267 


208.4 


15.08 


35.63 


96.7 


2.60 


20.3 kyr 


40 


B1916+14 


1.181 


212.4 


31.99 


33.71 


27.2 


2.10 


88.1 kyr 


41 


J1930+1852 


0.137 


750.6 


20.47 


37.08 


308 


5.00 


2.89 kyr 


42 


B1929+10 


0.227 


1.157 


1.034 


33.59 


3.18 


0.36 


3.10 Myr 


43 


B1951+32 


0.040 


5.845 


0.971 


36.57 


45.0 


2.00 


107 kyr 


44 


J2021+3651 


0.104 


95.60 


6.361 


36.53 


371 


10.0 


17.2 kyr 


45 


J2043+2740 


0.096 


1.270 


0.706 


34.75 


21.0 


1.80 


1.20 Myr 


46 


B2224+65 


0.683 


9.659 


5.187 


33.08 


36.1 


2.00 


1.12 Myr 


47 


B2334+61 


0.495 


191.7 


19.69 


34.79 


58.4 


3.10 


40.9 kyr 


48 
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1.3 Nonthermal X-ray radiation 

The nonthermal emission, which is generally observed from radio to 7-ray frequencies, 
should be generated by charged particles accelerated at the expense of rotational energy 
in the magnetosphere of the neutron star. Nonthermal X-ray radiation is characterised 
by highly anisotropic emission patterns, which give rise to large pulsed fractions. The 
pulse profiles often show narrow (often double) peaks, however, in many cases nearly 
sinusoidal profiles are observed. As the X-ray efficiency is strongly correlated with Lsd, 
the most X-ray luminous sources (among rotationally powered pulsars) are the Crab pulsar 
and two young pulsars in the Large Magellanic Cloud, which are the only pulsars with 



Lsd > lO^^ergs-i (Mereghetti, 2011). 



Becker and Truemper (1997) suggested that in the 0. 1-2.4 keV band ROSAT sources 
that are identified as rotation-powered pulsars exhibit an X-ray efficiency which can be ap- 
proximated as a linear function Lx = ^Lsb, where the total X-ray efficiency 
^ = ^BB + Ct ~ 10~^, here ^^^ and ^^^ are efficiencies of the thermal (without the 
cooling component) and nonthermal X-ray emission, respectively. The higher sensitivity 
of both the Chandra and XMM-Newton allows detection of less efficient (^ < 10^^) X-ray 



pulsars (see Figure [L2| ) . Becker (2009) suggested that for these faint pulsars the orient- 



ation of the magnetic/rotation axes to the observer's line of sight might not be optimal. 
We believe that the efficiency of spin-down energy conversion processes is mostly affected 
by the strength and structure of the surface magnetic field. The variation of C, is rather 
due to the nature of physical processes than the geometrical effects. Let us note that 



the nonthermal X-ray luminosities presented in Figure |1.2| are calculated assuming an 
isotropic radiation pattern. In general, the X-ray emission pattern differs quite essentially 
from the isotropic one. Thus, one should introduce a beaming factor as the ratio of the 
opening angle of the radiation cone to the full solid angle Air. Since a beaming factor is 
generally unknown, the actual X-ray efficiency may differ by up to an order of magnitude 
(or even more) than we have presented. 

Various fitting parameters and efficiencies of nonthermal X-ray radiation suggest that 
the efficiency of processes responsible for the generation of nonthermal X-ray radiation 



should highly depend on the pulsar parameters (see Figure 1.2). The fitting parameters 
for the data of all pulsars show a linear trend with C, ~ 10~^, however, if we divide them 
into two groups of less and more luminous pulsars, we can see that the fitting parameters 
for these two groups differ from one another. The efficiency of less luminous X-ray pulsars 
depends on Lsd to a lesser extent than is the case for more luminous pulsars. 

As we mentioned in the Introduction, there are two main types of models: the polar 
cap models and the outer gap models. The outer gap model was proposed to explain the 



bright 7-ray emission from the Crab and Vela pulsars (Cheng et al. , 1986a|[b). Placing a 



7-ray emission zone at the light cylinder, where the magnetic field strength is considerably 
reduced to -Blc = B^ (_R/_Rlc) > provides higher 7-ray emissivities that are in somewhat 
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Figure 1.2: Nonthermal luminosity within the 0.1 — lOkeV band (Lnt) vs spin-down luminosity 
(Lsd)- The black solid line corresponds to the linear fitting for all pulsars, while the blue dotted 
and red dashed lines correspond to the linear fit for less luminous (Lsd < lO^^ergs^^) and more 
luminous (Lgo > lO'^^ergs^^) pulsars, respectively. 

better agreement with the observations. The observational data can be interpreted with 
any of the two models, although under completely different assumptions about pulsar 
parameters. 



1.3.1 Observations 

Generally, the X-ray spectrum of relatively young (r < lOkyr) and middle-aged 
(r < lOkyr) pulsars is dominated by the nonthermal component. However, it is not 
possible to find an exact correlation between r and the type of spectra, i.e. which com- 



ponent, thermal or nonthermal, dominates the spectrum (see the left panel of Figure 1.1). 
As we mentioned above, it is quite often impossible to resolve the components. The Crab 
pulsar (r = 958 yrs) is the most characteristic example of a young pulsar. The upper limit 
for X-ray luminosity of the Crab pulsar (one of the strongest known X-ray radio pulsars) 
is about L^T = 8.9 x 10^^ erg s~^. This value is calculated assuming an isotropic radiation 
pattern, however, even if we assume an angular anisotropy of the radiation (beaming 
factor ?^ 1/47?), the lower limit of its luminosity L^!^ = 7.1 x lO^^ergs"^ continues to be 
very high. The luminosities calculated above correspond to the following X-ray efficien- 
cies: ^ = 10"^'''^ (isotropic radiation pattern) and ^™'" = lO"^'*^^ (anisotropic radiation 
pattern). Although ^^^ is quite small, the nonthermal component still obscures all the 
thermal ones. To obtain a similar efficiency of the thermal radiation from the entire stellar 
surface, its temperature should be Tg = 5.9 x 10^ K (assuming R = 10km), which vastly 
exceeds the upper limit (Tg < 2.3 x 10^ K). Furthermore, the temperature of the polar 
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caps should be about 2.5 x 10^ K to obtain a comparable luminosity. 

The Vela-like pulsars compose another characteristic group of pulsars. This group 
consists of pulsars with high spin-down luminosities but considerably low X-ray efficiencies 
^NT ~ 10^^. A characteristic age of the Vela is about 1.1 x lO^yrs (10 times older than 
the Crab), but it can still be classified as a very young pulsar. The nonthermal luminosity 
of a Vela pulsar is Lj^!^ = 4.2 x 10^^ and efficiency ^™'"' = 10~^'^^. Some of the Vela- 
like pulsars (like the Vela itself) also exhibit a thermal component, which in some cases 
can be comparable to the nonthermal component. The thermal efficiency of the Vela 
^gg = 10~^'^^ is quite similar to ^™'''' = 10~^'^^, but if we assume an anisotropic radiation 
pattern of the nonthermal component than ^""" = 10"^'^^, thus even less than ^^g. 

The third group includes pulsars with low spin-down luminosity Lsd ^ 10^^. In 
most cases, the X-ray spectra of such pulsars (e.g. PSR 9050+08, PSR B1929+10) have 
both thermal and nonthermal components, with similar efficiencies. Thus, the spectrum 
fitting procedure is more complicated. The nonthermal X-ray efficiencies of these pulsars, 
^^^ ~ 10^'^, are considerably higher than those of the Vela-like pulsars. Note that even 
when the observed spectra are dominated by nonthermal radiation, we cannot rule out 
a situation that the thermal component is stronger than the nonthermal one, but due to 
unfavourable geometry we cannot observe it. 

Even with the improved quality of X-ray observations performed by both the Chandra 
and XMM-Newton, the available data do not allow us to fully discriminate between the dif- 
ferent emission scenarios. However, these data can be used to verify whether the proposed 



model of X-ray emission meets all the requirements. Table [T72| presents the observed spec- 
tral properties of pulsars showing nonthermal components. 



Table 1.2: Observed spectral properties of rotation-powered pulsars with X-ray spectrum showing the nonthermal (power-law) component. The individual 
columns are as follows: (1) Pulsar name, (2) Additional information, (3) Spectral components required to fit the observed spectra, PL: power law, BB: 
blackbody, (4) Pulse phase average photon index, (5) Maximum nonthermal luminosity Lnt. (6) Maximum nonthermal X-ray efficiency C™. (7) Minimum 
nonthermal X-ray efficiency ^^^ , (8) Total thermal luminosity Lrb. (9) Thermal efficiency ^gg, (10) References, (11) Number of the pulsar. Both nonthermal 
luminosities and efficiencies were calculated in the 0.1 — lOkeV band. The maximum value was calculated with the assumption that the X-ray radiation is 
isotropic while the minimum value was calculated assuming strong angular anisotropy of the radiation {S,^^ « 1/ (4vr)-^j^^ ). Pulsars are sorted by nonthermal 
X-ray luminosity (5). 



Name 



Comment 



Spectrum 



Photon-Index 



loeL 



NT 



(ergs 1) 



lOff f log f 



logLBB log 4 

(crgs^i) 



Ref. 



No. 



CO 



o 

a 
i-t 

B 
X 

r 

i-i 

Bt 

rf- 

5' 
3 



B0540-69 

B0531+21 

J0537-6910 

B1509-58 

J1846-0258 

J1420-6048 

J202H-3651 

J1617-5055 

J 1747-2958 
J1811-1925 

J1930+1852 



N158A, LMC 

Crab 

N157B, LMC 

Crab-like pulsar 

Kes75 



Crab-like pulsar 

Mouse 

Gil. 2-0.3 

Crab-like pulsar 



PL 
PL 
PL 
PL 
BB + PL 

PL 
PL, BB 

PL 
PL, BB 

PL 

PL 



0.11 

-0.11 

0.07 

0.07 



i.8ot^:}° 



1.92 
1.63 



0.04 

-0.04 

0.10 

0.10 



1.19 
1.90 

1.60 
1.70 



1 14+0.06 

1 on+0.08 

-'-•oU_o.08 
-0.39 
-0.32 



0.40 

-0.40 

0.30 

0.20 



0.97: 
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-1.27 


-2.37 


35.95 


-2.71 


-3.81 


35.95 


-2.74 


-3.84 


35.24 


-2.00 


-3.10 


35.13 


-1.78 


-2.88 


34.77 


-2.25 


-3.35 


34.36 


-2.17 


-3.27 


34.25 


-2.95 


-4.05 


34.09 


-2.31 


-3.41 


33.97 


-2.84 


-3.94 



KaOl 



CuOl 



34.06 



33.78 



-2.85 



-2.75 



Ca08 



Be09 



Mi05 



De06 



Ng08 He03 



Be09 





RoOl 




Va08He04 


Ka09 


Be02 




Ga04 




Ro03 


Ro04 



1 20+°-20 

-'-•^'J-0.20 



33.92 



-3.15 



-4.25 



Lu07 



Ca02 



7 

4 

5 

29 
39 

26 
45 
30 
33 
37 

42 
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Name 


Comment 


Spectrum 


Photon-Index 


log Lnt 

(erg s-l) 






log Lbb 

(ergs-l) 


log^BB 


Ref. 


No. 


J1105-6107 

B1757-24 

B1951+32 

J0205+6449 

J1119-6127 

J1124-5916 

B1259-63 

B0833-45 

B1706-44 

J1357-6429 

B1853+01 
B 1046-58 
B1916+14 
J1509-5850 
B1823-13 

B1800-21 


Duck 

CTB 80 

3C58 

G292.2-0.5 

Vela-like pulsar 
Be-star bin 

Vela 
G343.1-02.3 

W44 
Vela-like pulsar 

MSH 15-52 
Vela-like 

Vela-like pulsar 


PL 

PL 
BB + PL 
BB + PL 
BB + PL 

PL 

PL 
BB + PL 
BB + PL 
BB + PL 

PL 

PL 

BB, PL 

PL 
BB + PL 

PL + BB 


i-OO-O.OS 
-r 70+0.02 
-'-•'^^-0.04 

1 ciO+0-30 

-'-•""-0.10 
i-OJ_0.04 

r, 7^^ + 0.40 

^•'"-0.40 

9 00+0-50 

^•""-0.50 
1 qO+0-20 

-r 00+0.48 
^-•^^-0.48 
-, 7Q+0.40 
-■-•'"-0.20 

cn+i.eo 

O-0U_0.70 

-'-•""-0.30 
-, 7Q+0.70 
-■-•'"-0.70 

1 40+0-60 

-'-•^"-0.60 


33.91 
33.46 
33.22 
33.10 
32.95 

32.91 
32.87 
32.62 
32.16 
32.15 

32.07 
32.04 
32.00 
31.80 
31.78 

31.60 


-2.48 
-2.95 
-3.35 
-4.33 
-3.42 

-4.17 
-3.05 
-4.22 
-4.37 
-4.35 

-3.57 
-4.26 
-1.71 
-3.92 
-4.67 

-4.74 


-3.58 
-4.05 
-4.45 
-5.43 
-4.51 

-5.27 
-4.15 
-5.32 

-5.47 
-5.44 

-4.66 
-5.36 

-2.81 
-5.02 

-5.77 

-5.84 


31.95 
33.60 
33.37 

32.12 
32.78 
32.50 

31.07 
32.19 


-4.62 
-3.83 
-3.00 

-4.72 
-3.76 
-3.99 

-2.63 

-4.27 
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Name 


Comment 


Spectrum 


Photon-Index 


log Lnt 

(ergs-i) 


, ^max 


1 ^min 


log Lbb 

(ergs-i) 


log^BB 


Ref. 


No. 


J1809-1917 

B2334+61 

J2043+2740 

B2224+65 

B0355+54 

B1055-52 

B0656+14 

J0633+1746 

B1929+10 

B0628-28 

B0950+08 

B1451-68 

B1133+16 

B0823+26 

B0943+10 

B0834+06 
J0108-1431 


Guitar 
Geminga 

Chameleon 


BB + PL 
BB + PL 
BB + PL 
PL, BB 
BB + PL 

BB+BB+PL 

BB+BB+PL 

BB+BB+PL 

BB + PL 

BB + PL 

BB + PL 
BB + PL 

BB, PL 
PL 

BB, PL 

BB + PL 
BB + PL 


1 90+0. 62 
-'-■^'-'-0.62 

2.80;^;°° 
2.20;°i° 
1.00^:^° 

1 70+010 

9 in+O.SO 
^•-'-'J-0.30 

-l--D8_o.o6 

1 7o-f0.46 
-'-•''J-0.66 

9 QO+0.91 
2-98_o,65 

1 01 +0.14 
-'-•'-'-'--0.14 

9 p; -1+0.36 
^•^-'--0.33 
1 c:o+0.43 
-'-■^°-0.33 

9 «n+o-7o 

1 n+0.50 
"J- -'-'-'-0.20 


31.57 
31.55 
31.41 
31.21 
30.92 

30.91 
30.26 
30.24 
30.23 
30.22 

29.99 
29.77 
29.52 
29.42 
29.38 

28.70 

28.57 


-4.68 
-3.24 
-3.34 
-1.87 
-3.73 

-3.57 
-4.33 
-4.27 
-3.36 
-1.94 

-2.76 
-2.56 
-2.42 
-3.23 
-2.64 

-3.41 
-2.19 


-5.78 
-4.34 
-4.44 
-2.97 
-4.83 

-4.67 
-5.42 
-5.37 
-4.46 
-3.04 

-3.86 
-3.66 
-3.52 
-4.33 

-3.74 

-4.51 
-3.29 


31.69 
32.06 
30.77 
30.51 
30.40 

32.63 
32.77 
31.67 
30.06 
30.22 

28.92 
29.27 
28.56 

28.40 

28.70 
27.94 


-4.56 
-2.73 
-3.98 
-2.57 
-4.25 

-1.85 
-1.81 
-2.84 
-3.53 
-1.94 

-3.82 
-3.06 
-3.38 

-3.62 

-3.41 
-2.82 








36 

48 

46 

47 

3 

18 

10 

9 

43 

8 

16 
27 
22 
12 
15 

14 
1 




Ka07 




Mc06 


Be04 


Hul2 


Hu07b 


Mc07 S194 










De05 




De05 


Ja05 


Mi08 


Te05 


Be05 










Za04 




P0I2 


Ka06 


Be04 


Zh05 Ka06 










Gi08 




Po 


12 


'!Pa09l 



w 



o 
tr 

CD 
i-i 

3 
X 

i-i 



O 
3 



00 



14 Chapter 1. X-ray emission from Radio Pulsars 

1.4 Thermal X-ray radiation 

1.4.1 Modelling of thermal radiation from a neutron star 

Thermal X-ray emission seems to be quite a common feature of radio pulsars. The black- 
body fit to the observed thermal spectrum of a neutron star allows us to obtain the 
redshifted effective temperature T°° and redshifted total bolometric flux F°° (measured 
by a distant observer). To estimate the actual (unredshifted) parameters, one should take 
into account the gravitational redshift, g,. = a/1 — 2GM/Rc'^, determined by the neutron 
star mass M and radius R, here G is the gravitational constant. Then the actual effective 
temperature and actual total bolometric flux can be written as: 



— It-ioo 



5 



T = g;'T 



F = g;^F°°. 



:i.2) 



Knowing the distance to the neutron star, D, we can use the effective temperature and 
total bolometric flux to calculate the size of the radiating region. If we assume that the 
radiation is isotropic (same in all directions, e.g. radiation from the entire stellar surface) 



then the radius of the radiating sphere (star) can be calculated as (Zavlin, 2007a) 



RT = D^^;^,=g;'R^, (1.3) 

where a ~ 5.6704 x lO^^ergcm^^ s^^ K~'* is the Stefan-Boltzmann constant. 



Knowing that Lbb = 4:7iD F and using Equations 1.2 and 1.3, we can write that 



^BB 



9;'L^B- (1-4) 



The modelling of thermal radiation is more complicated if we assume that it comes 
from the hot spot on the stellar surface. One should take into account such factors as: 
time-averaged cosine of the angle between the magnetic axis and the line of sight (cosi), 
gravitational bending of light, as well as whether the radiation comes from two opposite 
poles of the star or from one hot spot only. In general, the observed luminosity of the hot 
spot can be written as: 

L- = A-aT-^ (1.5) 

where A^ = 7ri?^^ is the observed area of the radiating region. 

The observed area of the radiating spot is also influenced by the geometrical factor /. 
This geometrical factor depends on following angles: ( between the line of sight and the 
spin axis, and a between the spin and magnetic axes, as well as on g^. and whether the 
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radiation comes from the star's two opposite poles or from a single hot spot only: 






(1.6) 



Finally, the hot spot luminosity can be calculated as 



^hs 



-J. f-lTCO 

9t J ^hs ■ 



(1.7) 



The luminosity of a radiating sphere with radius R±_ can be calculated as 
Lsp = AA^crT"^ = 4:TxR\(jT'^. On the other hand, if we assume that the radiation ori- 
ginates only from one hot spot we can calculate the luminosity as Lhs = AhgO'T'^. If the 
hot spot size is small compared to the star radius (i?hs ^ R) then the area of the spot 
can be calculated as Ahs ~ txR\. Thus, we have to remember that the luminosity calcu- 
lated assuming a spherical source will be four times higher than the actual luminosity of 
a radiating hot spot Lhs = 1/4 • Lgp (see the next section for details). 



1.4.2 Thermal radiation of hot spots 
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Figure 1.3: Coordinate system co-rotating with a star. The system was chosen so that the z-axis 
is along fi (the angular velocity) and o lies in the x-z plane (fiducial plane, i.e. at longitude zero). 
Here, /i is a unit vector in the direction of the magnetic axis and a is the angle between ft and jl, 
/3 is the impact parameter. 
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Let us consider a neutron star with two antipodal hot spots associated with polar caps 
of a stellar magnetic field. For simplicity's sake we assume that the spot size is small 
compared to the star radius R. If the magnetic axis p, is inclined to the spin axis by an 
angle a < 90°, the spots periodically change their position and inclination with respect 
to a distant observer. To compute the radiation fiuxes from the primary (closer to the 
observer) as well as the antipodal spot, we need to know their inclinations: cosii = n ■ o 
and cos 22 = n ■ o = — cosii, where n and n = — n are normal vectors to spots surfaces, 
and o is the unit vector pointing toward the observer. In the calculations we use a 
coordinate system co-rotating with a star. The z-axis is along f2 (the angular velocity) 



and o lies in the x-z plane (see Figure 1.3). 



In the chosen coordinate system we can write that the spherical coordinates of vectors 
have the following components: 

n = (n, 0, 0) ; 

o = (1, a + /3, 0); (1.8) 

p. = (1, a, fit) . 

Here the impact parameter /3 represents the closest approach of the line of sight to 
the magnetic axis. Note that ft = n and n = —p,; thus, we can write the following 
components of Cartesian coordinates: 

o = (sin(a + /3), 0, cos(a + /3)); 
n = (sin a cos fit, sin fit sin a, cos a) . 

Finally, the inclination angle for both primary and antipodal hot spots can be calcu- 
lated as 



cosii = sina ■ cosfit ■ sin (a -|- /3) -|- cosa ■ cos (a -|- /3) ; 

cos ^2 = — cos ii = — sin a ■ cos fit ■ sin {a + (3) — cos a ■ cos {a + (3) . 

We can estimate the contributions of the primary and antipodal spots to the observed 
X-ray fiux by calculating the time-averaged cosine of the angle between the magnetic axis 
and the line of sight. Note that we should take into account only positive values of cosi 
since for larger angles {i > 90°) the spot is not visible (at least in this approximation. 



see Section 1.4.3 for more details). Thus, the contribution of the primary spot can be 



calculated as follows: 



(cosii) 



p 

COS (ll) dt if tanatan(a+/3) < "^ ^^ tan«tan(a+/3) > 1' 

C0s(2i)dt+ / COs{ti)dt if - 1< tanatan(a+/3) < ^^ 



+ 



:i.ii) 
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where integration limits are 



P P 

t± = — ± — arccos 

2 2tt 



tan a tan {a + j3) 



:i.i2) 



On the other hand, the contribution of the antipodal spot can be calculated as 



(cos 22) 




if 



tan o tan(a+/3) 



< — 1 or 



cos (22) dt if — 1 < 



tan cj tan{a+/9) 



tanotan(a+/9) 
< 1. 



>1, 



:i.i3) 



Depending on the orientation of fi, o and /i, the thermal radiation may originate from: 



'!) both the primary and antipodal hot spots (see Figure 1.4); (2) mainly the primary spot 



but with a small contribution from the antipodal spot (see Figure 1.5); (3) the primary 



spot only (see Figure 1.6) 
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Figure 1.4: Cosine of the hot spots' inclination angle as a function of the pulsar phase for 
PSR B0950+08. The following parameters were used: a = 105.46°, /3 = 21.1°. For this geo- 
metry the thermal radiation of both primary and antipodal spots has a significant influence on the 
observed thermal flux. 
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Figure 1.5: Cosine of the hot spots' inclination angle as a function of the pulsar phase for 
PSR B1929+10. The following parameters were used: a = 35.97, /3 = 25.55. For this geometry 
there is only a small contribution from the antipodal spot. 
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Figure 1.6: Cosine of the hot spots' inclination angle as a function of the pulsar phase for 
PSR B0943+10. The following parameters were used: a = 11.58°, /3 = -4.29°. For this geo- 
metry only the primary hot spot is visible. 
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1.4.3 Gravitational bending of light near stellar surface 

The radius of a neutron star is only a few times larger than the Schwarzschild radius. The 
approach presented in the previous section does not include the gravitational bending 
effect, which is very strong in neutron stars. A strong gravitational field just above the 
stellar surface causes the bending of light. A photon emitted near a neutron star surface 
at an angle 6 with respect to the radial direction escapes to infinity at a different angle 
6' > 6. As a consequence, even when the spot inclination angle to the line of sight is i > 90° 
we can still observe thermal radiation from this spot. For a Schwarzschild metric we can 
calculate an observed flux fraction from the primary /i = Fi/Fq and antipodal /2 = F2/FQ 
spots. Here Fq is the maximum possible flux that is observed when the primary spot is 



viewed face-on. The primary and antipodal fluxes are given by (Beloborodov, 2002) 



/i = cos(^)(l-|) + ^, 

/2 = -COs(z)(l-^)+^, 

here r^ = 2GM/(? is the Schwarzschild radius. The primary spot is visible when 
coszi > — Tg/ (i? — Tg) and the antipodal spot when cosi2 > —fgl {R — rg). Consequently, 
both spots are seen when —rg/{R — rg) < cosi < r^/ {R — Vg), and then the observed 
flux fraction is 

2r 

/mi„ = /l + /2 = ^. (1.15) 

Hence, the blackbody pulse of primary and antipodal spots must display a plateau 
whenever both spots are in sight. Depending on the geometry of a pulsar we can distin- 



guish four classes (Beloborodov, 2002). Class I: when the antipodal spot is never seen 



and the primary spot is visible all the time (see the bottom right panel of Figure 1.7). 
For such pulsars the blackbody pulse has a perfect sinusoidal shape. Class II: when the 
primary spot is seen all the time and the antipodal spot is also in the visible zone for some 



time (see panels a, b and c of Figure 1.7). For these pulsars the sinusoidal pulse shape 
is interrupted by the plateau. Class HI: the primary spot is not visible for a fraction of 
the period and during this time only the antipodal spot is seen. The primary sinusoidal 
profile of such pulsars is interrupted by the plateau, and the plateau is interrupted by a 
weaker sinusoidal subpulse from the antipodal spot. Class IV: both spots are seen at any 
time. The observed blackbody flux of such pulsars is constant. 

The gravitational bending of light can significantly increase the visibility of a pulsar 



i.e. the observed flux, compare Figures 1.4 and 1.8). For some specific geometry the 



gravitational effects can also drastically change primary to the antipodal flux ratio (com- 
pare Figures 1.9 and 1.5). Our calculations show that for canonical values M = 1.4 M© 



and R = 10 km the gravitational effect is quite strong and the observed flux fraction is 
in the range of 0.85 — 1, while the geometric approach results in the 0.43 — 1 range (see 
Table Ol). 
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Table 1.3: Viewing geometry of pulsars. The individual columns are as follows: (1) Pulsar name, 
(2) Inclination angle with respect to the rotation axis a, (3) Opening angle p, (4) Impact parameter 
/3, (5) Total flux correction factor (including gravitational bending of light) (/) , (6) Flux correction 
factor of the primary spot (/i), (7) Flux correction factor of the antipodal spot (/2), (8, 9, 10) 
Time-averaged cosine of the angle between the magnetic axis and the line of sight: (cosi) (the total 
value), (cosii)(the primary spot), (003^2) (the antipodal spot), (10) Number of the pulsar. The 
gravitational bending effect was calculated using M = 1.4 M© and R = 10km. 
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Figure 1.7: Comparison of the observed flux fraction for geometric effect only (blue dotted line) and 
for geometric effect with the inclusion of a gravitational bending of light (red solid line). Individual 
panels correspond to the following pulsars: (a) PSR B1133+16, (b) PSR B1929+10, (c) PSR 



B0834+06 (d) PSR B0943+10. Parameters used in the calculations are presented in Table 1.3 
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Figure 1.8: Observed flux fraction / as a function of the rotation phase for PSR B0950+08. 
The following parameters were used: a = 105.46°, /? = 21.1°, M = 1.4M0, R = 10km. The 
gravitational bending of light increases the flux ratio of the antipodal to primary spots almost two 
times (0.85/0.5 = 1.7) and also increases the antipodal to the primary flux ratio (~ 1.3). 
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Figure 1.9: Observed flux fraction / as a function of the rotation phase for PSR B1929+10. 
The following parameters were used: a = 35.97, /3 = 25.55°, M = 1.4 M©, R = 10km. The 
gravitational bending of light increases the observed flux fraction two times ((/)/(cosi) = 1.98) 
and also increases the flux ratio of the antipodal to primary spots almost seven times (~ 6.7). 
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1.4.4 Observations 

As we have shown in the previous sections, the blackbody fit to the X-ray observations 
allows us to directly obtain the surface temperature Tg. Using the distance to pulsar D 
and the luminosity of thermal emission Lbb we can estimate the area of spot Abb- In 
most cases, Abb differs from the conventional polar cap area Ajp ^ 6.2 x 10^P~^ m^. We 
use parameter h = Adp/Abb to describe the difference between A^p and Abb- 
Entire surface radiation and warm spot component (b<l) 
In most cases the observed spot area Abb is larger than the conventional polar cap area 



(see Table 1.4). We can distinguish two types of pulsars in this group, with 6^1 and 

b<l- 

The first type is associated with observations of a thermal emission from the entire 
stellar surface and can be used to test cooling models. Although the entire surface radi- 
ation is strongest for young pulsars (r < 10 kyr ), observation of this radiation is very 
difficult due to the strong nonthermal component. A common practice is to separately 
fit the nonthermal (PL) and thermal (BB) components. However, the temperature ob- 
tained in such a BB fit (without the PL component) is most likely overestimated (e.g. 



see PSR J2021+3651 in Table 1.4). The nonthermal luminosity of an aging neutron star 
decreases proportionally to its spin-down luminosity Lsd, which is thought to drop with 
the star age as Lsd oc r^™, where m ~ 2 — 4 depends on the pulsar dipole breaking index 



(Zavlin, 2007b). As a pulsar becomes older, its surface temperature decreases. Depend- 
ing on the model, a predicted temperature decrease in the early stages is gradual (the 
standard model) or rapid (the accelerated cooling scenario). For a number of middle-aged 
(r ~ 100 kyr) and some younger (r ~ 10 kyr) pulsars the thermal radiation from the en- 
tire stellar surface dominates the radiation at soft X-ray energies (e.g. PSR J0633+1746, 
PSR B1055-52, PSR J0821-4300, PSR B0656+14, PSR J0205+6449, PSR J2021+3651). 
However, the sample of pulsars is not sufficient to unambiguously identify the cooling 
scenario. 

The second type is associated with observations of the warm spot area that is larger 
than the conventional polar cap area but still significantly less than the area of the star 
(6 < 1). The age of pulsars in this group varies from very young (r ~ 1 kyr) to middle- 
aged (r ~ 100 kyr) neutron stars. There is one exception, namely PSR J1210-5226, which 
is very old (r = 105 Myr) and can still be classified as a pulsar with the large warm spot 
component. Note, however, that the age of this pulsar is estimated using a characteristic 
value and if the pulsar period at birth is comparable with the current period then the 
age is highly overestimated (see, e.g. PSR J0821-4300). Furthermore, the fit to the X-ray 
spectrum was performed using only one thermal component and assuming no nonthermal 
radiation (PL). We believe that in many cases the size of the warm spot component and 
its temperature are overestimated by neglecting other sources of X-ray radiation, i.e. the 
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nonthermal component and the hot spot radiation. The small number of observed X- 
ray photons in some cases prevents a full spectrum fit with all thermal and nonthermal 
components. Therefore, we need observations with better statistics so that the spectrum 
fit can be extended using more spectral components. The non-dipolar structure of the 
surface magnetic field may cause significant deviations from the spherical symmetry of the 
transport processes in the crust. The magnetic field slightly enhances heat transport along 



the magnetic lines, but strongly suppresses it in the perpendicular direction (Greenstein 



and Hartke, 1983). Hence, the non-isothermality of the crust strongly depends on the 



geometry of the magnetic field (Geppert et al. , 2004). The drastic difference of the crustal 



transport process causes significant differences in the surface temperature distribution 



(Page et al. , 2006). Thus, the non-dipolar structure of the surface magnetic field can 



explain the existence of large warm spot components for young and middle-aged pulsars. 



We also suggested a mechanism of heating the surface adjacent to the polar cap (Szary 



et al. , 2011). The model of such heating is also based on the assumption that the pulsar 



magnetic field near the stellar surface differs significantly from the pure dipole one. The 
calculations show that it is natural to obtain such a geometry of the magnetic field lines 



that allows pair creation in the closed field line region (see Figure 1.10). 
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Figure 1.10: Cartoon of the magnetic field lines in the polar cap region. Red lines are open field 
lines and green dashed lines correspond to the dipole field. The blue arrows show the direction of 
the curvature photon emission. 



The pairs move along the closed magnetic field lines and heat the surface beyond 
the polar cap on the opposite side of the star. In such a scenario the heating energy 
is generated in lAR, and hence the luminosity of such a warm spot is limited by the 



power of the outflowing particles (for more details see Section 5.1.2.3). In most cases 



the large size of the emitting area and its high temperature make it unlikely that the 
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warm spot is related to the particles accelerated in lAR and is rather connected with the 
non-isothermality of the crust (e.g. PSR J1210-5226, PSR J1119-6127). 



The hot spot component (b > 1) 



In many cases the observed hot spot area Abb is less than the conventional polar cap area 
{b > 1). The temperature of the emitting area of these pulsars is usually higher than the 
temperature of the emitting area of pulsars with a warm spot component (6 < 1). The 
hot spot component is a natural consequence of the non-dipolar structure of the surface 



magnetic field (see Figure 1.10). In order to define an actual polar cap we need to follow 



the open field lines from the light cylinder up to the stellar surface by taking into account 
the non-dipolar structure of the surface magnetic field (see Figure 1.10[ ), which can be 
estimated by the magnetic flux conservation law as b = A^p/A]^^ = B^/ B^. Thus, if 6 ^ 1 
then Bs'> Ba- 

In neutron stars with positively charged polar caps (f2 • B < 0), the outflow of iron 
ions depends on the surface temperature and the surface binding energy (the so-called 



cohesive energy) (Cheng and Ruderman 1980 Jones 1986 Abrahams and Shapiro 1991 



Gil et al. , 2003). The cohesive energy of condensed matter increases with magnetic fleld 



strength (Medin and Lai, 2007). If for a given strength of the surface magnetic fleld the 



temperature is below the so-called critical temperature T^it the ions can tightly bind to 



the condensed surface and a polar gap can form (see Chapter^ for details). Medin and Lai 



(2008) calculated the dependence of the critical temperature (for a vacuum gap formation) 



on the strength of the surface magnetic fleld. In Figure 1.11 we present the positions of 
pulsars with derived surface temperature Tg and hot spot area Abb on the B^ — T^ diagram, 
where i?s is estimated as 5s = bB^. The red line represents the dependence of the critical 
temperature Tent on B^. We can see that in most cases the pulsars' positions follow the 
Bs — Tent theoretical curve. Note that the Figure includes only pulsars with a visible hot 
spot component (old pulsars). For younger pulsars (with warm spot components) it is 
not possible to estimate the surface magnetic fleld. There are a few cases which do not 
coincide with the theoretical curve. We believe that they correspond to the observations 
of warm spot component but with the area of radiation smaller than the conventional 
polar cap area (e.g. due to reheating of the surface beyond the polar cap, see Section 



1.4.4). 



According to our model the actual surface temperature is almost equal to the critical 
value Ts ^ Tent, which leads to the formation of the Partially Screened Gap (PSG) above 



the polar caps of a neutron star (Gil et al. , 2003). The hot spot parameters derived from 



X-ray observations of isolated neutron stars are presented in Table 1.4 
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Figure 1.11: Diagram of the surface temperature {Tq = Tg/ (lO^K)) vs. the surface magnetic field 



(i3i4 = Bs/ (lO^^ G)). The red line represents the dependence of Tent on B14 according to Medin 
and Lai ( 2008D and the dashed lines correspond to uncertainties in the calculations. The diagram 
includes all pulsars with 6 > 1 with the exception of PSR J2043+2740, for which the blackbody fit 
was performed using a fixed radius (estimation of the surface magnetic field is not possible). Error 
bars correspond to la. 



Table 1.4: Spectral properties of rotation-powered pulsars with detected blackbody X-ray components. The individual columns are as follows: (1) Pulsar 
name, (2) Spectral components required to fit the observed spectra, PL: power law, BB: blackbody, (3) Radius of the spot obtained from the blackbody fit 
-Rbb. (4) Surface temperature T^, (5) Surface magnetic field strength Bs, (6) b = ^^pMbb = B^/B^, A^^ - conventional polar cap area, Abb - actual polar 
cap area, (7) Bolometric luminosity of blackbody component Lbb, (8) Bolometric efficiency ^^g, (9) Maximum nonthermal luminosity -^^T ' (■'■^) Maximum 
nonthermal X-ray efficiency i^^ , (11) Best estimate of pulsar age or spin down age, (12) References, (13) Number of the pulsar. Nonthermal luminosity and 
efficiency were calculated in the 0.1 — lOkeV band. The maximum value was calculated with the assumption that the X-ray nonthermal radiation is isotropic. 
Pulsars are sorted by h parameter (6). 
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Chapter 2 

Model of a non-dipolar surface 
magnetic field 



2.1 The magnetic field of neutron stars 

Generally, the properties of pulsar radio emission support the assumption that the mag- 



netic field of pulsars is purely dipolar at least in the radio emission region (Radhakrish 



nan and Cooke, 1969). However, radio emission is generated at altitudes Rem of more 



than several stellar radii (e.g. Kijak and Gil (1997), Kijak and Gil (1998), Krzeszowski 



et al. (2009) and references therein). Thus, radio observations do not provide information 
about the structure of the magnetic field at the surface of the neutron star. On the other 
hand, strong non-dipolar surface magnetic fields have long been thought to be a necessary 
condition for pulsar activities, e.g. the vacuum gap model proposed by [Ruderman and 



Sutherland (1975) implicitly assumes that the radius of curvature of field lines above the 
polar cap should be about 10^ cm in order to sustain pair production. This curvature is 
approximately 100 times higher than that expected from a global dipolar magnetic field. 
Furthermore, to explain radiation from the Crab Nebula, the Crab pulsar should provide 
quite a dense stellar wind, as such a high particle multiplicity is not possible in a purely 
dipolar magnetic field. 

There are several theoretical studies concerning the formation and evolution of the non- 
dipolar magnetic fields of neutron stars (e.g. Blandford et al.|1983 Krolik|1991 Ruderman 
T99TJ |Arons|[T993J |Chen and Ruderman|[T993| [Geppert and Urpin|[T994l |Mitra et al.|[T999 



Page et aL]|2006 ). According to Woltjer (1964), the magnetic field in neutron stars results 
from the fossil field of the progenitor stars which is amplified during the collapse and 
remains anchored in the superfiuid core of the neutron star. Several authors also noted 
that during the collapse (or shortly after) there is possible magnetic field generation in the 



external crust, for instance, by a mechanism like thermomagnetic instabilities (Blandford 



et al. , 1983 ). Urpin et al. ( 1986 ) also showed that it is possible to form small-scale magnetic 



field anomalies in the neutron star crust with a typical size of the order of 100 meters. 
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The soft X-ray observations of pulsars presented in Chapter [T] show non-uniform sur- 
face temperatures which can be attributed to small-scale magnetic anomalies in the crust. 
Further observational arguments in favour of the non-dipolar nature of the surface mag- 



netic field can be found in many articles (e.g. Bulik et al. 1992 Thompson and Duncan 
1995 Bulik et al. 1995; Page and Sarmiento 1996 Thompson and Duncan] 1996 Becker 
and Truempi?lpM7l [Cheng et al.1p^98| |Rudak and Dyki|[T999| [Cheng and Zhang|pl99 



Murakami et al. 1999 Tauris and Konar 2001 Maciesiak et al. 2012). 



2.2 Modelling of the surface magnetic field 

In order to model a surface magnetic field we used the scenario proposed by |Cil et at 



(2002). In this scenario the magnetic field at the neutron star's surface is non-dipolar 



in nature, which is due to superposition of the fossil field in the core and crustal field 
structures. To calculate the actual surface magnetic field described by superposition of 
the star-centred global dipole d and the crust-anchored dipole moment m, let us consider 



the general situation presented in Figure 2.1 




Figure 2.1: Superposition of the star-centred global magnetic dipole d and crust-anchored local 
dipole anomaly m located at rg = (r^ ~ i?, ^ = Or) and inclined to the z-axis by an angle 9m- The 
actual surface magnetic field at radius vector r = (r, 9) is Bg = Bj + Bm, where Bd = 2d/r^ , 
B^ = 2?7i/|rrs|^, r is the radius and 9 - \s the polar angle. R is the radius of the neutron star and 



L is the external crust thickness. Gil et al. (2002) 
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The actual surface magnetic field is a sum of the global magnetic dipole and crust- 
anchored local anomalies 

B, = Bd + B^ + ... (2.1) 

Using the star-centred spherical coordinates with the 2;-axis directed along the global 
magnetic dipole moment we obtain: 

„ f 2d cos 6 dsm.6 \ , , 

Bd=(^:^,— ^,0 , (2.2) 



3(r-rs)(m- (r-r^)) -m|r-rs|2 
i^m = , r^ • (2.3j 

I s| 

Here rg = (r^, 6r, (j)r), m = (m, 9m, (pm) and the spherical components of Bm are explicitly 
given in Equation |2.7[ 



The global magnetic moment can be written as 

d = ^5pi?^ (2.4) 

where Bp = 6.4 x 10^^ ( PP ) ^/^ G is the dipole component at the pole derived from pulsar 
spin-down energy loss. 

The crust-anchored local dipole moment is 

m = ^B^AR\ (2.5) 

where AR ~ 0.05-R < L and L ~ 10^ cm is the characteristic crust dimension (for 
R = 10^ cm). For these values a local anomaly can significantly influence the surface 
magnetic field {B^ > B^) if m/d > 10""^. 

The system of differential equations for a field line of the vector field B = [B^., Bg, B^) 
in spherical coordinates can be written as 



de _ Be 
dr rB, 



B (2.6) 



dr r sin(9)Br ' 

The solution of these equations, with the initial conditions 9q = 6{r = R) and 0o=0(r=R) 
determining a given field line at the stellar surface, describes the parametric equation of 
the magnetic field lines. The spherical components of Bm can be written in the following 
form 
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Here 



and 



5™ 



BT 



D2.5 
1 



(3Tr^ - 3Tr + Dnir) , 
(3Tr^ + Dnie) , 



D2.5 



D 



D = Tg + r — 2rsr (sin 6r sin 6' cos 



v) + cos 9r cos ( 



T = mr?" — [rrirr^ + merg) . 



(2.7) 



(2i 



(2.9) 



According to the geometry presented in Figure [231 the components of the radius vector 
of the origin of the crust-anchored local dipole anomaly can be written as 



r^ = Vs (sin 9r sin 9 cos (0 — 0^) + cos 9r cos ( 
Tq = Tg (sin 6'r cos 9 cos (0 — 0^) + cos 6*^ sin ( 
r^ — r^ sin 9^ sin (0 — 0^) • 



The components of the local dipole anomaly are 



rrir = m (sin 9.m sin 9 cos (0 — 0m) + cos ft^ cos 9) , 
mg = m (sin ft^ cos 6^ cos (0 — 0^) + cos 6'^ sin 9) , 
m^ = —m sin ^^ sin (0 — 0^) . 



(2.10) 



(2.11) 



Finally, we obtain the system of differential equations from Equation 2.6| by substitu- 

and 



tions Br = B^ + S;", Be = Bf + B^ and B^ = Bl + 5™ (Equations 



2.2 



2.7) 



d^ 



B'g + B^ 



dr r {Bf + 5™) 



Bl 



©1, 



d0 

dr ~ r{B^ + B - r™) sin 6^ 



= $1. 



(2.12) 
(2.13) 



2.3 Curvature of magnetic field lines 

As Curvature Radiation (CR) may play a decisive role in radiation processes, it is im- 
portant to calculate the curvature (or curvature radius) for each field line. The curvature 
Pc = l/3f? of field lines (where 3f? is the radius of curvature) is calculated as (Gil et al. 



2002) 



pc 



ds\ 
dr J 



d^r ds dr d^s 



dr^ dr dr dr^ 



(2.14) 



2.3. Curvature of magnetic field lines 33 

where 

^ = ^J[l+r^ej + r^^lsm\e)]. (2.15) 

Thus, the curvature can be written in the form 

p, = (s,r'{j', + .p, + jiy^\ (2.16) 



where 



Ji —X2S1 — X1S2, 

Jz =2281 — ZIS2, 

Xi = sin 6 cos (p + rOi cos cos (p — r$i sin 6 sin 0, 
Yi = sin 6 sin + r0i cos ^ sin — r$i sin 6' cos 0, 
Zi = cos 6* — rOi sin 6', 
X2 = (201 + r62) cos 6* cos - (2$i + r$2) sin 9 sin < 

r (e^ + $1) sin 6^ cos + 2r9i$i cos 6 sin 0, 
F2 = (261 + r62) cos 6' sin + (2$i + r$2) sin 9 cos < 

r (61 + $1) sin 6^ sin + 2r9i$i cos 9 cos 0, 
Z2 = — ©1 sin ^ — Oi sin 9 — r02 sin 9 — rO^ cos 9, 



(2.17) 



^I=^l + r2e2 + r2$2sin2^, 

^2 =5f ^ (te? + r^BiGs + r$2 sin^ ^ + r2$i$2 sin^ 9 + r^Bi*? sin 9 cos 0) , 
dGi 



62 

$2 



dr 
d$i 

dr 



2.3.1 Numerical calculation of the curvature 



Let us note that when evaluating Equation 2.6 it was assumed that sin^ 7^ 0. Thus $1 
and $2 are undefined for ^ = 0. Figure 2.2 presents the first and second derivative ($1, 



$2) of the 0-coordinate of a magnetic field line with respect to the r-coordinate for a 
magnetic field structure with 5^ 7^ 0. 



The singularity in Equation 2.6 may result in an overestimation of curvature of field 
lines that cross the ^ = plane for a complex structure of the surface magnetic field. To 
solve this problem, and in addition to the analytical approach, the numerical calculation 
of the curvature of magnetic field lines was implemented. 

Let us consider three consecutive points of the given magnetic field line A, B, C (see 
Figure Q. 
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Figure 2.2: Plot of the first and second derivative of the (/)-coordinate of the magnetic field line 
with respect to the r-coordinate vs. the distance from the stellar surface. Values were calculated 



using the approach described in Section 2.3 Panels (a) and (b) show the first and second derivative 
of the (/)-coordinate while panel (c) shows the 6'-coordinate of the magnetic field line. 



magnetic field line 



C — (Ci,C2, Cs 




ne of nodes 



Figure 2.3: For any given three points {A, B, C) we can always find a common plane. We use the 
following transformations to achieve this: (I) shift the origin of the system to point A (prime), (II) 
rotate the shifted system by an angle c^y around the y'-axis and by an angle (^x around the x"-axis 
(double prime). After these transformations the 2;"-axis will be aligned with normal vector N and all 
points will lie in the x"y"-p\ane of such a system of coordinates. 
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To calculate curvature (or radius of curvature) in point B we can use the following 
procedure: 

• simplify the 3-D problem to 2-D by finding a common plane for all three points 

— move the origin of the coordinate system to point A 

— rotate the coordinate system to align the z-axis with the normal vector to the 
common plane of all three points {A, B, C) 

• calculate the radius of the circle passing through all three points (in a 2-D coordinate 
system) 

3-D to 2-D transition 

To simplify the calculations we shift the origin of the coordinate system so that point A 
will be the origin of the new system: 



A' = (0,0,0); 

B' = (51-^,52-^2,53-^3); (2.18) 

C = (Ci-A, C2-A2, C3-A3). 

The unit normal vector to the plane enclosing all three points (A', 5', C) can be 
calculated as 

N = ^^ = (^i> ^2, iV3) , (2.19) 

|b X c| 

where b = {B[, B'^, B'^) and c = {C[, q, C'^). 

The next step is to rotate the shifted coordinate system to align the z'-axis with normal 
vector N. In the new system all three points will lie in the a;"?/"-plane. In our calculations 
we rotate the shifted system by an angle <,y around the y'-axis, Ry {<,y), and a rotation by 
an angle q^ around the x"-axis, R^ {<ix)- The final rotation matrix can be written as 



^yx -fLy \^y) J^x \Sx) 



COS <,y sm <,nc sm <,y sm <,y cos <^x 

COS ?a; — sin q^ 

— sin qy cos ^y sin i,^ cos ^^ cos ^^ 

The Euler angles for these rotations can be calculated as 



(2.20) 



•^x 



atan2 {N^, N3) Q 



^where atan2 {y, x) equals: (1) arctan {y/x) if a; > 0; (2) arctan (y/x) + tt ii y > and x < 0; (3) 
arctan (y/o;) — tt if y < and a; < 0; (4) 7r/2 if y > and x — 0; (5) — 7r/2 if y < and a; = 0; (6) is 
undefined if y = and a; = 0. This function is available in many programming languages. 
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arctan 
arctan 

k 2 



TV] 

N3 
N2 



if A/'a = and N^ 



(2.21^ 







Finally, we can write the components of all three points in our new (shifted and 
double-rotated) system of coordinates as follows 



A" = Ry,A' = (0, 0, 0) ; 
B = RyxB = {B-^, i?2, 0) ; 

o — JrlyxLy — {^i, ^2 , uj . 



(2.22) 



Circle passing through 3 points (Bourne, 2012) 



Finding the radius of the circle passing through three consecutive points of a given mag- 
netic field line {A = (0, 0), B = {Bi, B2), C = (Ci, C2)) is an exact method for finding 
the radius of curvature 5R and hence the curvature p = l/3fJ of this line. Note that 
for simplicity's sake we hereafter describe points without double prime notation but 
they refer to coordinates in the shifted and double-rotated system of coordinates (e.g. 
B = (B,, B,) = {Bl B'^). 

Slope rrii of the line joining Aio B and slope 1112 of the line joining B to C (see Figure 



2.4) are given by 



mi 



m2 



Ay 

"Ax 
Ay 

"Ax 



B2 

C2 — B2 
Ci-Bi 



(2.23) 



In general, the centre of the circle passing through our points is given by 






mima {A2 - C2) + ma {Ai - Bi) - nii {Bi - Ci] 



2 {7712 — rrii) 



nil 






Ai + Bi 



+ 



A2 + B2 



(2.24) 
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Figure 2.4: The radius of curvature K of the magnetic field line at a given point B can be 
calculated as the radius of the circle passing through this and the neighbouring two points {A, C). 
The system of coordinates was moved and double- rotated so that A is in its origin and all points lie 
in the x"y"-plane. The slopes of the lines joining A to B and B to C are described by Equation 
[2:231 

Since point A is in the centre of the coordinate system we can simphfy these formulas 
as follows 






Vc 



2 B^B2 — 2 B1B2C1 + B2C1 — B2C2 — (-Bi — B2) C2 



2 {B,C2 - B2C1) 



Bl~{B2-2x,)B, 
2Bo 



(2.25) 



Finally, we can calculate the radius of curvature simply by finding the distance between 
the centre of the circle and any of the points on the circle (we have chosen point A) 



^ 



-p = V^ 



Vc- 



(2.26) 



In this thesis we consider complex structures of the surface magnetic field, thus the 
numerical method presented above was used in all the calculations of curvature. The 
analytical approach may result in an overestimation of curvature for points with 6^0 



(see Figure 2.5). 
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Figure 2.5: Curvature of the magnetic field lines vs. the height above the stellar surface cal- 
culated using the analytical approach described in Section |2.3| (green lines) and the numerical 
approach presented above (red lines). Panel (a) corresponds to the magnetic field line which has 
no (j) component, while panel (b) corresponds to a more general scenario i.e. the nonzero com- 
ponent of the magnetic field line. As can be seen, the analytical approach is not valid for every 
case. This is caused by the undefined value of the 4> derivative for 6* = (see Equation [2.6| ). Here 
/9„6 = i/^6 = p/ (10-6 cm-1) and ^q = K/ (lO'^cm). 



2.4 Simulation results 



In this section we model the surface non-dipolar magnetic field structure for some pulsars. 
Note that we can estimate the size of the polar cap and the strength of the surface magnetic 



field only for pulsars with an observed hot spot (see Section 1.4.4). Here we present only 



pulsars listed in Table 3.1 



We use spherical coordinates (r, 9, (p) to describe the location and orientation of crust- 
anchored local anomalies. The parameters of anomalies are as follows: Va = (r„ 



is a radius vector which points to the location of the anomaly and nia = {rria, Oa ,0a) 



IS 



its dipole moment. The value of iria is measured in units of the global dipole moment d, 
i.e. the moment which corresponds to the pulsar's global magnetic field. In the Figures 



showing a possible non-dipolar structure (e.g. Figure 2.6, 2.9, 2.12) the dashed lines 



correspond to the dipolar configuration of the magnetic field lines, while the solid lines 
correspond to the actual magnetic field lines (taking into account the crust- anchored 
anomalies). Green and red lines represent the open magnetic field lines for dipolar and 
non-dipolar structures, respectively. 



2.4. Simulation results 
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2.4.1 PSR B0628-28 



Pulsar B0628-28, a bright radio pulsar, was discovered by Large et al. (1969) during a 
pulsar search at 408 MHz. The pulsar period P ~ 1.24 s and its first derivative P_i5 ~ 7.1 
result in a dipolar component of magnetic field Sd = 6 x 10^^ G and a characteristic age 
Tc ~ 2.8 Myr, which makes it a typical, old pulsar. The large distance to this pulsar 
D = 1.44 kpc (evaluated using the Galactic free electron density model of Cordes and 



Lazio 2002) makes it impossible to use the parallax method to determine the distance 



with better accuracy. 

PSR B0628-28 is one of the longest period pulsars among those detected in X-rays. 
The pulsar was first detected in the X-ray band by ROSAT and then later observed 
with both the Chandra and XMM-Newton. Observations with the Chandra revealed 
no pulsations, while the XMM-Newton observations revealed pulsations with a period 



consistent with the period of radio emission (Tepedelenhoglu and Ogelman 2005). The 



inconsistency of the observations is a reflection of the fact that the pulsar is detectable just 
at the threshold of sensitivity of both the observatories. The two-component spectral fit 
(BB+PL) shows that both the nonthermal and thermal components have a comparable 
luminosity (at least if we assume that the nonthermal radiation is isotropic, see Table 



1.2|). PSR B0628-28 is characterised by one of the largest X-ray efficiencies among the 

10-2. 



observed pulsars ^bb ~ ^j— ~ 
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Figure 2.6: Possible non-dipolar structure of the magnetic field lines of PSR B0628-28. 
The structure was obtained using two crust anchored anomalies located at: 
ri = {0.95R, 4°, 0°), r2 = {0.95R, 5°, 180°), with the dipole moments 
mi = (4.5 X IQ-^d, 5°, 0°), m2 = (4.5 x IQ-'^d, 170°, 180°) respectively (blue arrows). The 
influence of the anomalies is negligible at distances D > 2R, where B^/B^ ~ m/d = 4.5 x 10^^ 



(top panel). For more details on the polar cap region see Figure 2.7 
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Figure 2.7: Zoom of the polar cap region of PSR B0628-28. See Figure 2.6 for a description 
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Figure 2.8: Dependence of a curvature of the open magnetic field lines on the distance from the 
stellar surface for PSR B0628-28. The distance is in units of the stellar radius ze = z/R and the 
curvature of the magnetic field lines is p_6 = 1/3^6 = p/ (lO^^cm^^). 
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2.4.2 PSR J0633+1746 



Geininga was discovered in 1972 as a 7-ray source by Fichtel et al. (1975). The visual 



magnitude of the pulsar was estimated by Bignami et al. (1987) to be of the order of 



25.5 . The pulse modulation was discovered in X-rays (Halpern and Holt 1992), in 



7-rays, and at optical wavelengths (Shearer et al. , 1998). Geminga has been determined 
to be a relatively old (r = 342 kyr) radio-quiet pulsar with a period P = 237 ms. The 
distance to the pulsar D = 0.16 kpc, evaluated using the parallax method, makes it the 
closest pulsar with available X-ray data. 

The pulsar exhibits one of the weakest radio luminosities known and a cutoff at fre- 



quencies higher than about 100 MHz. The model presented by Gil et al. (1998) explains 



this weak radio emission with absorption by the magnetised relativistic plasma inside the 



light cylinder. As the exact model of radio emission is still unknown (see Section 5.2.3) 



it is difficult to verify if this weak radio emission is a result of absorption or the absence 
of coherent radio emission. 



The three-component fit to the X-ray spectrum (PL+BB+BB, see Table 1.4) reveals 



the hot spot component with a size that is considerably smaller than the conventional 
polar cap size {b ~ 23). The entire surface temperature Tg = 0.5 MK is consistent with 
the theoretical value predicted by the cooling model. 
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Figure 2.9: Possible non-dipolar structure of the magnetic field lines of PSR J0633-I-1746. 
The structure was obtained using two crust anchored anomalies located at: 
ri = {0.95R, 3°, 0°), r2 = {0.95R, 6°, 180°), with the dipole moments 
mi = (5.5 X lO^^d, 10°, 0°), m2 = (5.5 x lO^^fi, 160°, 0°) respectively (blue arrows). The 
influence of the anomalies is negligible at distances D > 3.1R, where B^/B^ ~ m/d = 5.5 x 10^^ 



(top panel). For more details on the polar cap region see Figure 2.10 
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Figure 2.10: Zoom of the polar cap region of PSR J0633-I-1746. See Figure 2.9 for a description 
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Figure 2.11: Dependence of a curvature of the open magnetic field lines on the distance from the 
stellar surface for PSR J0633-I-1746. The distance is in units of the stellar radius zq = z/R and the 
curvature of the magnetic field lines is p_6 = 1/% = p/ (lO^^cm^^). 
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2.4.3 PSR B0834+06 

The bright radio emission of PSR B0834+06 shows frequent nulls (nearly 9% of the 



pulses is absent, see Rankin and Wright, 2008) . With a relatively long rotational period 
P = 1.27s and P_i5 ~ 7.1 (Taylor et al. , 2000), its inferred physical properties, e.g. 
i?d = 3 X 10^^ G, are close to the average. The characteristic age Tc = 2.97 Myr implies that 
the pulsar should be categorised as an old pulsar. The distance to the pulsar, estimated 
a.s D = 0.64 kc, was derived from its dispersion measure using the Galactic free-electron 



density model of Cordes and Lazio, 2002 Weltevrede et al. (2006) suggest a drift of 



subpulses, but the estimated value of a subpulse separation is larger than the pulse width. 
Despite the fact that the geometry based on the carousel model could be fitted to the 
observations, there is no clear evidence for a drift of emission between the components of 



the pulsar (Rankin and Wright 2007). 



The pulsar was detected in X-ray by Gil et al. (2008) with a total of 70 counts from 



over 50 ks exposure time. Because of the low statistical quality of the X-ray data, it was 
not possible to constrain the absorbing column density Nh- The two-component spectral 
fit (BE -|- PL), as presented in this thesis, was performed using the assumption that both 
the thermal and nonthermal fluxes are of the same order. 
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Figure 2.12: Possible non-dipolar structure of the magnetic field lines of PSR B0834+06. 
The structure was obtained using two crust anchored anomalies located at: 
ri = {0.95R, 2°, 0°), r2 = {0.95R, 5°, 180°), with the dipole moments 
mi = (3 X 10-3(i, 15°, 0°), m2 = (3 x lO^^d, 150°, 0°) respectively (blue arrows). The 
influence of the anomalies is negligible at distances D > 3.2R, where B^/B^ ss m/d = 3 x 10^^ 



(top panel). For more details on the polar cap region see Figure 2.13 
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Figure 2.13: Zoom of the polar cap region of PSR B0834+06. See Figure 2.12 for a description 
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Figure 2.14: Dependence of a curvature of the open magnetic field lines on the distance from the 
stellar surface for PSR B0834+06. The distance is in units of stellar radius {zq = z/R) and the 
curvature of the magnetic field lines is p_6 = 1/% = p/ (lO^^cm^^). 
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2.4.4 PSR B0943+10 



Pulsar B0943+10 is a relatively old pulsar with a characteristic age of Tc = 4.98 Myr. The 
pulsar period P = 1.1 s and its first derivative P_i5 ~ 3.5 result in the dipolar component 
of a magnetic field B^ = 4.0 x 10^^ G. Using the Galactic free-electron density model of 
Cordes and Lazio, 2002, we can estimate the distance to the pulsar D = 0.63 kpc. 



PSR B0943-I-10 is a well-known example of a pulsar exhibiting both the mode changing 
and subpulse drifting phenomenon. Strong, regular subpulse drifting is observed only in 
radio-bright mode, and only hints of the modulation feature have been found in the radio- 



quiescent mode. Very recent results presented by Hermsen et al. (2013) show synchronous 



switching in the radio and X-ray emission properties. When the pulsar is in a radio-bright 
mode, the X-rays exhibit only an unpulsed component. On the other hand, when the 
pulsar is in a radio-quiet mode, the fiux of X-rays is doubled and a pulsed component is 
also visible. 
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Figure 2.15: Possible non-dipolar structure of the magnetic field lines of PSR B0943+10. 
The structure was obtained using two crust anchored anomalies located at: 
ri = (0.96i?, 0°, 0°), r2 = (0.96i?, 15°, 0°), with the dipole moments 
mi = (2.0 X IQ-'^d, 180°, 0°), m2 = (6 x lO^^d, 20°, 180°) respectively (blue arrows). The 
influence of the anomalies is negligible at distances D > A.5R, where B^/B^ ss m/d = 2 x 10^^ 



(top panel). For more details on the polar cap region see Figure 2.16 



46 



Chapter 2. Model of a non-dipolar surface magnetic field 



-100 



-150 



-200 



-250 



-300 




2250 2300 2350 2400 2450 2500 2550 2600 2650 

X [m] 



Figure 2.16: Zoom of the polar cap region of PSR B0943+10. See Figure 2.15 for a description 
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Figure 2.17: Dependence of a curvature of the open magnetic field lines on the distance from the 
stellar surface for PSR B0943+10. The distance is in units of stellar radius {z% = z/R) and the 
curvature of the magnetic field lines is p_6 = 1/% = p/ (lO^^cm^^). 
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2.4.5 PSR B0950+08 

Pulsar B0950+08 is one of the strongest pulsed radio sources in the metre wavelength 
range. The pulsar radiation also exhibits an interpulse located at 152° from the main 
pulse (Smirnova and Shabanova, 1992). Based on the period P = 1.1s and its first 
derivative P_i5 ~ 3.5, we can estimate the pulsar's characteristic age Tc = 17.5 Myr. PSR 
B0950+08 has a relatively weak dipolar component of magnetic field B^ = 0.5 x 10^^. For 
this pulsar the distance D = 0.26 kpc was estimated using the parallax method. 



PSR B0950+08 was detected in the ultraviolet-optical range (2400 - 4600 ) by Pavlov 



et al. (1996) with the Hubble Space Telescope. Further observations suggest that the 



optical radiation of the pulsar is most likely of a nonthermal origin ( Mignani et al. , 2002 



Zharikovet al. 2004). 



X-rays from PSR B0950+08 were first detected with the ROSAT by [Manning and 
Willmore (1994[) (~ 55 source counts). Further X-ray observations revealed pulsations 



of the X-ray flux at the radio period of the pulsar (Zavlin and Pavlov 2004). The X- 



ray spectrum manifests two components (thermal and nonthermal). Which of the two 
components dominates the spectrum depends on the radiation pattern of the nonthermal 
component (isotropic or anisotropic). Due to the poor quality of the X-ray data, the 
connection of the optical and X-ray spectra remained unclear. 
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Figure 2.18: Possible non-dipolar structure of the magnetic field lines of PSR 60950+08. The 
structure was obtained using two crust anchored anomalies located at: ri = (0.95-R, 4°, 0°), 
ra = (0.95i?, 5°, 180°), with the dipole moments mi = (5.9 x lO^^d, 15°, 0°), ma = 
(5.9 X lO^^d, 140°, 0°) respectively (blue arrows). The influence of the anomalies is negligible 
at distances D > 5.0R, where B^/ B^ ^ m/d = 5.9 x 10^^ (top panel). For more details on the 



polar cap region see Figure 2.19 
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Figure 2.19: Zoom of the polar cap region of PSR B0950+08. See Figure 2.18 for a description 
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Figure 2.20: Dependence of a curvature of the open magnetic field lines on the distance from the 
stellar surface for PSR B0950+08. The distance is in units of stellar radius {z% = z/R) and the 
curvature of the magnetic field lines is p_6 = 1/% = p/ (lO^^cm^^). 
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2.4.6 PSR B1133+16 

Pulsar B1133+16 is one of the brightest pulsating radio sources in the Northern hemi- 
sphere (Maron et al. , 2000). The relatively long pulse period P = 1.19 s and its first de- 
rivative P_i5 ~ 3.5 result in the following inferred physical properties: B^ = 4.3 x 10^^ G, 
Tc = 5.04 Myr. The pulsar profile exhibits a classic double peak along with the usual 
S-shaped polarisation-angle traverse. The pulsar also shows the phenomenon of drifting 
subpulses but only for some finite time-spans, outside of which the behaviour of individual 



pulses is chaotic (Honnappa et al. 2012). 



PSR B1133+16 is located at a high galactic latitude, thus implying a low interstellar 



extinction (Schlegel et al. , 1998). Zharikov et al. (2008) suggested a possible optical 
counterpart with brightness B = 28™"^. 



X-ray observations performed by Kargaltsev et al. (2005) with the Chandra result in 



a small number of counts (33 counts from over 17 ks), thus the X-ray spectrum can be 
described by various models. The photon statistics are so low that they allowed only 
separate fits for the thermal (BB) and nonthermal (PL) components. 
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Figure 2.21: Possible non-dipolar structure of the magnetic field lines of PSR B1133+16. 
The structure was obtained using two crust anchored anomalies located at: 
ri = (0.95i?, 2°, 0°), r2 = (0.95i?, 5°, 180°), with the dipole moments 
mi = (8 X 10-3(i, 20°, 0°), m2 = (8 x lO^^d, 170°, 0°) respectively (blue arrows). The 
influence of the anomalies is negligible at distances D > 4.2i?, where B^/B^ ss m/d = 5 x 10~^ 



(top panel). For more details on the polar cap region see Figure 2.22 
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Figure 2.22: Zoom of the polar cap region of PSR B1133+16. See Figure 2.21 for a description 
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Figure 2.23: Dependence of a curvature of the open magnetic field lines on the distance from the 
stellar surface for PSR B1133+16. The distance is in units of stellar radius {zq = z/R) and the 
curvature of the magnetic field lines is p_6 = 1/% = p/ (lO^^cm^^). 
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2.4.7 PSR B1929+10 

With a pulse period of P = 0.23 s and a period derivative of P-15 ~ 1.2, the pulsar's 
characteristic age is determined to be r^ = 3.1 Myr. These spin parameters imply a dipolar 
component of the magnetic field at the neutron star magnetic poles B^ = 1.0 x 10^^. The 
distance to the pulsar D = 0.36 kpc was estimated using the parallax. 



Pavlov et al. (1996) identified a candidate optical counterpart of PSR B1929+10 with 



brightness U ~ 25.7'"'''^, which was later confirmed by proper motion measurements per- 



formed by Mignani et al. (2002). 



The X-ray pulse profile of PSR B1929+10 consists of a single, broad peak which is 



in contrast with the sharp radio one of Misanovic et al. (2008). The two-component 



spectral fit (BB+PL) suggests that both the thermal and nonthermal luminosities are of 
the same order. The derived surface temperature Tg = 4.5 MK and the surface magnetic 
field -Bs = 1.3 x 10^^ G do not coincide with the theoretical curve To — B^ of the critical 



temperature calculated by Medin and Lai (2008). We believe that this inconsistency can 



be removed by adding an additional blackbody component (the whole surface or the warm 
spot radiation). 
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Figure 2.24: Possible non-dipolar structure of the magnetic field lines of PSR B1929+10. 
The structure was obtained using two crust anchored anomalies located at: 
ri = {0.95R, 14°, 180°), r2 = {0.95R, 0°, 0°), rs = {0.95R, 14°, 0°), with the dipole moments 
mi = (1 X IQ-'^d, 20°, 180°), m2 = (2 x IQ-'^d, 180°, 0°), ms = (3 x IQ-'^d, 10°, 0°) respect- 
ively (blue arrows). The influence of the anomalies is negligible at distances D > 4.5R, where 



Bm/Bd ~ m/d = 3 X 10 (top panel). For more details on the polar cap region see Figure 2.25 
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Figure 2.25: Zoom of the polar cap region of PSR B1929+10. See Figure 2.24 for a description 
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Figure 2.26: Dependence of a curvature of the open magnetic field lines on the distance from the 
stellar surface for PSR B1929+10. The distance is in units of stellar radius {z% = z/R) and the 
curvature of the magnetic field lines is p_6 = 1/% = p/ (lO^^cm^^). 
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Chapter 3 

Partially Screened Gap 



The charge-depleted inner acceleration region above the polar cap can be formed if a local 



charge density differs from the co- rotational charge density (Goldreich and Julian, 1969). 



We assume that the crust of the neutron stars mainly consists of iron (leF^) formed at the 



neutron star's birth (e.g. Lai, 2001 ). Depending on the mutual orientation of f2 and ^, the 
stellar surface at the polar caps is either positively {ft ■ /j, < 0) or negatively {ft ■ /j, > 0) 
charged. Therefore, the charge depletion above the polar cap depends on the binding 
energy of either the positive geFe ions or electrons. In this thesis we consider the case of 
positively charged polar caps {ft- fi < 0). We assume that due to the high cohesive energy 
of iron ions, the positive charges cannot be supplied at a rate that would compensate for 



the inertial outflow through the light cylinder (see Medin and Lai (2006, 2007); Gil et al. 



(2007b)). This is actually possible if the surface temperature Tg is below the critical value 
Tcrit. Since the number density of the iron ions in the neutron star crust is many orders 
of magnitude larger than the co-rotational charge density (the so-called Goldreich- Julian 
density) pcj = ^ ■ B/ (27rc), then a thermionic emission from the polar cap surface is not 
simply described by the usual condition ej ~ kTg , where e^ is the cohesive energy and/or 
work function, Tg is the actual surface temperature, and k is the Boltzman constant. The 



outflow of iron ions can be described in the form (Gil et al. 2003 and references therein) 



Pi 

Pgj 



a 






(3.1) 



where pi < pcj is the charge density of the outflowing ions. As soon as the surface 
temperature Tg reaches the critical value 



T. 



crit 






(3.2) 



the ion outflow reaches the maximum value p\ = pcj- The numerical coefficient 
Ci = 30 ± 3 is determined from the tail of the exponential function with an accuracy 
of about 10%. Thus, for a given value of the cohesive energy, the critical temperature Tent 
is also estimated within an accuracy of about 10%. The cohesive energy is mainly defined 



by the strength of the magnetic field and was calculated by Medin and Lai (2006, 2007). 
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3.1 The Model 



As it follows from the X-ray observations (see Section 1.4), the temperature of the hot 
spot (which is associated with the actual polar cap) is more than 10^ K. As we mentioned 
above, in order to sustain such a high temperature bombardment by the backstreaming 
particles is required. But particle acceleration (and therefore the surface heating) is 



possible only if Tg < Tent- Gil et al. (2003) introduced the model of the Partially Screened 



Gap to describe the polar gap sparking discharge specifically under such circumstances. 

The PSG model assumes the existence of heavy iron ions (26?^) with a density near but 
still below the co-rotational charge density (pcj), thus the actual charge density causes 
partial screening of the potential drop just above the polar cap. The degree of screening 
can be described by screening factor 



77 = 1 - Pi/pGJ- 



(3.3) 



where pi is the charge density of the heavy ions in the gap. The thermal ejection of ions 
from the surface causes partial screening of the acceleration potential drop 



/W = rjAVm 



(3.4) 



where A\4iax is the potential drop in a vacuum gap. We can express the dependence of 
the critical temperature on the pulsar parameters by fitting to the numerical calculations 



of Medin and Lai (2007) 



T, 



crit 



1.6 X 10^ 



0.5 
PP-15 ) b 



1.1 



17.7 



(3.5) 



or Tcrit = 1.1 X 10^ (Bll + 0.3), where Bu = BJ {W^ G) , B^ = bB^ is a surface magnetic 
field (applicable only if hot spot components are observed, i.e. 6 > 1). 

The actual potential drop AV should be thermostatically regulated and a quasi- 
equilibrium state should be established in which heating due to the electron/positron 



bombardment is balanced by cooling due to thermal radiation (see Gil et al. 2003 for 
more details). The necessary condition for this quasi-equilibrium state is 



aT^ = rjeAVcncj, 



(3.6) 



where a is the Stefan-Boltzmann constant, e - the electron charge, and 
hq] = pGj/e = 1.4 X lO^^feP^fgP"'^'^ is the co-rotational number density. The Goldreich- 
Julian co-rotational number density can be expressed in terms of B14 as 



ncj = 6.93 X WBuP 



)-i 



(3.7) 
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Here we assume that the density of backstreaming relativistic electrons is tjugj- 



By using Equations 3.6 3.5 and |3.7| we can express the acceleration potential drop 



that satisfies the heating condition (Equation 3.6) as follows 



AV = 7.3 X 10 



,{Bli + 0.3)' P 



r]B 



14 



(3.8) 



The above equation may suggest that the acceleration potential drop is inversely pro- 



portional to the screening factor. In fact, it is just the opposite (see Equations 3.4 and 



3.26). 



Knowing that AV = 7max'"^c^/e, where m is the mass of a particle (electron or 
positron), we can calculate the maximum Lorentz factor of the primary particles in PSG 
as 



7n 



450 



(Bii + o.syp 



r]B 



(3.9) 



14 



3.1.1 Acceleration potential drop 

As the actual polar cap is much smaller than the conventional polar cap (see section 



1.4.4), we cannot use the approximation proposed by Ruderman and Sutherland (1975) 



that the gap height is of the same order as the gap width [h ^ h±). On the contrary, the 
small polar cap size and subpulse phenomenon suggest that in the PSG model the spark 
half- width is considerably smaller than the gap height {h± < h). For such a regime we 
need to recalculate a formula for the acceleration potential drop AV. 

Let us consider a reference frame co-rotating with a star and with the z-axis aligned 



with the star's angular velocity f2 (see Figure 3.1). 




C X 



Figure 3.1: Co-rotating frame of reference with the z-axis aligned with the angular velocity CI. 
The magnetic dipole moment /i is constant in this frame of reference, thus dB/dt = 0. 
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Let us underline that we will neglect the effects of non- inert iality of the co-rotating 
system. Thus, we assume that in any given moment we have a system moving with a 
constant velocity. 

In this co-rotating frame of reference we can write the spherical components of an 
angular velocity as follows 

n = {ncos9, -nsm9,0). (3.10) 

Gauss's law in the co-rotating frame (after Lorentz transformations) takes the form 

V ■ E = 47rp (r) - 47r . (3.11) 

\ 27rc / 

While Faraday's law of induction can be written as 

VxE = 0. (3.12) 

Note that if we consider a drift of plasma in the Inner Acceleration Region (lAR), we 



should expect temporal variations of the magnetic field (V x E = —d'B/{cdt)) (Schiff 



1939), but as was shown by van Leeuwen and Timokhin (2012), even if we consider 



fluctuations of the electric current of the order of the Goldreich- Julian current pcjc, the 
resulting variation of the magnetic field is so small that V x E = with a high accuracy, 
and circulation of the non-co-rotational electric field along a closed path is zero. 



Equation 3.11 in the spherical system of coordinates has the following form 



-Er + ^ + —.Ee + -^ + —n^ = 47rp r,^,0 -47r —— . 3.13 
r or rsmt7 r OU rsmt/ 0(p \ lire J 

The PSG model assumes the existence of ions in the lAR region that affects the charge 
density. Using the screening factor, rj, we can write that 

ri B 

P {r, 9, 0) = (1 - r]) pgj (r, 9, 0) = (1 - r/) ^^• 

In general, t] depends on the curvature and strength of the magnetic field, thus it 
varies across the polar cap, but we can still assume that rj is approximately constant at 
least for a given spark. Then 



2^ dEr cos9 ^ IdEe 1 dE^ ^ / Br^ cos 9 - Bg^ sin 9\ , ^, 

-i?, + — ^ + —.Eg + -^ + —a^ = -4vrr7 — . 3.14 

r or rsm.9 r o9 rsin.9 ocp \ Zirc J 

Let us change the variables as follows: r = R+ z and 9 = a + ^. Here R is the stellar 
radius and a is the inclination angle between the rotation and the magnetic axis. 
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2 ^ dEr 



COS (a + {}) 



-Ee + 



dEe 



+ 



dE. 



R + z dz {R + z)sm{a + ^) R + z d^ {R + z) sin {a + ^) 

/ (S^f) cos 61 - fief] sin 6*; . 

=-'"H ^^^^ '• ^ ^ 

Assuming that R^ z, which is correct as the gap height is less than the stellar radius 
{h '^ R), a ^ -&, and Bj. ^ Bg, which is correct for the polar cap region, we can write 



Equation 3.15 in the first approximation (R — ?■ oo) as follows 



dEr 



+ 



1 dEf, 



dz R d^ 



-Attt] 



Br^ cos 6 

27rc 



(3.16) 



Note that for spark widths considerably smaller than the stellar radius 



1 dEg ^, cot(a+tf) 



R 



Eg. 



h±<t:R (A^ ^ h±/R) we can write that ;^^ > 

Let us now consider Faraday's law (Equation 3.12). The curl of an electric field in 
spherical coordinates can be written as 



(V X E), 
(V X E), 
(V X E), 



1 



-{E.sme)- — 



rsinO \d9 

1 / 1 dEr _ d 

r V sin^ (90 dr 



(rE, 



1 fd_ 
dr 



[rEg 



dE, 
~d9 



0. 



(3.17) 



Using the same change of variables we performed above {r = R + z and 9 = a + d), 



the third equation of System |3.17| can be written as 

/?^^^-^^^ r3i8i 

From this equation in the zeroth approximation we can estimate the variations of the 
electric field components as 



R/\EgM ^ AEr-Az. 



(3.19) 



Since hi'^Rwe can write that 



{h±Eg) = {hEr) = AV. 



From Equation |3.16| we can also briefly estimate that 

.Bi.r2cos0 



AEr AEg 

-ir^ — 



/l, 



—Airr] 



2nc 



(3.20) 



(3.21) 



Using Equations 3.16 and 3.20 we can write that 
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(hEr) {h±Ee) _ AV AV 
h^ ^ hi ~~hF^ ~hl' 

Finally, we can estimate the potential drop in a spark region 



AV AV _ 27]BrVL cos (a + ^) 



(3.22) 



h^ 



h' 



(3.23) 



If we use the same assumption as Ruderman and Sutherland (1975), i.e.: (1) the spark 
half-width is of the same order as the gap height h± = h, (2) there is no ion extraction 
from the stellar surface (r/ = 1), and (3) the pulsar magnetic and rotation axes are aligned 



[a 



0°), we get: 



AV 



Bri'i 



RS 



h\ 



(3.24) 



Note that the potential drop defined by Equation |3.23| differs from that used in the 
Standard Model by the screening factor (as the presence of ions screens the gap) and by 
the factor of cos {a + •&) which also takes into account non-aligned pulsars. In our case the 
polar cap size is much smaller than the conventional polar cap size. It seems reasonable 
to also consider sparks with widths much smaller than the gap height {h_j_ <^ /i), in that 
case the potential drop can be calculated as 



Al^ 



2r]B,n cos {a + ^) 



h' 



(3.25) 



Even for a relatively small inclination angle between the rotation and magnetic axis, 
we can still write ^ <^ a, thus 



AV 



AnriBr cos a , n 
cP ^- 



(3.26) 



3.1.2 Acceleration path 

Since the exact dependence of the electric field on z is unknown we use the same linear 



approximation that Ruderman and Sutherland (1975) used. In the frame of the PSG 



model a^s h± < h or even h± ^ h, we can use Equations 3.20 and 3.26 to describe the 
component of the electric field along the magnetic field line: 

STTrjBs cos a h'^_ 



E 



cP 



K^ 



(h-z), 



(3.27) 



which vanishes at the top z = h. The Lorentz factor of particles after passing distance 



/ace can be calculated as follows 



7a 



TTIC^ 



^^ ^ , SnriBse cos a h?^ , 

Edz ^ ^^; -^ (-22 



21 



nic 



:3P h^ 



Zl)[h- 



Zl + Z2 



(3.28) 
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where m is the mass of a particle (electron or positron) and Z2 — zi = /ace- Then we can 
approximate zi + Z2 ^ h, thus 



7accy»C^-P h 

ATTrjBse cos a h'] 



''acc,ap . J-) J 2 ' yu.^uu) 



Assuming that a non-relativistic particle is accelerated from the stellar surface (zi = 0, 
7o = 1) we can calculate the distance /ace which it should pass to gain a Lorentz factor 

Tacc- 



/ill-Wl-^l, (3.30) 



where d = Sttt] Bseh\ cos (a) / (Pern). Although the approximate formula 3.29 is much 



more readable, in the calculations we use the exact value (see Equation 3.30 ) as for Lorentz 



factors that are considerably smaller than the maximum value, the discrepancy is about 

a factor of two, /acc.ap ~ 2/acc- 

3.1.3 Electron/positron mean free path 

The mean free path of a particle (electron and/or positron) /p can be defined as the mean 
length that a particle passes until a 7-photon is emitted. In the case of the CR particle, 
mean free path can be estimated as a distance that a particle with a Lorentz factor 7 



travels during the time which is necessary to emit a curvature photon (see Zhang et al 



1997) 



PcR \ ^ Qh'^c 



IcK-ci-^ =j-^, (3.31) 

V^7,CR/ 4 76^ 

where Pqr = 27^e^c/33f?^ is the power of CR, E^qr = Sh'y^c/2^ is the photon character- 
istic energy, and 3fJ is the curvature radius of the magnetic field lines. 

For the ICS process calculation of the particle mean free path lies is not as simple as 
that of the CR process. Although we can define lies in the same way that we defined Iqr, 
it is difficult to estimate the characteristic frequency of emitted photons. We have to take 
into account photons of various frequencies with various incident angles. An estimation 



of the mean free path of an electron (or positron) to produce a photon is in Xia et al. 

a' (e, fi) (1 - (3^i) nph (e) dedfi 



(1985) 



ICS 



Ml 



1 -1 



(3.32) 
'mo ■'o 

Here e is the incident photon energy in units of mc^, /i = cosip is the cosine of the photon 
incident angle, /3 = v/c is the velocity in terms of speed of light, a' is the cross section of 
ICS in the particle rest frame. 

An e2 

nph (e, T) de = — — — -de (3.33) 

A;? exp (e/0) — 1 
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represents the photon number density distribution of semi-isotropic blackbody radiation, 
O = kT/mc^, k is the Boltzmann constant, and Ac = h/mc = 2.424 x 10~^° cm is the 
electron Compton wavelength. A detailed description of how to calculate a' can be found 
in Section 14.4.11 



We should expect two modes of ICS: resonant and thermal-peak (see Section 4.4.3 
for more details). The Resonant ICS (RICS) takes place if the photon frequency in the 



particle rest frame is equal to the electron cyclotron frequency. As shown in Section 4.4.4 



the particle mean free path strongly depends on the distance from the polar cap. Both 
the photon density and incident angles (/io and /xi) change with increasing altitude. In 
our calculations we take into account both of those effects, thus we replace nph (e, T), /iq 
and yUi with n^p (e, T, L), fi^i^{L) and /imax (L), respectively (for more details see Section 
4.4). Here, L is the location of the particle, ngp (e, T, L) is the photon density at 
location L, and /Xmin (L) and /imax(-^) correspond to the highest and lowest angle between 
the photons and particle at a given location L. Thus, just above the polar cap for RICS 
the mean free path of outflowing positrons is: 



/ 



RICS 



max 

"res 



mm 

res 



1^1 — /3/i) a' (e, /i) Usp (e, T, L) dedjji 



(3.34) 



'/»min(i) 

where the limits of integration over energy, e""" and g™'"' are chosen to cover the resonant 



energy (for more details see Section [4.4.3) 



The thermal-peak ICS (TICS) includes all scattering processes of photons with fre- 
quencies around the maximum of the thermal spectrum. As an example we adopt 
e™'" ^ O.OSe^^, and e"^"" ~ 2ey^ where e^^ = 2KlkT j {m,(?^ is the energy, in units of 
m(? ^ at which blackbody radiation with temperature T has the largest photon number 
density. The electron/positron mean free path for the TICS process is 



/' 



TICS 



Qaxl-t^) 
/imiii(i) 



max 

'th 



mm 
th 



/3/i) o' (e, /i) nph (e, T, V) dedfi 



(3.35) 



3.1.4 Photon mean free path 

The photons with energy E^ > 2mc^ propagating obliquely to the magnetic field lines can 
be absorbed by the field, and as a result, an electron-positron pair is created. To describe 
the strength of the magnetic field we use /3g = B/Bg, where Bg = m^c' /eh = 4.413 x 10^^ G 
is the critical magnetic field strength. 



For strong magnetic fields (/3g > 0.2, see Section 4.2.3) the photon mean free path can 



be calculated as (see Section 4.2.4 for more details) 

2mc^ 



ph 



^- 



E^ 



(3.36) 
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while for weaker magnetic fields (/3q < 0.2) we can use an asymptotic approximation 



derived by Erber (1966) 



/ h = ^-^ ^ ^^ (±] 

^^ (e^ / he) mcB sin "^ V^X/ 



(3.37) 



X 



E^ 5sin\E' 
2mc^ Ba 



ix « 1) 



(3.38) 



where \E' is the angle of intersection between the photon and the local magnetic field. 

3.2 Gap height 

By knowing the acceleration potential drop in PSG AV^ we can evaluate the gap height 
h and the screening factor rj, which actually depends on the details of the avalanche pair 
production in the gap. First, we need to determine which process, Curvature Radiation 
(CR) or Inverse Compton Scattering (ICS), is responsible for the 7-photon generation in 
the gap region. In order to identify the proper process we need the following parameters: 
lace - the distance which a particle should pass to gain the Lorentz factor 7acc, ^p - the 
mean length a particle (electron and/or positron) travels before a 7-photon is emitted, 
and /ph - the mean free path of the 7-photon before being absorbed by the magnetic field. 



As mentioned above, PSG can exist if Equation |3.6| is satisfied. On the other hand, 
in order to heat the polar cap surface to high enough temperatures the high enough flux 
of back-streaming particles is required. By using Equations 3.8 and |3.26 we can find the 
relationship between the screening factor, the spark half-width and pulsar parameters 



Tjhj 



4.17 



{Bii + o.syp 



-Bi4V I 



cos a 



(3.39) 



Thus, for specific pulsar parameters we can define a product of the two main para- 
meters of PSG, namely the screening factor rj and the spark half-width h±. 



3.2.1 Particle mean free paths, CR vs. ICS gap 



The Figure |3.2 shows the dependence of particle mean free paths on the Lorentz factor 
7 for some pulsar parameters (the dependence on pulsar parameters will be discussed 



in Section 3.3). Let us note that these free paths do not depend on the gap height h 



see Equations 3.31, 3.32 and 3.36). The results presented in the Figure do not allow to 



define the gap height unambiguously. However, we can find which process is responsible 
for generation of the 7-photon in PSG. For narrow sparks the acceleration potential drop 
decreases, and as a result the Lorentz factor of the primary particles is about 7 ~ 10^ — 10^. 
In this regime lies ^ ^cR) so the gap will be dominated by ICS. Thus, ICS will dominate 
the gap if deceleration due to Inverse Compton Scattering prevents further acceleration 
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by the electric field. Let us remember that ICS is not efficient for particles with the 
Lorentz factor 7 ^ 10^. If the sparks are wider ot rj ^ 1, the acceleration potential drop 
increases and the Lorentz factors of primary particles reach values about 7 ~ 10^ — 10^. 
In this regime the 7-photon emission is dominated by CR. Let us note that the condition 
^ics ^ ^CR is satisfied for particles with 7 ~ 10^ — 10^, as one can see from Figure 3.2 
(panel a), but this does not mean that the ICS event happens. Since /ace ^ ^icS) the 
particles will be accelerated to higher energies (7 ~ 10^ — 10^) before they upscatter 
the X-ray photons. Thus, the particles start emission of 7-photons (via CR) as soon as 
condition /ace ~ /cr is met. 
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Figure 3.2: Dependence of the mean free path of the primary particle on Lorentz factor 7 for both 
the CR and ICS processes. Panel (a) corresponds to calculations for a relatively higher potential drop 
(e.g. a wider spark with h± = 3m. and t] = 1), while panel (b) corresponds to calculations for a 
relatively lower potential drop (e.g. a narrow spark with h± = Im. and t] = 0.1). The acceleration 
paths on both panels were calculated for the same pulsar parameters {B14 = 3.5, Tq = 4.4, JRg = 1, 
P = 1, a = 10°). Note that for the RICS process the particle mean free paths were calculated for 
optimal conditions (just above the polar cap). 



3.2.2 Possible scenarios of the gap breakdown: PSG-on and 
PSG-off modes 



As is seen from Figures |3.3| and |3.4[ in the CR-dominated gap the primary particle should 
travel a distance comparable with gap height /aec ~ h/2 in order to gain an energy 
corresponding to the characteristic Lorentz factor 7^, . On the other hand, the primary 

ICS 

particles in the ICS-dominated gap reach a characteristic value 7^, at altitudes that 
are considerably smaller than gap height /ace ^ h (see Figures 3.3 and 3.5). Thus, 



7, 



7n 



is about three orders of magnitude higher than 7^ 



10^ < 7max, here 



7max is the value of the Lorentz factor after the particle travels a distance h. Furthermore, 
the characteristic energy of CR photons is considerably smaller than the energy of emitting 
(primary) particles, e.g. for 7 = 10^, ^q = 1, 7sec ~ 10^. On the other hand, RICS photons 
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upscattered in an ultrastrong {B > i?crit) magnetic field gain a significant part of the 
energy of the scattering (primary) particle. Therefore, the electron/positron pair created 
by the RICS photon has energies comparable with the energy of the scattering (primary) 
particle. This will essentially influence the multiplicity Mpr in the ICS gap, as all the 
newly created particles will participate in further cascade pair-production. Additionally, 
RICS in ultrastrong magnetic fields produces approximately the same amount of photons 
with II and _L polarisation (see Section 4.4.2), while most of the photons produced by CR 
are ||-polarised (see Section 4.1). Splitting of the ±-polarised photons will increase the 
photon mean free path, but it will also increase the multiplicity in the ICS gaps. 

Figure [373 presents a sketch of a cascade formation for CR- and ICS-dominated gaps. 
The CR photons are emitted in the upper half of the gap. Most of these photons produce 
pairs at about the same height, in the region where the acceleration potential is almost 
equal to zero, hereinafter we will call this region the Zero-Potential Front (ZPF). The 
newly created particles have much lower Lorentz factors as compared with the primary 
particle, thus they are not able to emit CR photons. 



z ~ 100m fs hcR p+ 




z ~ 10m <S hi, 




Figure 3.3: Sketch of differences in a cascade formation for the CR-dominated gap (left panel) 
and the ICS-dominated gap (right panel). In order to increase readability, only a few points (filled 
circles) are shown which correspond to altitudes where 7-photons are emitted. The unfilled circles 
correspond to places where 7-photons are also emitted, but those photons (and their evolution) are 
not included in the diagram. Note that for the ICS-dominated gap we plot only the bottom (active) 
part of the gap {z <C /iics). furthermore, points of radiation are tracked only for the first population 
of newly created particles. The avalanche nature of the ICS-dominated gap will result in a much 
higher multiplicity and continuous backflow of relativistic particles. 
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Figure 3^ presents the primary particle evolution and photon mean free paths of 
7-rays produced in the CR-dominated gap. As can be seen, the first 7-photon produces 
a pair approximately at the same time (and same place) as the primary particle reaches 
ZPF. Thus, the multiplicity in a gap region (the number of particles created by a single 
primary particle) in the CR scenario is strictly related to the number of photons produced 
by the primary particle Mcr ~ 2 x A^p^ • 
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Figure 3.4: Cascade formation for a CR-dominated gap. Blue lines represent the mean free path 
of 7-photons. The filled circles correspond to places of 7-photon emission. Panel (a) includes the 
free paths of 7-photons which produce pairs below ZPF (red circles) while panel (b) includes the 
free paths of 7-photons which produce pairs above the acceleration gap (blue circles). The results 
were obtained using the following parameters: N^^ = 50, B^ = 2.3 x lO^'^G, Bd = 2.9 x IQi^G, 
r= 3MK, P = 1.3s, Ke = 0.5, and a = 60.7°. 



The energy of 7-photons produced by ICS depends on the Lorentz factor of the primary 
particles and on the strength of the magnetic field. In ultrastrong magnetic fields the 
energy of newly created particles is comparable with the energy of the scattering particle 
7new ~ 7c/2. Figure 3.5 shows schematically the locations at which 7-photons are emitted 



by ICS. The first 7-photon is produced already at altitudes of about a few metres and then 
converted to an electron-positron pair well below ZPF. Note that already at relatively 



low altitudes {z > 100 m) the photon density decreases rapidly (see Section 4.4.4.1) 



furthermore, the small size of the polar cap entails a rapid change of the particle-photon 
incident angles (see Section 4.4.4.2[ ). Those two effects make the ICS process significant 
only in the lower parts of the gap {z < 20 m). On the other hand, the multiplicity in the 
ICS-dominated gap is enhanced by all newly created particles which are created in the 
lower part of the gap. Furthermore, the ICS is more effective for backstreaming particles 



see Figure 3.2 ), thus most of the 7-photons in the gap region will be created by scatterings 



on electrons. For the ICS scenario it is not possible to evaluate a simple expression for 
the multiplicity produced by a single primary particle in a gap region. Furthermore, it is 
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not possible to determine the actual value of iVph required to break the gap (both for CR 
and ICS) without a full cascade simulation. 
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Figure 3.5: Cascade formation for an ICS-dominated gap. Blue lines represent the mean free path 
of 7-photons. The filled circles correspond to places of 7-photon emission. The results were obtained 



.ICS 



using the following parameters: N^^ = 15, B^ = 2.3 x lO^^G, Ba = 2.9 x lO^^G, T^ = 3MK, 
P = 1.3s, Ke = 0.5, and a = 60.7°. 

The differences between the CR and ICS gaps that we mention above have drastic con- 
sequences on the cascade formation process. Since the cooling time of the hot spot is very 
short (tcooI ^ 10^* s , see Gil et al. , 2003), to sustain the hot spot temperature just below 



the critical temperature a continuous backflow of relativistic particles is required. An 
energetic enough flux of backstreaming particles can be produced only in ICS-dominated 
gaps. The heating of the surface will sustain the outflow of iron ions from the crust, 
maintaining 77 < 1, hence we call this mode the PSG-on mode. As the temperature of 
the polar cap is in quasi equilibrium with the backstreaming particles (temperature is 
close to the critical value) the gap can break only due to production of a dense enough 
plasma Up ^ fjUQj in the gap region. The multiplicity in the PSG-on mode is much higher 
than the multiplicity of CR-dominated gaps. Moreover, in the gap dominated by CR the 



particles are created in a cloud-like fashion (see Figure 3.4). The successive clouds heat up 
the surface once per tq ~ 2h/c, which for a typical gap height h ^ 100 m is much longer 
than the time needed for the surface to cool down tq ~ 6 x 10^^ ^ TcooI- Therefore, in the 
CR-dominated gaps the backstreaming particles cannot sustain the temperature that is 
close to the critical value during ri ^ tq ^ TcooI, thus for most of the time the screening 
factor is ?7 ^ 1 and we call this mode the PSG-off mode. The low multiplicity of a cascade 
in the PSG-off mode can cause that the gap to breakdown only due to overheating of 
the surface, but not due to production of a dense enough plasma. The growth of particle 
density will continue to the point when the backstreaming particles heat up the surface 
to a temperature equal to or higher than the critical temperature, Xheat ^ tq- Let us note 
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that the primary particles in the PSG-off mode are very energetic 7 ~ 10^, and hence the 
density of particles required to close gap pc is much lower than the Goldreich- Julian dens- 
ity. To describe this difference we use the overheating parameter k = Pc/pcj- Knowing 
that in the PSG-off mode 77 ?^ 1, we use Equation 3.6 and the relation AV = 7acc^™c^/e 
to calculate the overheating parameter: 



aT^ 



^GJ 7max mc" 



(3.40) 



3.2.3 PSG-ofF mode 



Curvature emission by a primary particle is effective for Lorentz factors 7 > 10^ (when 
^CR < 4cc)- An equilibrium between acceleration and deceleration (by reaction force) 
would be established if the CR power were equal to the "electric power". In our case 
{^Q ^ 1, 7j, Ri 10^), the reaction force is not high enough to stop acceleration by the 
electric field. In the PSG-off mode the spark region is free from ions [r] ^ 1), thus the 



heating condition (Equations 3.6 and 3.39) is no longer satisfied. Taking into account 



the curvature of magnetic field lines just above the stellar surface, we can estimate the 



dependence of the minimum spark half- width on the gap height (see Figure 3.6): 



^-VW^^¥. 



(3.41) 



Figure |3]7] presents the minimum spark half-width for three different radii of curvature: 
^Q = 0.1, ^Q = 0.5, Sf^e = 1- Note that as long as the gap height does not exceed some 
specific value [h ^ 40m, h ^ 100 m, h ^ 140 m, respectively for the given curvature radii) 
the minimum spark half-width is well below 1 m. 




Figure 3.6: Diagram of the minimum spark half-width h^ for a given gap height h and a radius 
of curvature R. 
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Figure 3.7: Minimum spark half-width vs. gap height calculated for three different radii of 
curvature: Kg = 0.1 - red solid line, Kg = 0.5 - green dashed line, and Kg = 1 - blue dotted 
line. 

On the other hand, we can estimate the acceleration potential AV^ (and thus the spark 
half- width h^^ ) required to produce a specified number of photons A^pj^ within a gap. 



Figure 3.8 presents the dependence of both /i™'" and hj^ on the gap height. As results 



^CR 



and 



from the Figure, the gap height in PSG-off does not change drastically with A^ j^ , 
for these specific parameters of a pulsar it is /i ~ 240 m. For historical reasons, hereafter 
unless stated otherwise, we will use A^pj^ = 50 to calculate the gap parameters of the 
PSG-off mode. Note that in order to find the gap height, we assume /i™'" = h^^"^ , which 
results in a gap that allows both overheating of the entire spark surface by backstreaming 
particles and the creation of the required number of photons A^ j^ . 




H 6 



?i"™ (from curvature 
h^ (to meet the N^i^ = W condition) 
h^ (to meet tine JVph=50 condition) 
/i^ (to meet thie TVp,, =100 condition) 



Figure 3.8: Dependence of a spark half-width on the gap height for the PSG-off mode. The results 
were obtained using the following pulsar parameters: B14 = 2.3, Tq = 3.0, P = 1.3s, JRg = 1-0, 
a = 60.7°. 
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In our calculations we use the algorithm presented in Figure 3.9| to find the gap height 
in the PSG-off mode for given pulsar parameters: a pulsar period P, a pulsar inclination 
angle a, a surface magnetic field strength Bg, and a curvature radius of field lines 3fJ. 



Initial parameters: 
P, a, Bs, K 



Shielding factor: 77 = 1 



Initial gap height: h = /i;, 




Spark width: /ix = hj^ (h) 



Calculation of number of scatterings: A^pi, 



Higher gap: h — h + Ah 




Final gap parameters: h, h±, n 



Figure 3.9: Flowchart of the algorithm used to estimate the gap height in the PSG-ofF mode. The 
initial gap height from which we begin our calculations is an arbitrary set to /linit. = 10m, while the 
step Ah depends on the required accuracy. The number of 7-ray photons created in a spark by a 
single primary particle is set to A'pj^ = 50 (see text for more details). 



Figure 3. 10| presents the result of finding the gap height in the PSG-off mode for PSR 
B0943-I-10. The presented solution corresponds to the magnetic field structure presented 
in Section 2.4.4] The average radius of curvature in the gap region is relatively high, 
3f?6 = 0.7, hence the inclination of the gap region. The polar gap conditions, the strength 
of magnetic field Bi^ = 2.4 (i?bb = 17m) and the polar cap temperature Tg = 3.0 



were restrained to follow the observed values (see Table 1.4). The presented solution 
corresponds to the following PSG parameters: gap height h = 166 m, spark half-width 



1.9 m, T] = 1 (fixed), k = 7 x 10" 



7c 



1.4 X 10^. Note that the primary particles 



will gain ■j„ 



1.9 X 10^ as the CR efficiency is not high enough to stop the acceleration. 
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Figure 3.10: Gap structure in the PSG-off mode for PSR B0943+10. Filled columns represent the 
locations and sizes of the active regions of sparks. Here we assumed that the active region of a spark 
(the place where acceleration is high enough to produce a cascade) has a size comparable with the 
spark half-width. The iron ions extracted from the surface (due to a high surface temperature) are 
represented by circle-plus symbols. 



3.2.4 PSG-on mode 

In the PSG-on mode, radiation of the surface just below the spark is in quasi-equihbrium 
with the flux of backstreaming particles. When the surface temperature rises, the density 
of iron ions increases, thus resulting in a decrease in the potential drop, which in turn, 
reduces the flux of backstreaming particles. On the other hand, when the surface tem- 
perature decreases it entails the drop of iron ion density and, consequently, an increase 
in the flux of backstreaming particles. Thus the polar cap temperature is maintained 
slightly below the critical value. This quasi-equilibrium state prevents the gap breakdown 
due to surface overheating. However, a high multiplicity in the PSG-on mode leads to a 
production of dense plasma. When the density of the plasma n^ ^ ^^gj, the acceleration 
potential drop will be completely screened due to charge separation. 

Alongside the pulsar parameters the gap height in the PSG-on mode also depends 
on the spark half-width h± and on the number of scatterings by the first population of 

ICS 

newly created particles A^pj^ . For a sample of pulsars we can use drift information to put 



constraints on the spark half-width (see Section 3.4). Figure 3.11 presents the procedure of 
finding the gap height in the PSG-on mode for the following pulsar parameters: a pulsar 
period P, a pulsar inclination angle a, a surface magnetic field strength B^, a surface 
temperature Tg, a curvature radius of magnetic field lines 3fJ, and a spark half- width h±. 



First we use Equation |3.6| to estimate the screening factor r] which defines the electric 
field, and thus the particle acceleration. Then we estimate the number of scatterings for a 
single outflowing particle N^l^ for the initial gap height. The initial gap height from which 
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we begin our calculations is an arbitrary set to /imit. = 10 m. We track the propagation of 
7-photons produced by ICS on a primary particle to find the location L-^^w where pairs 
are created. Then we calculate their propagation through the acceleration region and we 
estimate the number of scatterings by every newly created particle of the first population 
A^p^". If the total number of scatterings by the first population (including the primary 

ICS 

particle) is A^ph < A^p^ , we resume our calculations assuming a higher gap until the 
A^ph > A^ph is met. 



Initial parameters: 



Shielding factor: 77 {Eqn. 3.39) 



Initial gap height: h = hinit. 



Primary particle progression: N^^^-, 7sec, ^st 



i = Q 



First population progression: N-^^ (7scc,-^scc) 




Final gap parameters: 

h, h^, r], 7c, 7maxi 7i: 



Figure 3.11: Flowchart of algorithm used to estimate the gap height in PSG-on mode for a given 
spark half-width (see text for more details). 
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As a result we obtain the gap parameters: the gap height /i, the screening factor t], 
the characteristic Lorentz factor of a particle at the moment of ICS photon emission 7^., 
the maximum value of the Lorentz factor 7max, and the characteristic Lorentz factor of 
iron ions 7i. In our calculations, if not stated otherwise, we use iV, 



ph 



25 to calculate 



the gap parameters of the PSG-on mode. Note that in this approximation we take into 
account only the first population of newly created particles. In fact, the avalanche nature 
of the ICS-dominated gap will result in a much higher multiplicity than in the PSG-off 
mode Mies ^ Mqr. For details of particle/photon propagation, see Chapter |4j Subpulse 
drift observations are available only for a few X-ray pulsars with the hot spot component. 
Thus, to find the approximate gap parameters for pulsars without the predicted spark 
half-width we use h± = 2in. 



Figure 3.12 presents the result of finding the gap height in the PSG-on mode for PSR 
B0943-I-10. In this model the gap parameters, such as the magnetic field strength Bg 
and the surface temperature Tg, were restrained to follow the observed values (see Table 



1.4). The result was obtained for the non-dipolar structure of a surface magnetic field 
presented in Section 2.4.4 and for the predicted value of a spark half-width h_j_ ^ 2Ta (see 

I L ICS 

Table 3.3). The height required to produce N^^ = 25 photons by the first population of 
particles was estimated as /i ~ 92 m. Other gap parameters for this solution can be found 
in Table O 




Figure 3.12: Gap structure in the PSG-on mode for PSR B0943+10. Filled columns represent the 
locations and sizes of the active regions of sparks. Here we assumed that the active region of a spark 
(the place where acceleration is high enough to produce a cascade) has a size comparable with the 
spark half-width. The iron ions extracted from the surface (due to a high surface temperature) are 
represented by circle-plus symbols. Note that iron ions are present in both the non-active and active 
regions. The density of ions in the non-active regions is so high that it prevents cascade formation 
of pairs. 
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3.2.5 Results 



In Table 3.1 we present the results of finding the gap height for the sample of pulsars. 
For the PSG-on mode we show the estimated PSG parameters found using the predicted 
spark half-width and the spark half-width h±^ = 2 m. The only exception is Geminga 
(PSR J0633+1746), for which drift information is not available and we can only present 
calculations for h^ = 2 m. For PSR B0628-28 the predicted spark half-width is large 
{h± = 3.9 m), which entails a high acceleration potential. For such wide sparks it is not 

ICS 

possible to find the PSG-on solution with the required number of scatterings N^ . We 
believe that for this specific pulsar the predicted spark half-width is overestimated. Ac- 
tually, if a spark is narrower {h± = 2 m), it can operate in the PSG-on mode (see Table 



3.1). This result may suggest that for this specific pulsar the parameters of the sub- 
pulse phenomenon could be overestimated (e.g. due to aliasing). On the other hand. 
X-ray observations of Geminga suggest a relatively low temperature of the hot spot 
(Ts ~ 1.9 MK). The low density of the background photons requires the formation of 
narrow sparks {h_j_ = 1 m) to allow the gap to operate in the PSG-on mode. We believe 
that the relatively large hot spot {Rpc = 44.5 m) of Geminga causes the width of the 
sparks to grow so fast that it can operate only in the PSG-off mode. We believe that this 
can explain the very weak radio luminosity of the Geminga pulsar. 



Table 3.1: Estimated parameters of PSG for the sample of pulsars. The conditions in the polar 
cap region: surface temperature, magnetic field strength, polar cap radius, and curvature radius 
of the field lines are given the in headers next to the pulsar name. The individual columns are as 



follows: (1) PSG mode (see Section 3.2.2 ), (2) Gap height, (3) Spark half-width, (4) Screening factor, 
(5) Overheating parameter, (6) Characteristic Lorentz factor of scattering particles , (7) Maximum 
Lorentz factor of primary particles, (8) Lorentz factor of iron ions (if they are relativistic), (9) Particle 
mean free path, and (10) Photon mean free path. The results are presented for two different gap 



breakdown scenarios: the PSG-off and PSG-on modes (see Section 3.2 for more details) 



"• The modes correspond to calculations using the predicted spark half-width (see Table 3.3) 
^ The modes correspond to calculations with a spark half-width h± = 2m. 
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Table 3.1 - continued from previous page 



mode 



h 

(m) 



hi. 

(m) 



K 



7c 



O'max Ti 



(m) 



(m) 



PSR B0628-28 Tr = 2.t 



off 



on 



on 



Bu = 2.2 /2pc = 21.3m ^e = 0.6 



198.3 


3.2 




0.007 


- 


3.6 


- 


- 


78.6 


2.0 


0.15 


— 



1.3 X 10^ 



6.1 X 10' 



1.6 X 10^ 



8.9 X 10^ 23 



1.4 



1.9 



58.9 



1.3 



PSR J0633+1746 Tg = 1.9 B^ = 1-5 



Rpc = 44.5 m 



3fifi = 2.1 



off 



on 



252.1 



1.5 
2.0 



0.0002 



2.9 X 10 



6 



3.5 X 10 



6 



2.2 



66.6 



PSR B0834+06 Tg = 2.4 5^ = 1.9 i?pc = 22.7m % = 0.6 



off 


172.5 


3.2 




0.0027 


1.2 X 10*^ 


1.4 X 10^ 




1.3 


52.1 


on*^ 


82.0 


1.8 


0.12 


- 


5.3 X 10=^ 


7.0 X 10^ 


18 


2.3 


1.5 


on^ 


102.9 


2.0 


0.11 




4.9 X 10^ 


7.8 X 10^ 


20 


2.4 


1.8 



PSR B0943+10 T6 = 3.1 ^u = 2.4 i?pc = 17.6 m ^f^g = 0.7 



off 



on 



168.3 
71.6 



1.9 
2.0 



0.1 



0.0068 



1.6 X 10*^ 
8.8 X 10^ 



2.0 X 10^ - 
2.5 X 10^ 63 



1.4 
1.1 



46.5 
1.1 
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Table 3.1 - continued from previous page 



mode 



(m) 



hi. 

(m) 



K 



7c 



7n 



7i 



(m) 



(m) 



PSR B0950+08 76 = 2.6 Bu = 2.0 i?„c = 14.0 m ^q = 



"-pc 



off 



on 
on 



a 



172.8 


1.9 




0.0009 


16.9 


0.7 


0.09 


- 


61.7 


2.0 


0.03 


— 
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B 
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Rpc = 17.9 m 



3fifi = 0.6 



off 
on" 
on* 



167.4 

95.9 

54.3 



2.4 
2.9 
2.0 



0.08 
0.11 



0.0076 
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1.3 X 10^ 33 
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PSR B1929+10 
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B 



14 



2.4 



R 



pc 



20 m 



^f, = 0.6 



oflF 


112.7 


1.0 




0.0012 


1.8 X 10^ 


2.1 X 10^ 




1.1 


28.0 


on'* 


50.5 


1.6 


0.02 


- 


9.8 X 10^ 


1.2 X 10^ 


31 


1.6 


1.0 


on'' 


75.1 


2.0 


0.02 
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8.4 X 10^ 


1.5 X 10^ 


39 


1.5 


1.5 
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3.3 PSG model parameters 

We can distinguish two types of PSG parameters: observed and derived. As we have 
mentioned above, in some cases when X-ray observations are available we can directly 
estimate the surface magnetic field Bg. On the one hand, Bg can be calculated using the 
size of the hot spot Abb, and on the other hand we can find Bg by using the estimation of 
the critical temperature and the assumption that Tg = Tent- One of the most important 
requirements for the PSG model is that these two estimations should coincide with each 



other. As is clear from Figure 1.11, in most cases when the hot spot parameters are 
available this requirement is fulfilled. Thus, we can assume that the characteristic values 
of Bg vary in the range of (1 — 4) x 10^^ G, which corresponds to the critical surface 



temperature in the range of (1.3 — 5) x 10 K (see Table 1.4). By using these values we 
can estimate the derived parameters of PSG, such as the gap height h, the screening 
factor rj (or the overheating parameter k in the PSG-ofT mode) and the characteristic 
Lorentz factor of primary particles 7c. Let us note that these parameters also depend on 
the curvature radius of the magnetic field lines 3fJ. The curvature can be neither observed 
nor derived, but modelling of the surface magnetic field (see Chapter |2]) indicates that 
the curvature radius varies in the range of (0.1 — 10) x 10^ cm. Below we will discuss the 
influence of pulsar parameters, such as the magnetic field, the curvature of field lines and 
the period on derived PSG parameters. 



3.3.1 Influence of the magnetic fleld 



The conditions in PSG are mainly defined by the surface magnetic field. In Figure 3.13 



panel (a) we present the dependence of the gap height on the surface magnetic field 



calculated according to the approach described in Section 3.2 It is clear that in the PSG- 
off mode the gap height decreases as the surface magnetic field increases. In the PSG-on 
mode, on the other hand, the gap height shows a minimum at a specific value of the 
magnetic field strength (for a given pulsar's parameters it is B14 ^ 3). This behaviour is 
the result of an increasing potential acceleration drop with an increasing surface magnetic 
field. When the magnetic field strength exceeds the optimal value, which corresponds to 
acceleration when ICS is more effective, the increase in the acceleration potential results 
in less effective scattering. Panel (b) shows the dependence of the screening factor (or 
the overheating parameter k in the PSG-off mode) on the surface magnetic field. We 
can see that for stronger magnetic fields both rj and k, increase, which means that: (1) 
the density of heavy ions above the polar cap in the PSG-on mode decreases, (2) the 
density of particles required to overheat (and thus to close) the polar cap increases. Let 
us note that the surface temperature Tg stays very near to the critical temperature Tent, 
which is shown on the top axis of the Figures. In panel (c) the red, solid and dotted 
lines correspond to characteristic and maximum Lorentz factors (7c, 7max) in the PSG-on 
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mode, while the blue, dashed and dashed-dotted lines correspond to 7c and 7max in the 
PSG-off mode. We see that especially for the PSG-off mode 7c does not depend on the 
magnetic field strength. Note also that in the PSG-off mode (CR-dominated gap), the 
characteristic Lorentz factor (the Lorentz factor for which most of the gamma photons 
are produced) slightly differs from the maximum value, 7c ~ 7max- On the other hand, 
in the PSG-on mode 7c <^ 7max; which reflects the fact that most of the scatterings take 
place in the bottom part of the gap. 
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Figure 3.13: Dependence of the gap height (panel a), the screening factor or the overheating 
parameter (panel b), and the particle Lorentz factor (panel c) on the surface magnetic field. Solid 
red lines correspond to the PSG-on mode (ICS-dominated gaps) while dashed blue lines correspond to 
the PSG-off mode (CR-dominated gaps). Calculations were performed using the following parameters: 
P = 0.23, Ke = 0.6, B^ = 1.2 x lO^^ q and a = 36°. The actual polar cap radius was calculated 
separately for a given surface magnetic field as Rpc = Rdp\/Bd/Bs. In panel (c) the red solid and 
dotted lines correspond to characteristic and maximum Lorentz factors (7c, 7max) in the PSG-on 
mode while blue dashed and dashed-dotted lines correspond to 7c and 7max in the PSG-off mode. 
Corresponding critical temperature is shown on top axis of the figures. 



3.3.2 Influence of the curvature radius 

The curvature of the magnetic fleld lines signiflcantly affects the gap height in the PSG-off 



mode (see Figure 3.14, panel a). In the case of the CR-dominated gap, the curvature of 
the magnetic fleld lines affects not only the photons' mean free path (for higher curvature 
the magnetic fleld will absorb photons faster), but also the particle mean free path and, 
more importantly, the energy of photons generated in the gap region. The higher energy of 
photons further reduces the photon mean free path, thus resulting in lower heights of the 
PSG. In contrast, the gap height in the PSG-on mode is only slightly affected by changes 
in the curvature of the magnetic fleld lines. In this case the most important parameter 
which determines the cascade properties is the primary particle mean free path which 
does not depend on the curvature of the magnetic fleld lines. 

The overheating parameter in the PSG-off mode inversely depends on the radius of 



curvature of the magnetic fleld lines (see Figure 3.14, panel b). The higher the curvature 
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the higher the overheating parameter, which means that the sparks are narrower. This 
is consistent with the expectation that for a higher curvature of the magnetic field hues, 
the gap breakdown is easier to develop and takes place before the sparks manage to grow 
in width. On the other hand, the screening factor in the PSG-off mode does not depend 
on the curvature of the magnetic field lines. 

With an increasing radius of curvature the Lorentz factor of primary particles (both 7c 
and 7max) required to close the gap in the PSG-off mode also increases. This reflects the 
fact that in order to produce a sufficient number of photons in the gap region, the primary 
particles should be accelerated to higher energies (if the curvature is lower). Higher 
energies of the primary particles will increase the emitted 7-photon energy, thereby they 
will partly inhibit the growth of the photon mean free path due to the lower curvature. As 
mentioned above, the gap height in the PSG-on mode very weakly depends on the photon 
mean free path, thus both 7c and 7max are not affected by the increase in the radius of 
curvature. 





Figure 3.14: Dependence of the gap height (panel a), the screening factor or the overheating 
parameter (panel b), and the particle Lorentz factor (panel c) on the curvature radius of magnetic 
field lines. Calculations were performed using the following parameters: P = 0.23, B^ = 1.2xlO^^G, 

36° 



Bs = 2.4x 10^^ G, a 



For a more detailed description see Figure 3.13 



3.3.3 Influence of the pulsar period 



As we can see from Figure 3.15 panel (a) and panel (c), in the PSG-on mode the gap 
height and the Lorentz factor of primary particles do not depend on the pulsar period. 



The increase in the screening factor (see Figure 3.15 d) compensates the increase in the 



acceleration potential drop (see Equation 3.8). Thus the particles in the gap region are 
accelerated in the same way independently of the pulsar period. On the other hand, 
the gap height in the PSG-off mode increases with the increasing pulsar period. This 
reflects the fact that in the PSG-off mode the acceleration potential, and hence 7c and 
7max, decreases with longer periods (see Equation 3.26). Longer pulsar periods entail 



an increase in the screening factor (in the PSG-on mode, see Equation 3.39) and in the 



overheating parameter (in the PSG-off mode). Note that for periods longer than some 
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specific value (for a given pulsar's parameters it is Pmax ~ 9 s), the screening factor in 
the PSG-on mode would exceed unity. This means that the PSG-on mode cannot be 
responsible for the gap breakdown for pulsars with such long periods. 
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Figure 3.15: Dependence of the gap height (panel a), the screening factor or a overheating 
parameter (panel b), and the particle Lorentz factor (panel c) on the pulsar period. Calculations 



were performed using the following parameters: P = 0.23, Ba = 1.2 x lO^^'G, B^ = 2.4 x 10^* G, 
a = 36°, Kg = 0.6. The actual polar cap radius was calculated separately for a given pulsar period 
as i?pc = Rdp^/Bd/Bg, where R^p = ^2-kR? j (cP). For a more detailed description see Figure 

EH 



3.4 Drift model 

The existence of lAR in general causes a rotation of the plasma relative to the NS, as 
the charge density differs from the Goldreich- Julian co-rotational density. The power 
spectrum of radio emission must have a feature due to this plasma rotation. This feature 
is indeed observed and is called the drifting subpulse phenomenon. 



3.4.1 Aligned pulsars 



An explanation for drifting subpulses was offered by Ruderman and Sutherland ( 1975 ) as 



being due to a rotating carousel of sub-beams within a hollow emission cone. According 
to this model a pair cascades may not occur simultaneously across the whole polar cap 
but is localised in the form of discharges of small regions in the polar gap. Such sparks 
may produce plasma columns that stream into the magnetosphere to produce the observed 
radio emission. The location of the discharges on the polar cap determines the geometrical 
pattern of instantaneous subpulses within a pulsar's integrated pulse profile. 

In the PSG model the stable pattern of subpulses is due to heating of the inactive 
part of the spark (the place where no cascade forms due to a low acceleration potential) 
by all the neighbouring discharges. The lifetime of a single spark is very short. On the 
other hand, an inactive region is continuously heated by all the neighbouring sparks. Even 
when one of them dies, the temperature is still high enough (high ion density) to prevent 
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spark formation in this region. As the discharges do not exchange information (they are 
not synchronised) and their hfetime is very small, the geometrical pattern of sparks on 
the polar cap should be stable. 

For pulsars with an aligned magnetic and rotation axis the sparks circulate around the 
rotation axis. Note that this circulation is not related with the magnetic axis but with 
the direction of the co-rotational velocity. Namely, the drift velocity is opposed to the 
co-rotational velocity. 

We can calculate the drift velocity of an aligned pulsar using the following approxim- 



ation (see Figure 3.16) 



where Rpc is the actual polar cap size, P2 is the characteristic spacing between subpulses 
in the pulse longitude, P3 is the period at which a pattern of subpulses crosses the pulse 
window (in units of the pulsar period), (3 is the impact angle, and p is the opening angle. 




a«0' 



Figure 3.16: Top view of a polar cap region of an aligned pulsar. Small circles represent sparks, 
while the red line corresponds to the line of sight. If we neglect the transition from a non-dipolar 
structure of the magnetic field on the stellar surface to a dipolar structure in the region where radio 
emission is produced, we can assume that the observed subpulse separation P2 also describes spark 
separation Qs (angular separation between the adjacent sparks on the polar cap). 

In such an approximation the assumption that only half of the spark is active can be 
written as 



2/1, 



PI ^_ 

360°'^27rPp4 



(3.43) 



Finally, we can define the observed drift velocity of aligned pulsars as 



Vdr 



2hi_ 

pK 



(3.44) 
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3.4.2 Non-aligned pulsars 



Most observed pulsars are non-aligned rotators. It is very common to apply the carousel 



model to interpret observations of the drifting subpulses of non-aligned pulsars (Ruder- 



man and Sutherland 1975[ ). Despite the fact that the carousel model can explain some 
properties of subpulses (for example the change in intensity), we believe that this model 
is not suitable for describing the spark's behaviour on the polar cap. There is no phys- 
ical reason for a spark to circulate around the magnetic axis. The circulation in aligned 
pulsars is caused by a lack of coronation with respect to the rotation axis. For non-aligned 
pulsars, the co-rotation velocity in the polar cap region has more or less the same direc- 
tion: what is more, if we assume circulation around the magnetic axis we will get plasma 
with a velocity that is higher than the co-rotational velocity, which is difficult to explain 
in a region where the charge density is lower than the co-rotational density. 

As in our model, the drift is caused by a lack of charge in lAR, thus the plasma should 
drift in approximately the same direction, i.e. in the direction opposite to the co-rotation 
velocity. We believe that the change in subpulse intensity is caused by the observation 
of a different part of a spark and/or different conditions across the polar cap at which 
the spark is formed (e.g. magnetic field strength, curvature of the magnetic field lines, 
background photon fiux). 

For pulsars with a relatively high inclination angle a we can calculate the drift velocity 
using the following approximation 



VdT 



'^Rpc^.P^ _2i?pe P^ 



PP^ W PP:, Wr 



(3.45) 



where W is the profile width and VF/30 ~ W/ \ 1 — ( ) is the profile width calculated 
assuming (3 = (see Figure 3.17). Using the assumption that only half of the spark is 



active, we can write that 



and the drift velocity 



2/1, 



2 — 

W ^ 2R -^ 



P° 

-'2 



W, 



130 






(3.46) 



Vdr 



2h± 

pK 



(3.47) 



The spark half-width can be calculated using Equation 3.46 as follows 



P° 

h± = Rn 



"'w, 



(3.48) 



/30 



3.4. Drift model 



81 



w = ir 



P2 ^ Qs 



p = -4.3° 




a 



11.6° 



Wf,n = 45° 



Figure 3.17: Top view of the polar cap region in the case of a non-aligned pulsar. Small circles 
represent sparks, the red line corresponds to the line of sight. In general, the observed subpulse 
separation Pg <^06S not describe the actual spark separation Qs (the angular separation between the 
adjacent sparks on the polar cap). In order to calculate the distance between the sparks we use an 
approximation from Equation |3.45[ 



3.4.3 Screening factor 

In our model the drift is caused by a lack of charge in lAR, thus we can write the equation 
for the drift velocity as follows 



v± = Vdr 



cAE X B 
52 ' 



(3.49) 



where AE is the electric field caused by the difference of an actual charge density from 
the Goldreich- Julian co-rotational density. We can a use calculation of the circulation of 
an electric field, Equations 3.20 and 3.26[ to find the dependence of the drift velocity on 
the screening factor: 



-^dr 



EgBr ATTrjhx COS a 



Bl 



P 



(3.50) 



Finally, by using Equations |3.47| and |3.50| we can find the dependence of the screening 
factor on the observed drift parameters 



1 



27rPcj cos a 



(3.51) 



3.4.4 Profile width and subpulse separation 

The key parameters in the above calculations are the pulse width W (or W^q), the char- 
acteristic spacing between subpulses P2 5 ^^cl the period at which a pattern of subpulses 
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crosses the pulse window P3. Of these three only P3 is easy to apply, both W and P2 
need serious study before they can be used. 

In general, the profile width depends on the frequency at which we observe the pulsar, 
and most normal pulsars show a systematic increase in pulse width and the separation of 
profile components when observed at lower frequencies. The model known as radius-to- 
frequency mapping explains this effect as a direct consequence of the emission at higher 
frequencies being produced closer to the neutron star surface than at lower frequencies. For 
this reason both the pulse width and the spacing between subpulses should be measured 
at the same frequency. Note that P3 is not affected by this effect since its determination 
involves analyses of many pulses and does not depend on the pulse width. The observed 
pulse width W, measured in longitude of rotation, can be calculated by applying simple 



spherical geometry (Gil et al. , 1984): 



4 sin a sin [a + p) 



In the above ealeulations we a.e using the W,, ^W/^1- (f ) approximat.o,, where 
Wi3o is the pulse width calculated assuming /3 = 0. In the first approximation we can 
assume that W^o corresponds to the distance 2Rpc at the polar cap which, is valid for 
non-aligned pulsars with a relatively high inclination angle. A more accurate value can be 



calculated using formulas presented in Gil et al. (1984). The running polar coordinates 



along the line of sight trajectory can be expressed in the form 



2 ^ • • ^ , «A , • 2 /5 



p (99) = 2 arcsin y sin — sinasin (a + /3) + sin — , (3.53) 



/ X / sin oj sin a sin (a + /3) , ,^ ^ .^ 

a {if) = arctan , ^ ^, ^ ^Vt . 3.54 

^ cos [a + p) — cos a cos p [ip) 



In numerical calculations of a {(f) it is convenient to use the "atan2" function which 
takes into account the signs of both components and places the angle in the correct quad- 



rant (see the footnote on page 35). Figure 3.18 presents the geometry of the emission 



region for pulsars with available radio observations of the subpulse drift and X-ray obser- 
vations of the hot spot. By knowing the actual polar cap radius Rpc we can determine 
the transverse size of the region responsible for the generation of plasma clouds in lAR 
(the spark half- width). 
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Figure 3.18: Top view of the polar cap region of pulsars with radio drift observations and X-ray 
hot spot radiation. Red lines correspond to the line of sight while green dashed lines correspond to 
the theoretical lines of sight calculated with an assumption that /? = 0°. The geometry of pulsars 
can be found in Table [L3l 

In our model the motion of sparks and the progressively different positions of the 
associated plasma columns are responsible for the observed drift of subpulses. For some 
pulsars it is possible to measure directly the subpulse separation using a single pulse. In 
most calculations it is assumed that the observed subpulses correspond to the adjacent 
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sparks. In general, this is not necessarily true. The distribution of sparks on the polar cap 
is unknown and it is very likely that the line of sight does not cross the adjacent sparks 
but it omits some sparks in between. Therefore, the observed value of P2 should be 
considered rather as an upper limit for spark separation. Furthermore, for many pulsars 
the observed value Pg > ^5 which means that it is not related to any structure at the 



polar cap but that it corresponds to some other periodicity. We can use Equations 3.39 



and 3.51 to calculate the spark half-width as follows 



26.2 



(SiV + 0.3)'PP3V|cosa| 



B 



(3.55) 



14 



Finally, using Equation |3.46| we can determine the predicted value of the subpulse 
separation 



P.° 



26.2 {B\-l + 0.3)^ PP3 v^^ai 



B\aR 



-W, 



14-n-pc 



m- 



(3.56) 



3.4.5 Heating efficiency 

The spin-down energy loss is 



LsD = 3.9 X 10 



31-P-15 



pz 



(3.57) 



We can use Equations |3.26[ |3.48| and |3.51| to calculate the dependence of the acceler- 
ation potential drop on the parameters of drifting subpulses: 



AV ^ 2.824 X 10 



10 



0.5 

P-15 \ 1 



P^ I P3 \Wpo^ 

The power of heating by backstreaming particles can be calculated as follows 



(3.58) 



^heat 



r^ncjc (AVe) ttRI^. 



(3.59) 



The number density of the Goldreich- Julian co-rotational charge can be calculated 
using 



riGJ 



OB, 



(3.60) 



2 Tree 
where Q = 27t/P is an angular velocity, Pg = 6Pd is the surface magnetic field, 

1.5 



b = Rjp/Rl^, Pd = 2.02 X lO^^^PP-isG, and Pdp = ^J2tiR^/ (cP) ^ 1.45 x 10^P-°' 
The Goldreich- Julian density in terms of observed parameters can be written as 



nGj = 2.9x 10^np-3p_i5 



1/2 cos a 



P2 



(3.61) 



Finally, using Equations 3.51 3.58 and 3.61 we can estimate the dependence of the 
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heating power and thus the X-ray luminosity of the hot spot radiation on the parameters 
of radio observations as follows 



exio^M^mrT^^ • (3-62) 



P-15 \ f I P^ 

The heating efficiency by the backstreaming particles can be calculated as 

e.e. = y^ = 0.15f^-^y. (3.63) 

3.4.6 Ion luminosity 

In the PSG-on mode the bulk of energy is transferred to the iron ions which shield the 
acceleration potential drop. Similar as for the backstreaming particles, we can estimate 
the power of ion acceleration as 

Lion = {1-V) riGiC (A^gion) vri?^^, (3.64) 

where gion = 26e = 1.25 x lO"*^ erg^-^cm*^'^ is the ion charge. Using the same approach 
as for electron, we can calculate the dependence of energy transformed to the ions per 
second on the parameters of the radio observations as follows 

Lion = 9.75 X 10^2 (1 - n) ^ (^] cos a. (3.65) 

P-^P-S \W,30j 

It is clearly visible that if the screening factor is low 77 ^ 1, most of the energy in lAR 



is transferred to the iron ions. Using Equations 3.59 and 3.64 we can show that 



Lion 26 (1 - ri) ^ 26 

= ^ — . (3.66) 

Lheat V 7] 

Finally, the ion acceleration efficiency can be calculated as 

e.„ = ^-25^f-^ycosa. (3.67) 

-t>SD -fS V^/30/ 

It may seem that ion luminosity exceeds the spin-down luminosity, but note that both 
P3 > 1 and P2 < WpQ. The predicted values of heating efficiency ^heat and ion acceleration 
efficiency ^ion are presented in the next section. 

3.4.7 Observations 

In this section we confront the values of the subpulse drift and X-ray radiation as es- 
timated by other authors with predicted values estimated using the approach presented 



above. In Table 3.2, alongside our predicted value of P2 we present two other estimates: 



(1) the subpulse separation estimated using the carousel model developed by Ruderman 
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and Sutherland (1975), -P2RSJ (2) the subpulse separation found using the analysis of 



the Longitude-Resolved Fluctuation Spectrum (Backer, 1970) and the integrated Two- 



Dimensional Fluctuation Spectrum (Edwards and Stappers, 2003), performed by Wei 



tevrede et al. (2006), P2°w- ^^ have found that the subpulse separation estimated using 



the fluctuations spectra is overestimated (in most cases -P2°w ^ ^)- -^y definition P2 
should correspond to the structure within a single pulse, thus if the geometry is not ex- 
treme {q > 1°) it should comply with P2 < W. For this specific sample of pulsars -P2°w 
should not be interpreted as the actual subpulse separation. On the other hand, P2 is in 
good agreement with -P2RS- '^^^ predicted values P2 for B0834+06 and B0943+10 sug- 
gest that -P2°RS fo^ those pulsars could be overestimated due to the aliasing phenomenon 
(P2°Rs ~ 2P2°)- For B1929+10 we do not hst P2°rs as its value presented in 



Gil et al. 



(2008) does not comply with the P2RS — ^ condition. We believe that the overestim- 



ated value of P2°RS ^^"^ B1929-I-10 is a result of using the fluctuations spectra presented in 



Weltevrede et al. (2006) to calculate the number of sparks in the carousel model. Note 



that the coincidence of Pg ^^'^ -P2°rs i^ Y^t to be clarified, as in our model there is no 
physical reason for sparks to circulate around the magnetic axis. In fact, the PSG model 
assumes the non-dipolar structure of the magnetic field lines in the gap region and the 
actual position of the polar cap is not necessarily coincident with the global dipole (e.g. 



see Figures 2.6, 2.18, 2.24). 



Table 3.2: Details of a subpulse drift for pulsars with X-ray hot spot radiation. The individual 
columns are as follows: (1) Pulsar name, (2) Predicted characteristic spacing between subpulses in 
the pulse longitude, P2 i (3) Spacing between subpulses, found in the literature, estimated using the 
carousel model, P2RS' ('^) Spacing between subpulses estimated using fluctuations spectra, -P2 W' 
(5) Period at which a pattern of subpulses crosses the pulse window (in units of the pulsar period), 
P3; (6) Number of sparks estimated using the carousel model, N; (7) Profile width at 10%, W; (8) 
Profile width calculated assuming /3 = 0, Wpo; (9) Angular width of the observed region on the polar 
cap Q (see Figure 3.18[); (10) References; (11) Number of the pulsar. 
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Table [3^3] presents the observed and derived parameters of PSG for pulsars with avail- 
able radio and X-ray observations. In the calculations we used the predicted value of 
subpulse separation P2. Note that we consider only pulsars with a visible hot spot com- 
ponent since only for these pulsars we can estimate the size of the polar cap. The low 
value of the estimated screening factor [r] <^ 1) suggests that when the drift is visible, the 
pulsar operates in the PSG-on mode. If the pulsar operates in the PSG-ofF mode, r/ ^ 1, 



the drift velocity is much higher (see Equation 3.50) and the drift phenomenon should be 
more chaotic and much difficult to identify. Observations of PSR 0943+10 show a strong, 
regular subpulse drifting in the radio-bright mode, with only a hint of modulation in the 
radio-quiescent mode. Based on this fact we believe that the two different scenarios of 
the gap breakdown (PSG-on and PSG-off modes) can explain the mode changing. 



Table 3.3: Derived parameters of PSG for pulsars with available radio observations of the subpulse 
drift and X-ray hot spot radiation. The individual columns are as follows: (1) Pulsar name, (2) 
Screening factor, r/; (3) Predicted heating efficiency, ^hoatl (4) Observed bolometric efficiency, ^g^; 
(5) Predicted ion acceleration efficiency, ^.^^^; (6) Surface temperature, T^; (7) Strength of the surface 
magnetic field, B^, (8) Observed polar cap radius, Rpc, (9) Estimated spark half-width, h±; (10) 
Number of the pulsar. Tg, Rpc, b were chosen to fit la uncertainty. Note that in the calculations 
P2 was used. 



Name 
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log Cat 


log^BB 


loge 
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Ts 


Bs 


Rpc 


h± 


No. 
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(x - ray) 
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(iQi^G) 
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0.07 


-4.03 


-3.61 
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B0943+10 
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3.3 


2.5 


17 
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15 


B0950+08 


0.09 
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-2.62 
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2.1 
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0.7 


16 


B1133+16 


0.09 
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-0.81 
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2.9 


22 
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-4.17 


-1.98 
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22 
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Chapter 4 



Cascade simulation 



In this chapter we present the approach of calculating the pair cascades developed by Medin 



and Lai (2010) which has been applied to cases with non-dipolar structure of magnetic 



field. The original approach was adapted to perform full three-dimensional calculations 
and extended with effects that may have a greater importance for non-dipolar configuration 
of surface magnetic fields (e.g. aberration). Additionally, to perform a thorough analysis 
of the Inverse Compton Scattering we present the detailed description of calculating the 



ICS cross section originally developed by Gonthier et al. (2000). 



Following the approach presented by Medin and Lai (2010) we can divide the cascade 
simulation into three parts: 



propagation of the primary particle (including photon emission), 



photon propagation in strong magnetic field (pair production, photon splitting), 
propagation and photon emission of the secondarjQ^articles. 



We use the "co-rotating" frame of reference (the frame which rotates with the star) to 
track both photons and particles. In calculations we consider regions far inside the light 



cylinder. Thus, following Medin and Lai (2010), we ignore any bending of the photon 



path due to rotation of the star. Furthermore, we also ignore effects of general relativity 
on trajectories of photons and particles. 



Figure [471] presents a summary flowchart of the algorithm used to calculate the proper- 
ties of secondary plasma and the spectrum of radiation for a given structure of a neutron's 
star magnetic fleld and gap parameters. 

^In this thesis the term "secondary" refers to any newly created particle except for the primary particles 
accelerated in lAR, e.g. the third generation of electrons and positrons are all considered as "secondary" 
particles. 
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Initial parameters: 
ro = {r,e,4>), 7c 



Calculation of primary particle 
propagation 



7i, ^ = z + 1 




Particle position: 
n = {ri,Oi,(l)i) 



,--'' J sec 

....?.X '"sec \Xsec-) ^sec^ H^s 



t 



Curvature Radiation 



Inverse Compton Scattering 



Spectrum of CR 




ICS photons 



Photon propagation 
(including sphtting effect) 




Photon propagation 
(including splitting effect) 




Secondary particles creation 
Synchrotron Radiation 



^Spectrum of radiation 





/ Secondary plasma \ 




Figure 4.1: Flowchart of algorithm used to calculate a cascade simulation. 
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4.1 Curvature Radiation 



As we have shown in Chapter M an ultrastrong surface magnetic field {Bs > 10^^ G) is 
accompanied by high curvature (curvature radius 3f?6 ~ 0.1 — 10). This suggests that one 
of the important processes of radiation which should be considered is Curvature Radiation 
(CR). 

CR is quite similar to ordinary synchrotron radiation (radiation of ultrarelativistic 
particles in the magnetic field), the only difference being that the radius of circular motion 
(the gyroradius) is in fact the curvature radius of magnetic field lines. Due to beaming 
effects the radiation appears to be concentrated in a narrow set of directions about the 
velocity of the particle. The angular width of the cone of emission is of the order ~ I/7, 



where 7 is a Lorentz factor of an emitting particle (for more details see Rybicki and 



Lightman||1979[ ). 

We track the primary particle above the acceleration zone (the gap region) as it moves 
along the magnetic field line. The length of the step As is chosen so as to achieve sufficient 
accuracy even for large curvature of the magnetic field line, As ^ 0.019f?mm) where 3f?min is 
the minimum radius of curvature. The distribution of CR photon energy can be written 



as (see Equation 14.93 in Jackson, 1998) 



dA^ E 9^3 , r. ^ ^. .. .^ 

K,/3it)dt, (4.1) 



where E = 47re^7^/33fJ is the total energy radiated per revolution, e^^j^ = 37^/ic/(2$R) is 
the characteristic energy of curvature photons, and A's/s is the n = 5/3 Bessel function 
of the second kind. The total energy radiated by a particle after it passes the length As, 



E. , can be written as 

As " 



E^ =E^. (4.2) 



Thus, by using Equations 4.1|and 4.2 we can write the formula for the distribution on 



CR photon energy after a particle passes length As 

dA^ As VSe'^-f 



noo 

/ K,/,it)dt. (4.3) 



de 27r3f^ he /,/, 

It is convenient to divide the spectrum of photon energy into discrete bins. Then, the 
number of photons in each energy bin can be calculated as 

K = / ^de, (4.4) 

where e^ is the lowest energy for the z-th bin and Ae is the energy bin width. Our 
simulation uses 50 bins with an energy range of e^ = 4 x 10^^ keV (soft X-ray) to 
e^g = 4 X 10^ MeV (hard 7-rays). 
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Depending on the photon frequency the polarisation fraction of CR photons is between 



50% and 100% polarised parallel to the magnetic field (see Jackson, 1998, Rybicki and 



Lightman, 1979). Therefore, using similar approach as Medin and Lai (2010) we randomly 
assign the polarisation in the ratio of one photon ±-polarised to every seven || -polarised 
photons, which corresponds to 75% parallel polarisation. 



4.2 Photon propagation 

To explain some of the properties of pulsars and their surroundings (e.g. nebulae radi- 
ation), large magnetospheric plasma densities exceeding the Goldreich- Julian density (see 



Equation 3.7) by many orders of magnitude are required. In order to simulate the process 



of generation of such a dense plasma it is necessary to check the conditions of photon 
decay into electron-positron pairs. 

A photon with energy E^ > 2mc^ and propagating with a nonzero angle \1/ with respect 
to an external magnetic field can be absorbed by the field and, as a result an electron- 
positron pair is created. The concurrent process is photon splitting 7 — )■ 77, which may 
occur even if the photon energy is below the pair creation threshold {E^ < 2m(?). 

In the cascade simulation the photon is emitted (or scattered in the case of ICS) from 
point Pph in a direction tangent to the magnetic field line Asy. The direction vector is 
calculated as the value of the magnetic field at the point of photon creation (see Equations 
2.1, 2.2 and 2.7) normalised so that its length is equal to the desired step Asy = BAs/_B. 



However, the direction of the magnetic field at the point of photon emission does not take 
into account the randomness of the emission direction due to the relativistic beaming 



effect. In Section 4.2.1 we describe a procedure to include the beaming effect in the 
emission process which alters Asy — ?► Asph. Finally, we can write that at the point of 
curvature emission photons are created with energy eph, polarisation || or ±, weighting 
factor N^ (number of photons), and with both optical depths (for pair production r and 
for photon splitting Tgp) set to zero. Since we neglect any banding of the photon path we 
assume that from the point of emission it travels in a straight line. In each following step 
the photon travels a distance Asph. In the co-rotating frame of reference in every step we 
need to take into the account aberration due to pulsar rotation. In order to do so, in every 
step we alter the photon position according to the procedure described in Section |4.2.2[ 



As stated by Medin and Lai (2010) we can calculate the change in the pair production 
optical depth , Ar, and in the photon splitting optical depth Axgp, at the new position 
as: 

Ar ~ Asphi?j|,±, (4.5) 



Arsp - Asph/?;P^, (4.6) 

where -Riu and RT^ are the attenuation coefficients for || or _L polarised photons for pair 
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production and photon splitting, receptively. 

4.2.1 Relativistic beaming (emission direction) 

Due to relativistic beaming the emission direction should be modified by an additional 
emission angle of order ~ I/7. We use the following steps to include the beaming effect 



in our simulation (see Figure 4.2) 



(I) The first step is rotation of the xyz frame of reference in order to align the z-axis 
with As||. In our calculations we used rotation by angle <^y around the y-axis, Ry (<jj^), and 
rotation by angle ^x around the x-axis, Rx {'^x)- The final rotation matrix can be written 
as 



^^x -f^y \Sy) ^x \Sx) 



COS <;y sm <,x sm ^y sm ^y cos <;x 

cos a — sin a 

- sin (^y cos (,y sin c^x cos c^y cos <,x 



(4.7) 



magnetic field line 




Figure 4.2: Relativistic beaming effect of photon emission (for both CR and ICS). In the simulation 
we include the beaming effect by performing three steps: (I) rotation of the xyz frame of reference 
in order to align the z-axis with Asn, (II) transformation of the step vector from a Cartesian to 
a spherical system of coordinates and alteration of the 6 and components with random values 
1/7 cos A and 11, respectively, (III) transformation of the step vector from a spherical to a Cartesian 
system of coordinates and rotation back to the original system of reference. Note that after these 
steps we get a new vector Asph inclined to the primary one, Asph, at an angle ranging from to 
1/7. 



The Euler angles for rotations can be calculated as 
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<^x 



cLLcLIlZ ( fS-y, Sz) 5 



arctan ( — ^ cos ^x) if s^ 7^ 
arctan ( — ^ sin q^) if Sj^ 7^ 

if Sa; = and Sy 



(4.8) 



k 2 



0. 



Note that in order to increase readability, the A symbol and || index were discarded 
(e.g. Sx = As||,^.)- 

(11) The second step is the transformation of the step vector's coordinates in the 
double rotated frame of reference ASpj^ = (s^, Sy, s") to spherical system of coordinates 
and alteration of the 6 and (p components as follows 



J/2 



s"2 



s7 + si' + s'^ , 



arccos 



arctan 



v/^f + 



5y -r 5 



.112 



+ n, 



1 

H — cos A, 

7 



(4.9) 



where A and 11 are random angles between and 27r. The inverse tangent denoted in the 
0-coordinate must be suitably defined by taking into account the correct quadrant (see 



the "atan2" description in the footnote on page 35 ) 



(111) The last step is the transformation of vector components to the Cartesian sys- 
tem of coordinates, Spj^ = [s" sin (sg) cos (s'^) , s'/ sin (sg) sin (s'^) , s^' cos (sg)] and rotation 
back to the original coordinate system Asph = {Ryx)~ Spj^. 

The rotation matrix of this transformation can be written as 



{R, 



yx) 



(Rya 



yx, 



COS ^y — sin qy 

sin q^ sin ^y cos ^^ sin ^j,. cos <^y 
sin <;y cos (^x — sin <;x cos ^x cos ^y 



(4.10) 



4.2.2 Aberration due to pulsar rotation 

Note that in our frame of reference (co-rotating with a star) the path of the photon 



should be curved (see Harding et al. 1978). In the dipolar case the angular deviation 



increases approximately as Sp^^/c = Sph/RLC- When the configuration of magnetic field 
in non-dipolar inclusion of an aberration is even more important. Therefore, the location 
of photon decay should be modified to include the growth of the photon-magnetic field 
intersection angle. 

In our simulation we include the aberration effect by alteration of photon position Pph 



in every step Asph (see Figure 4.3). 
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magnetic field line 




y 



Figure 4.3: Aberration due to pulsar rotation. We use the following procedure to include 
the aberration effect: (I) rotation around the y-axis to align Q. with fx, (II) rotation by angle 
oj = 27rAsbm/ (ci^) around the z-axis (which reflects the pulsar rotation), (III) rotation back to 
the original frame of reference (in which ^ is aligned with the z-axis). 



We use the three-step procedure to alter the photon position. 

(I) Rotation of the xyz frame of reference around the y-axis by angle a, 



p: 



ph 



Ry {a) Pph- Note that here a refers to the inclination of the magnetic axis with 
respect to the rotation axis and we assume that the pulsar's angular velocity vector f2 
lies in the xz-plane. The rotation matrix of this transformation can be written as 



Ry{a) 



cos a sin a 

1 
— sin a cos a 



(4.11) 



(II) After step (I) the 2;-axis is aligned to fi, and in order to include the rotation of the 
pulsar we need to again rotate the frame of reverence by angle u = 27r Asbm/ (cP) around 
the z-axis, Ppj^ = Rz (cu) Pph- We use the following rotation matrix 



RJoo] 



cos uj — sinw 

sinw cosw 

1 



(4.12) 
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(III) The final step is a rotation back to the original frame of reference, 
■Pph — i^y ('^))~ Pph; using the following rotation matrix 



{Ry{a))-' = {Ry{a)y 



cos a — sin a 

1 
sin a cos a 



(4.13) 



4.2.3 Pair production attenuation coefficient 



The pair production attenuation coefficient can be written as (Medin and Lai, 2010) 



R\\ 



R' sin \l/, 



(4.14) 



where R' is the attenuation coefficient in the frame where the photon propagates per- 
pendicular to the local magnetic field (the so-called "perpendicular" frame), \E' is the 
intersection angle between the propagation direction of the photon and the local mag- 
netic field. To increase readability we suppress the subscripts || and ±, but R' has to be 
calculated for both polarisations separately. 



As stated by Medin and Lai (2010) the total attenuation coefficient for pair production 



can be calculated as R' = ^jkR'j,k, where R'^j^ is the attenuation coefficient for the 
process producing an electron in Landau level j and a positron in Landau level k. For 
the electron-positron pair the sum is taken over all possible states (j and k). Note that 
production of electron-positron pairs is symmetric R'j/^ = R'f^j. Thus, to represent the pair 
creation probability in either the {jk) or (kj) state we will use i?';,. For a given Landau 



levels j and k, the pair production threshold condition is (Medin and Lai, 2010) 



E'^ > Ej + E'^., 



(4.15) 



where E' = E^ sin \Ef is the photon energy in the perpendicular frame and 
E'j^ = m(?^\ + 2e^ra is the minimum energy of a particle (electron or positron) in Landau 
Level n. This condition can be written in a dimensionless form as 



X 



1m(? 



E^ 



1 



— \- sin \E' > - 

2mc2 2 



v/l + 2e,j + V^l + 2e,fc 



(4.16) 



where e^ = fieBj {mc) is the cyclotron energy of a particle (electron or positron) in 
magnetic field B in units of mc^. 

The first nonzero pair production attenuation coefficients for both polarisations (_L 



and II) are (Daugherty and Harding, 1983 Medin and Lai, 2010) 



-Rii 



11,00 



2ao x'^Vx'^ - 1 



,-^^y^B 



x> {xo = 1) 



(4.17) 
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^'^'°i " ^ "" 2ao 2x2 



2x2 -e 



E p-2a^'As 



x2- 1-e, + ^ 



-s ^ 4a;2 



x> (xi= (l + v/l + 2e^)/2^ 

(4.18) 



i?; 



01 



2 X 



"s 4x2 



2an 



-2x2/e^ 



X > Xi, 



(4.19) 



^^-l-e. + i^ 



-^11,02 — 2 X 



1 2x2 l + e^-4 



^"""'^ ^/x2-l-2e, + j^ 



ze-^'^'/'is. 



X > ( X2 = ( 1 + Vl + 4e^ ) /2) , 

(4.20) 



R\ 



,02 



2 X 



1 



x2 — e. 



2ao 



x2-l-2e, + ^ 



,e-2^'Afl. 



X > X2, 



(4.21) 



where ao is the Bohr radius (let us note that -R^ go ~ 0)- ^"^ ^^^ above equations for all 
channels except 00 the pair production attenuation coefficients are multiplied by a factor 



of two (see the text above Equation 4.15). 



The pair production optical depth is defined as (Medin and Lai 2010): 



>ph 



T= I R{s)ds= I i?'(s)sin^ds. 
'0 Jo 



'ph 



(4.22) 



We can assume \E' <^ 1, because all high-energy photons (x > 1) will produce pairs 
much earlier than \E' reaches a value near unity. In this limit sin \I' ^ Sph/3^, so the relation 
between x and Sph can be expressed by 



X 






(4.23) 



Equation 4.22 can be rewritten as 



r = Ti + r||,2 + ^±,2 



Tl 



r 

J so 



i?iioods, T\ 



11,2 



rs2 rs2 

/ (i?|j,oo + R\\,oi) ds, rx,2 = / R±flids, (4.24) 

J Sl J SI 



where sq and si are distances which the photon should pass in order to have energy xq 
and Xi, respectively (in the perpendicular frame of reference). Let us note that Sq, Si and 
S2 are of the same order, and if s < sq the attenuation coefficient is zero. 
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The pair production optical depth to reach the second threshold is 

dx 



Sl 



dsR 



ll.oolsj 



so 



2\ 2 px2 

I ^ 

2ao \ E^ J J^^ 



eg f2mc 



x\/x^ — 1 



2x^/e 



(4.25) 



where sq is the distance travelled by the photon to reach the threshold xq = 1, and si is 
the distance travelled by the photon to reach the second threshold xi = (1 + ^/Y+~2e^) /2. 
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Figure 4.4: Panel (a) presents the dependence of the pair production optical depth on the magnetic 
field strength (/3q = B/Bq). Panel (b) presents the dependence of the optical depth on photon energy 
in the perpendicular frame of reference (x = esin^/ {2mc^^). On both panels the photon energy is 
e = 500 MeV, while panel (b) was obtained for magnetic field strength /3q = 1. 



4.2.4 Photon mean free path 



As was shown in the previous section (see Figure 4.4) for strong magnetic fields (e.g. 
(3q > 0.2), Ti, Ty 2, and r_L^2 are much larger than one. Therefore, the pair production 



process takes place according to two scenarios (see also Medin and Lai, 2010). If (3q > 0.2 
pairs are produced by photons almost immediately upon reaching the first threshold, the 
created pairs will be in the low Landau levels (ra < 2). If (5q < 0.2, the photons will travel 
longer distances to be absorbed and the created pairs will be in the higher Landau levels. 
Thus, for strong magnetic fields (/3g > 0.2) the photon mean free path can be approx- 
imated as 



ph ~ So 



^ 



2mc 

E^ 



(4.26) 



while for relatively weak magnetic fields {f3q < 0.2) we can use the asymptotic approxim- 
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ation derived by Erber (1966): 



tph 



GXD I 1 

[e^ / he) mc B sin^i \^x) 



(4.27) 



X 



E^ i?sin\l> 
2mc^ Bn 



(X « !)• 



(4.28) 



4.2.5 Photon-splitting attenuation coefficient 

In our calculations we include photon splitting by following the approach presented by 



Medin and Lai (2010 ). Since only the -L— J-H || process is allowed, for ||-polarised photons the 



photon splitting attenuation coefficient is zero R,F = (Adler 



and Harding 2001). To calculate the splitting attenuation coefficient in the perpendicular 



1971 



Usov 



2002 Baring 



frame for _L-polarised photons we use the formula adopted from the numerical calculation 



of Baring and Harding (1997) 



R 



/sp 



607r2 (315) (^^) 



[e| exp (-0.6x3) + g.OS] [0.25e3 exp (-0.6x3) + 20] ' 



(4.29) 



For photon energies x < 1 this expression underestimates the results of Baring and Hard- 
ing (1997) at /3q = 1 by less than 30%, while at both f3q < 0.5 and f3q ^ 1 the discrepancy 
is less than 10%. As can be seen from Figure 4.5, the attenuation coefficient R'^ drops 



rapidly with the magnetic field strength for /3„ < 1, thus photon splitting is unimportant 



foi' /5g ~ 0-5 (e.g. Baring and Harding 2001 Medin and Lai 2010). 
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Figure 4.5: Dependence of the photon-splitting attenuation coefficient on the energy of the photon 
in the perpendicular frame [x = esin^/ (2?7ic^), vertical axis) and on the strength of the magnetic 
field (/3g = B/Bq, horizontal axis). 
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4.2.6 Pair creation vs photon splitting 



As noted by Medin and Lai (2010), even though the photon sphtting attenuation coef- 



ficient above the first threshold [x > Xq) is much smaller than for pair production (see 



Figure 4.6), in ultrastrong magnetic fields {/3g > 0.5) the ±-polarised photons split before 



reaching the first threshold (see Figure 4.7). On the other hand, the || -polarised photons 
produce pairs in the zeroth Landau level. 




Figure 4.6: Attenuation coefficients of pair production and photon splitting in the perpendicular 
frame of reference. Panel (a) was obtained using photon energy Ej = lO^MeV and magnetic field 
strength B = Bq = 4.414 x 10^^ G (/3g = 1). Panel (b) presents calculations for photon energy 
E^ = lO^MeV and magnetic field strength 5 = 2.5 x lO^-^G (/3q = 5.7). 
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Figure 4.7: Optical depth for pair production and photon splitting for _L-polarised photons. Panel 
(a) presents results for E^ = lO^MeV and B = Bq = 4A x 10^^ G (/3g = 1), while panel (b) was 
obtained using the same photon energy but a stronger magnetic field B = 2.5 x 10^^ G {/3q = 5.7). If 
/3g = 1 the photon creates an electron-positron pair, while in an ultrastrong magnetic field (/3q = 5.7) 
the photon splits before it reaches the first threshold, x = xq. 
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4.2.7 Secondary plasma 



Following the approach presented by Medin and Lai (2010) whenever r > 1 and the 
threshold for pair production is reached {x = xq for || -polarised photons and x = xi 
for ±-polarised photons), the photon is turned into an electron-positron pair. Whereas 



if Tsn > 1 the photon is turned into two photons. Following the results of Baring and 



Harding[ |1997| we assume that the energy of parent photon is equally distributed between 
both newly created photons. A new ||-polarised photon is created with an energy O.Seph 
and weighting factor 2AN^. We assume that the newly created photon travels in the same 
direction as the parent photon, Asph. Note that the photon should split with probability 



1 — e'^, but as shown by Medin and Lai (2010) for cascade results this effect is negligible 



For f^q ~ 0-1 5 the particles are produced in high Landau levels with energy equal to 



half of the photon energy each (see Daugherty and Harding, 1983). In our calculations 



we assume that the newly created particles (electron-positron pairs) travel in the same 
direction as the photon. When Pg > 0.1, on the other hand, we choose the maximum 
allowed values of j and k for the newly created electron and positron. Note that for 
/3g > 0.1 the particles are created in low Landau levels. 



Figure 48 presents the spectrum of Curvature Radiation for a dipolar and non-dipolar 
structure of magnetic field lines. Note the characteristic three peaks in the CR distribution 
for the non-dipolar structure. 
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Figure 4.8: Distribution of CR photons produced by a single prinnary particle for a dipolar (blue 
line) and non-dipolar (red line) structure of the magnetic field. The minimum radius of curvature in 
thee dipolar case is 3f?™ ss 50, while in the non-dipolar case 3*?™ ss 2. In both cases the radiation 
was calculated up to a distance of D = lOOi?, and with an initial Lorentz factor of the primary 
particle 7c = 3.5 x 10'^. 
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Formation of the peaks is caused by the fact that the particle passes regions with three 
different values of curvature: (I) just above the stellar surface, z ^ 1 km, where curvature 
is the highest; (II) at altitudes where the influence of anomalies is comparable with the 
global dipole, z ^ 2.5 km, also with strong curvature; (III) and at altitudes where the 
influence of anomalies is negligible, z > 3.1R, with approximately dipolar curvature (see 



Figure 2.11 ). Hence, the spectrum is a sum of radiation generated in a highly non-dipolar 



magnetic field (high energetic and soft 7-rays) and with radiation at higher altitudes 
where the magnetic field is dipolar (X-rays). The primary particle loses about 63% and 
1% of its initial energy in the non-dipolar and dipolar case, respectively. As can be seen 
from the Figure, to get high emission of CR photons and, thus, a significant density of 
secondary plasma, a non-dipolar structure of the magnetic field is required. 

The high energetic photons produced in a strongly non-dipolar magnetic field will 
either split or create electron-positron pairs. Figure [49] presents the distribution of particle 
energy created by CR photons. Note that for /3g < 0.1 the pairs are created in high 
Landau levels and in order to get the final distribution of secondary plasma energy we 



should consider the loss of particle energy due to Synchrotron Radiation (see Section 4.3). 
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Figure 4.9: Distribution of particle energy created by CR photons calculated for a non-dipolar 
structure of the magnetic field. For this specific magnetic field configuration and initial parameters 
(see the caption of Figure 4.8) the secondary plasma multiplicity is Msec ~ 6 x 10'^. Note that this 
result does not include Synchrotron Radiation and the actual energies of the created pairs are lower 



as they lose their transverse momenta (see Section 4.3) 



4.3 Synchrotron Radiation 

When pairs (electrons and positrons) are created in high Landau Levels they radiate 
away their transverse momentum through Synchrotron Radiation (SR). The secondary 
positron (or electron) is created with energy 'ymc^ and pitch angle \1/, which corresponds 
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to a specific value of Landau Level n. Following Medin and Lai (2010) we choose the 



frame in which the particle has no momentum along the direction of external magnetic 
field. In such a frame of reference the particle propagates in a circular motion transverse 
to the magnetic field (the so-called "circular" frame). The relation of the energy of the 
newly created particle in the circular frame of reference {Ej_ = 'yj_mc'^) with the particle 



energy in the co-rotating frame can be written as (Medin and Lai, 2010) 



7± 



7^ sin^ \& -|- cos^ \E' 



l + 2e^n. 



(4.30) 



The power of synchrotron emission, Psr, can calculated as follows 

2e2 



PsK = ^ (7i - 1) 



2 2 



(4.31) 



In the circular frame E±, is radiated away through synchrotron emission after a particle 
travels a distance 



/ 



SR 



E, 



R 



SR 



7_i_r7ic 



(4.32) 



The particle (electron or positron) mean free path for SR is much shorter than for 



other relevant cascade processes (see Section 4.1 for Curvature Radiation, and Section 



4.4 for ICS). In fact, it is so short that in our calculations we assume that before moving 



from its initial position the particle loses all of its perpendicular momentum p_j_ due to SR 
(see Daugherty and Harding 19821 Medin and Lai 2010). Once the particle reaches the 



ground Landau level {n = 0, p± = 0) its final energy can be calculate as 



711 



(l-Z^^cos^^) 



-1/2 



7/7J 



(4.33) 



here 13 = v/c = a/1 — 1/7^ is the particle velocity in units of speed of light. 



Following the approach presented by Medin and Lai (2010), to simplify the simulation 



we assume that in the circular frame synchrotron photons are emitted isotropically in the 
plane of motion such that there is no perpendicular velocity change of the particle (the 
Lorentz factors 7 and 7_l decrease but 7|| is constant). Thus, the Equation 4.33 remains 



valid until the particle reaches the ground state. In order to simulate the full SR process 
the following procedure was adopted: the particle Lorentz factor in the circular frame 7_l 
drops from its initial value to 7_l = 1 (i.e., n = 0) in a series of steps. Each step entails 
emission of one synchrotron photon, with energy e^ depending on the current value of 7^. 
After the photon emission the energy of the particle is reduced by e^, A7_l = e^/mc^. 
Subsequently, the particle with reduced energy emits a photon with a new value of e^. 
This process continues until the particle is at n = Landau level. Depending on the 
particle's Landau level n, the SR photon energy e^ is chosen in one of three ways. 
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(I) When the particle is created in a high Landau Level {n > 3), we choose the 
energy of the photon randomly but according to a probability based on the asymptotic 



synchrotron spectrum (e.g. Sokolov and Ternov (1968), Harding and Preece (1987)): 



d'^N \/3a/e. 



dtdej 



27r ei_ 



f-F 



_e_L_ 



+ 



e± 



7_L mc^ 



G 



f^s 



where 



;7i^e. 



(4.34) 



(4.35) 



is the characteristic energy of the synchrotron photons, / = 1 — ej_/ {'y±mc^) is the 
fraction of the electron's energy after photon emission, F (x) = x J -ft's/s (t) dt, and 
G (x) = XK2/3 (x). The functions -R's/s and 7^2/3 correspond to modified Bessel functions 



of the second kind. The expression in Equation 4.34 differs from the classical synchrotron 
spectrum (e.g. Rybicki and Lightman, 1979) by a factor of / = 1 — ej_/ {'y±mc'^) which 
appears in several places in Equation 4.34 and by a term with the function G (x). Note 



that in the classical expressions for the total radiation spectra these terms cancel out. 



However, as noted by Medin and Lai (2010) when the quantum effects are considered 



there is asymmetry between the perpendicular and parallel polarisations such that term 
G (x) remain. 

(II) If n = 2, the photon's energy is either that required to lower the particle energy 
to its first excited state (n = 1) or to the ground state (n = 0). The probability of each 
process depends on the local magnetic field strength. We use the simplified prescription 



based on the results of Herold et al. (1982) to calculate the transition rates (see also 



Harding and Preece, 1987). If /3„ < 1 the energy of the photon is set to lower the particle 



energy to the first excited state, e 



mc 



{,/TTWg - v/1 + 2,9g) • If /^g > 1 the photon's 



energy is randomly chosen to be that which is required to lower the particle energy to either 
the first excited state, or the ground state (e_L = mc^ (a/1 + 4/3g — l)), with probability 
50% each. 



(HI) When n 
ground state, e^^ = 
in the ground state, 'j± is recalculated and a new energy of photon is chosen. 



1, the photon's energy is chosen to lower the particle's energy to its 
{^yl + 2/3g — 1) . If after emission of SR photon the particle is not 



mc 



The photon energy in the co-rotating frame can be calculated as 



711' 



(4.36) 



The weighting factor of the emitted photon is the same as the secondary particle that 
emitted it (AA''^ ). In the circular frame the photon is emitted in a random direction 
perpendicular to the magnetic field. Hence, in the co-rotating frame the emission angle 
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can be calculated using Equations 4.30 and 4.33 as follows 



^ 



arcsm 



^-2 — -cosn, 



1 



717| 



(4.37) 



where 11 is a random number from to 2tt. In our simulation we include this emission 



angle by using the same approach as presented in Section 4.2.1 but as the maximum 
value we use \l/ instead of I/7. 

The polarisation fraction of SR photons is the exact opposite of the CR case and it 
ranges from 50% to 100% polarised perpendicular to the magnetic field. Following the 



approach presented by Medin and Lai (2010) in our calculations the photon polarisation 
is randomly assign in the ratio of one || to every seven ± photons, which corresponds to 
a 75% perpendicular polarisation. 

Figure |4.10| presents the distribution of SR produced by a single secondary particle. 
To show the nature of the distribution, a relatively high pitch angle was used. Note that 
when a particle is created at a distance where the magnetic field is relatively weak (e.g. 
Pg = 10~^ for 7 = 10^) then most of the energy is radiated in the range of 1 — 10 keV. 
Thus, we believe that if a strong enough instability forms (that increases the particle's 
pitch angle), the SR process could be responsible for the production of a non-thermal 
component of the X-ray spectrum. 
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Figure 4.10: Distribution of SR produced by a single secondary particle with Lorentz factor 
7 = 10^. We have assunned that the particle was created in a region where the magnetic field 
strength was B = 4.14 x 10^ (/3g = 10"^) and with a pitch angle ^ = 7°. For such a relatively high 
pitch angle the particle loses most of its energy ending with Lorentz factor 7cnd ~ 6. 
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Figure 4.11 presents the final spectrum produced by a single primary particle with 
an initial Lorentz factor of 7c = 3.5 x 10^ for a non-dipolar configuration of the surface 



magnetic field of PSR J0633+1746 (see Section 2.4.2). Due to CR the particle loses about 



/o of its initial energy (Ae = 2.2 x lO^mc^), which is radiated mainly in close vicinity 
of a neutron star, where curvature of the magnetic field is the highest. As the 7-photons 
propagate they will split (only if the magnetic field is strong enough) and eventually most 
photons will be absorbed by the magnetic field - as a result electron-positron pairs emerge. 
These pairs radiate away their transverse momenta through SR, producing mainly X-ray 
photons (at larger distances) and only a few 7-photons (in a strong magnetic field just 
above the stellar surface). Note that at the end (after pair production) only 14% of the 
primary particle's energy (Aeph = 4 x lO^mc^) is converted into photons and the bulk 
of its energy, 54% (Ae 



pairs 



1.8 X lO^mc^), is allocated into secondary plasma. The 
multiplicity for this specific simulation is of the order Msec = 10^. Note that we use 
Msec to describe the multiplicity of secondary plasma in contrast to Mpr which describes 
particle multiplicity in the JAR. 
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Figure 4.11: Final photon distribution produced by a single primary particle. The blue line corres- 
ponds to the initial CR photons distribution for a non-dipolar structure of the magnetic field, while 
the red line presents the final distribution with the inclusion of photon splitting, pair production and 
SR. 



Figure 4.12 presents the distribution of particle energy created by CR photons but 
with the inclusion of SR emission (red line). Note that synchrotron emission both lowers 
the particle energy (after SR maximum at 7 ~ 5 — 8, while without SR at 7 ~ 15 — 20) 
and increases the multiplicity of secondary plasma Msec ~ 10^. 
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Figure 4.12: Distribution of particle energy created by CR photons calculated for a non-dipolar 
structure of the magnetic field. For this specific magnetic field configuration and initial parameters 
(see the caption of Figure 4.8) the secondary plasma multiplicity is Msec ~ 10"^. Note that this result 



does not include Synchrotron Radiation and the actual energies of the created pairs are lower as they 



lose their transverse momenta (see Section 4.3) 



4.4 Inverse Compton Scattering 

The Inverse Compton Scattering (hereafter ICS) process in the neutron star vicinity has 



been studied extensively by Xia (1982); Kardashev et ah (1984); Xia et ah (1985); Daugh- 



erty and Harding (1989); Dermer (1989 1990); Bednarek et aL (1992); Chang (1995) 



Sturner (1995); Zhang and Qiao (1996); Zhang et al. (1997); Zhang and Harding (2000) 



Harding et aL (2002), etc. According to these studies, the ICS process may play a signific- 



ant role in the physics of a neutron star's magnetosphere. Relativistic particles (positrons 
and electrons) can Compton-scatter thermal radiation from the neutron star surface. As a 
particle with a certain relativistic velocity scatters the thermal photons with a blackbody 
distribution, it will produce radiation in quite a wide energy range. However, we can dis- 
tinguish two characteristic frequencies of upscattered photons: one is the frequency due 
to resonant scattering, another is the range of frequencies contributed by the scattering 
of photons with frequencies around the "thermal-peak". The Resonant Inverse Compton 
Scattering (RICS) corresponds to a scenario when the scattering cross section is largest. 
On the other hand. Thermal-peak Inverse Compton Scattering (TICS) corresponds to 
interactions with photons with the maximum number density. These two modes are very 
different when it comes to the nature of the process. The photons' energy in RICS de- 
pends on the strength of the magnetic field, thus at low altitudes (where the field is very 
strong), it can power pair cascades, while TICS can be responsible for magnetospheric 
radiation at much higher altitudes. Note that for some specific combinations of magnetic 
strength and distribution of background photons, RICS and TICS are indistinguishable 
as the resonance frequency falls into the thermal peak range. 
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4.4.1 The cross section of ICS 



Due to the rapid time scale for synchrotron emission (see section 4.3), a particle in an 



excited Landau level almost instantaneously de-excites to the ground level. The particle 
motion is therefore strongly confined to the magnetic field direction. In our calculations we 



consider the geometry illustrated in Figure 4.13 In the observer's frame of reference (OF) 



a particle with Lorentz factor 7 travelling along the magnetic field line scatters a photon. 
Let ip = arccos/i be the angle between the magnetic field line (particle propagation) and 
the direction of photon propagation in OF and ip' = arccos fi' in the particle rest frame 
(PRF). The energy of the photon in PRF is given by 



e' = 7e(l-/3/x) 



(4.38) 



After scattering, the photon energy is denoted by e'^ in PRF and e^ in OF. The angle 
between the direction of propagation of the scattered photon and B (which describes the 
direction of particle propagation) is denoted by ips = arccos /i^ in OF and ip'^ = arccos /i'^,, 
where /i', = (/i, - /3) / (1 - /3/i,) in PRF (|Dermer| [l990|. 



Observer's Frame (OF) 




Particle Rest Frame (PRF) 




Figure 4.13: Reproduction of the Figure from |Dermer ( |1990 ). Geometry of the ICS event in the 
observer's frame (left) and the particle rest frame (right). A particle with Lorentz factor 7, beamed 
along the direction of the magnetic field, scatters a photon with energy e directed at angle i/j with 
respect to the magnetic field line. After scattering, the energy and angle of the photon are denoted 
by es and ips, respectively. Quantities in the particle rest frame are denoted by a prime. 



4.4.1.1 ICS cross section in the Thomson regime 

Restriction to the Thomson regime requires that 7e (1 — ;u) ^ 1. In the particle rest 
frame, the angle ip' = arcsin{7~^ [sinip/ (1 — fScosip)]}, and when 7^1, |/i'| — )■ 1. In 
the Thomson regime the only important Compton scattering process involves transitions 



between ground-state Landau levels. Daugherty and Harding (1989) and Dermer (1989) 



calculated the differential cross section (after summing over polarisation modes and in- 
tegrating over azimuth) for a photon scattered from ip' = 0° into angle ip'^ = arccos ^'^ as 
follows 
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da' 3(7^ 



d/i'. 



16 



(1 + /^") 



,/2 



.12 



+ 



L(e' + eJ2 (e'-eJ2 + (r/2)2j 



(4.39) 



where T = Aaje /3 is the resonant width (Daugherty and Ventura 



1978 



XiaetaL 1985) 



0"^ is the Thomson cross section, and af = e^ /he is the fine-structure constant. In the 
Thomson hmit e' <^ 1, and thus the scattered photon energy in PRF can be approximated 
as 

(4.40) 



e;^e' + e'^(/.'-^;)V2^e'. 



Equations 4.39 and 4.40 show that a differential magnetic Compton scattering cross 
section when 7 ^ 1 is similar in form to a nonmagnetic Thomson cross section. The 
important difference is that the magnitude of the cross section is enhanced when e' ap- 
proaches e^ and is depressed at energies e' < e^. The total cross section for magnetic 



Compton scattering, obtained by integrating Equation 4.39 over /i'^, was calculated by 



Dermer (1989); Zhang et al. (1997) and is given by 



a 
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IC 
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u 



u 



where a,, 



{u + l)"^ (M-l)2 + a 
cr^, a^ is the Thomson cross section, u = e'/e^, a 



(4.41) 



ftt/e^. 



4.4.1.2 ICS cross section in the Klein-Nishina regime 

The Klein-Nishina regime includes quantum effects due to the relativistic nature of scatter- 
ing, and it requires that 7e (1 — /x) > 1. The principal effect is to reduce the cross section 
from its classical value as the photon energy in PRF becomes large. In the Klein-Nishina 
regime instead of a^^^ = a^ we can use the following relationship 



3 r 1 + e' 



2e' (1 + e') 
l + 2e' 



In (1 + 2e') 



+ — ln(l + 2e') 



l + 3e' 
;i + 2eO' 



(4.42) 
In an extreme relativistic regime e' ^ 1 the Klein-Nishina formula can be simplified 



to 



cr,. 



/-I 



In (2e') + 



(4.43) 



The above formula clearly shows that Inverse Compton Scattering is less efficient for 
photons with energy in PRF significantly exceeding particle rest energy. 



4.4.1.3 QED Compton Scattering cross section 

Previous studies on upscattering and energy loss by relativistic particles have used the 
non-relativistic, magnetic Thomson cross section for resonant scattering or the Klein- 



Nishina cross section for thermal-peak scattering. As noted by Gonthier et al. (2000) 



this approach does not account for the relativistic quantum effects of strong magnetic 
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fields {B > 10^^ G). When the photon energy exceeds mc^ in the particle rest frame, 
the strong magnetic field significantly lowers the Compton scattering cross section below 



and at the resonance. Gonthier et al. (2000) developed expressions for the scattering of 



ultrarelativistic electrons with 7^1 moving parallel to the magnetic field. Because of the 
large Lorentz Factor of particle 7, the photon incident angle ip gets Lorentz concentrated 
to tp' ~ V'/27 ~ 0° in the PRF. The differential cross section in the rest frame of the 
particle can be written as 
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(4.46) 



The imaginary terms and the 0' dependence are isolated in the polarisation components 
and in the phase exponentials, leading to elementary integrations over the azimuthal angle. 



' (Gonthier et al. 2000) 



The differential cross section depends on the final Landau state /, thus a sum must be 
calculated over all the contributing Landau states. The energy of the scattered photon is 



given by (Gonthier et al. , 2000) 



2 U' - /e. 



1 ■> 
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(4.47) 



1 + e' (1 - cosV^'J + {[1 + e' (1 - cos7/;'J]^ - 2 (e' - It^) sin^ V^'J- 
where / is the final Landau level of the scattered particle. Each final state has an energy 
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threshold of le^ , thus the maximum contributing Landau state /max can be expressed as: 
e'/e^ — 1 < /max < e'/^s- To obtain the energy-dependent cross section, the Romberg's 
method can be used to numerically integrate the differential cross section over ip'^. For this 
particular case (scattering of relativistic particles) there is only one resonance appearing 
at the fundamental cyclotron frequency e^ = f3q = eB / (mc). 

The values of B and C can be expressed as: 
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Although the expressions presented above describe the exact cross section for ICS in 
strong magnetic fields, due to their complexity their usage in cascade simulation is limited. 
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4.4.1.4 Approximate cross section (final states 1=0) 

An approximation to the exact / = differential cross section can be given by assuming 
that the scattering is significantly below the resonance, where e' < e^ and also e' < 1. 



Gonthier et al. (2000) showed that by keeping only terms to first order in e' and e'^ in the 



region of validity, it agrees very well with the exact / = cross section. The approximation 
overestimates the exact / = cross section above the region of validity e' > e^. However, 
the approximation is close to the total cross section for both energy regions (e' < e^ and 



e'>e. 



even for high magnetic field strengths (see Figure 4.14) 



According to Gonthier et al. (2000), the polarisation-dependent and averaged, approx- 



imate cross section can be calculated as: 
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Figure 4.14 presents the total approximate cross section of Compton scattering, the 
exact QED cross section (summed over all contributing final electron/positron Landau 
states) and the exact cross section for final Landau state Z = as a function of energy 
of the incident photon in PRF (in units of cyclotron energy, e'/e^). As mentioned above, 
the approximation is valid in the region below the resonance, e' < e^. Although the 
approximation overestimates the cross section for / = final Landau state in the regime 
of high energetic photons (e' > e^), it can be used in this regime as the approximation of 
the total cross section. In our simulation we use this approach to calculate the total ICS 
cross section in both regimes, e' < e^ and e' > e^. Calculation of the cross section for the 
resonance frequency (e' = e^) is presented in the next section. 
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Figure 4.14: Total cross section of Compton scattering (in Thomson units) as a function of an 
incident photon energy in PRF (in units of the cyclotron energy) calculated for a magnetic field 
strength B\i^ = 3.5. The exact QED scattering cross section, summed over all contributing final 
electron/positron Landau states, is indicated as the red dotted curve. The cross section for final 
Landau states / = is plotted as a blue dashed line. 



4.4.2 Resonant Compton Scattering 

This section describes an approach used to calculate the RICS cross section for ultrastrong 
magnetic fields {B > 10^^ G). For weaker fields the calculations are much simpler and 



resonance is already included in Equation 4.41 



The trend as Pq increases is for the magnitude of the cross section to drop at all 
energies. For weaker magnetic fields (/3g < 1) the width of the resonance increases with 
increasing f3g, but for (3g > 1 this width actually declines. Since the resonance is formally 



divergent, the common practice (see Xia et al. 1985 Latal 1986 Daugherty and Harding 



1989[ |Dermer[ |1990t [Harding and Daugherty[ [19911 [Baring et al.[|20d5| [Harding and Lai 
2006 Baring and Harding, 2007, 2008) is to truncate it at e' = e^ by introducing a finite 
width r. The procedure is to replace the resonant (e' — e 



4.44 



and 



4.51) by [(e' 



denominator (see Equations 
+ r^/4]. In the /3g <^ 1 regime, the cyclotron decay width 
assumes the well-known result T ^ 4a/e^/3 in dimensionless units. When Pq ^ 1, 
quantum and recoil effects generate T ^ o^f^g (1 — 1/e) where e is Euler's number (e.g. 



see 



Baring et al. , 2011 ). These widths lead to areas under the resonance being independent 



of e^ in the magnetic Thomson regime of /3g ^ 1 and scaling as e^^^ when /3g ^ 1. These 
results can be deduced using the / = approximation derived in Equation 4.51 By using 
this approach the averaged, approximate cross section can be written as 



avg 






[9{e') + f{e')~{l + 2e')h{e')] 



(e' - e J^ + rV4 + (e' + e J 



(4.55) 
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The common practice to calculate a resonant cross section in an ultrastrong magnetic 



fields is to use the Dirac delta function as follows (e.g. Medin and Lai, 2010) 



a — 

^res 



27r^ 



,e^h 



5U[ 



mc 



(4.56) 



This simplified approach, however, does not include scatterings of photons whose en- 
ergy in a particle rest frame is not equal but very close to the resonance frequency. The 
relativistic quantum effects of strong magnetic fields that are included in the approximate 
solution increase the cross section, and thus the efficiency of the ICS process in previous 
estimates could be underestimated. 



According to Medin and Lai (2010) in ultrastrong magnetic fields the ICS polarisation 



fraction is about 50% (approximately 50% of the photons are slightly above resonance and 
50% are slightly below). Therefore, the polarisation of ICS photons is randomly assigned 
in the ratio of one _L (perpendicular to the field) to every || photon. 



4.4.3 Particle mean free path 

For the ICS process the calculation of the particle mean free path /ics is not as simple as 
that of the CR process. Although we can define lies in the same way as we defined Iqr, 
it is difficult to estimate a characteristic frequency of emitted photons. We have to take 
into account photons of various frequencies with various incident angles. An estimation 



of the mean free path of a positron (or electron) to produce a photon is (Xia et al. , 1985) 



/ 



ICS 



Ml poo 



MO "'0 



'1 - /3/i) a' (e, fi) Uph (e) ded/i 



(4.57) 



where (as before) /3 = v/c is the velocity in terms of speed of light, riph represents the 
photon number density distribution of semi-isotropic blackbody radiation (see Equation 



3.33). Here a' is the average cross section of scattering in the particle rest frame (see 



Equation 4.55). We should expect two modes of the ICS process, i.e. Resonant ICS and 



Thermal-peak ICS. 



4.4.3.1 Resonant ICS 



The RICS takes place if the photon frequency in the particle rest frame is equal to the 



cyclotron electron frequency. Using Equation 4.38 we can write that the incident photon 
energy is e = eg/[7(l — f3fi)]. For altitudes of the same order as the polar cap size 
we use /io = 1, yUi = as incident angle hmits for outflowing particles, and /io = 0, 
yUi = — 1 as incident angle limits for backfiowing particles. Thus, for outflowing particles 
the electron/positron mean free path above a polar cap for the RICS process is 
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(4.58) 



where limits of inteeration, e™'" and g""''' are chosen to cover the resonance. In our 

*-' ' res res ' 

simulation we use such limits to include the region where the integrated function decreases 
up to about two orders of magnitude from its maximum: 



min/max 



e,±yilT 



Here F is the finite width introduced in Section [4.4.2| to describe the decay of an excited 
intermediate particle state. 

Figure |4.15 presents the dependence of the integrand from Equation 4.58 on the in- 
cident photon energy for a given incident angle. The maximum of the integrand shows a 
significant decline for stronger magnetic fields. This is due to both the drop of the cross 



section at all energies with an increasing magnetic field (see Section 4.4.2) and due to the 
fact that for this specific incident angle resonance is in a different range of photon energy. 
In stronger magnetic fields resonance occurs not only for higher energetic photons but 



also the width of the resonance is wider (see the right panel of Figure 4.15). 
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Figure 4.15: Dependence of the integrand from Equation 4.58 on the energy of the incident photon. 
Both panels were calculated for surface temperature T = 3 x 10® K, cosine of the incident angle 
fjL = 0.1 and Lorentz factor of particle 7 = 10'^. The left panel corresponds to resonance in magnetic 
field B = lO^'' G, while the right panel was obtained using 5 = 3 x lO^"^ G. 



Note that both plots do not include the dependence of the photon density on distance 
from the stellar surface. Depending on whether the radiation originates from the whole 
stellar surface or from the polar cap only, the dependence of the photon number density 



on the height above the surface can differ significantly (see Section 4.4.4). 
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4.4.3.2 Thermal-peak ICS 

TICS includes all scattering processes of photons with frequencies around the maximum 
of the thermal spectrum. In our simulation we adopt e™'" ~ O.OSe^^, and e™'"' ~ 2e^^ where 
e^j^ = l.^lkT I {mc^) is the energy, in units of mc^, at which blackbody radiation with 
temperature T has the largest photon number density. The electron/positron mean free 
path for the TICS process can be calculated as 



TICS 



Ml 



MO 



-\ -1 



f3fi) a (e, yu) riph (e) ded/i 



(4.60) 



Figure 4.16 presents the dependence of the integrand from Equation 4.60 on photon 
energy for two different incident angles of background photons. As the number density 
depends exponentially on the photon energy, TICS is important only for small incident 
angles (/i ~ 1). Note that for some specific combination of magnetic field strength, 
the Lorentz factor of the primary particle and the incident angle of background photons 
the resonance is in region of thermal peak. In such a case the resonant component is 
dominating (a much higher cross section) and the particle mean free path should be 
calculated using the approach described in the previous section. 




Figure 4.16: Comparison of the integrand from Equation 4.60 with photon number density. The 
bottom panels present the dependence of the photon number density on photon energy in OF. The 
red dashed lines correspond to limits used to calculate the particle mean free path for TICS. The top 
panels present the dependence of the integrand on photon energy in PRF. Both panels were obtained 
using surface temperature T = 3 x 10^ K, Lorentz factor of the particle 7 = 10^ and magnetic field 
strength B = 10^^ G. The cosine of the incident angle, ^ = 0.975 and ^ = 0.96, was used for the 
left and right panel, respectively. 
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4.4.3.3 Calculation results 

For ultrastrong magnetic fields quite a wide range of the particle Lorentz factor falls into 



the peak of background photons (see Figure 4.17). In such a case RICS is enhanced by 



the fact that it involves photons with very high density. Furthermore, the RICS process 



for such particles is indistinguishable from the TICS (see Figure 4.16). For particles with 
Lorentz Factor 7 > 10^, the dominant process of radiation is CR. The exact value of 
this limit depends on conditions such as: density of background photons, incident angles 
between particles and photons, and curvature of magnetic field lines (1/3^?). 
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Figure 4.17: Dependence of a particle mean free path on its Lorentz factor for three different 
processes: CR, RICS and TICS. The calculations were performed for magnetic field strength B\/^ = 2, 
radius of curvature of magnetic field lines Ke = 1 (for the CR process) and hot spot temperature 
Tg = 3 (for RICS and TICS). Both RICS and TICS were calculated for a full range of incident angles 
(/io = 0, /xi = 1). Note that for a Lorentz factor in the range of 7 ss 2 x 10'^ — 10^ the particle 
mean free paths of RICS and TICS are equal as the resonance falls into the peak of the background 
photons. 



Figure 4.18 presents the dependence of a particle mean free path on the magnetic 
field strength and the particle Lorentz factor for RICS. The minimum of the mean free 
path for relatively weak magnetic fields (i?i4 = 0.5) is for particles with Lorentz factor 
7 Ri 2 X 10^, while for relatively stronger magnetic fields {Bn = 3.5) the RICS is most 
efficient for particles with energy an order of magnitude larger (7 ^ 2 x 10^). This is 
a natural consequence of the fact that resonance takes place when the photon energy in 
PRF is equal to the electron cyclotron energy, which in stronger fields is higher. As can 
be seen from the Figure, the particle mean free paths for RICS in stronger magnetic fields 
increase. This is due to the decreasing resonant cross section with increasing magnetic 



field strength (see Figure 4.15). Note, however, that this behaviour does not include the 



fact that photon density in regions with weaker magnetic fields is considerably smaller. 
In fact, the results of the cascade simulation presented in Chapter [5] show that RICS is 
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efficient only in the immediate vicinity of a neutron star since photon density at higher 
altitudes drops rapidly. 
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Figure 4.18: Dependence of a particle mean free path on magnetic field strength (-B14) and the 
Lorentz factor of a particle (7) for the RICS process. The particle mean free path was calculated for 
semi-isotropic blackbody radiation {fiQ = 0, /.ii = 1) with temperature Tg = 2.5. 



4.4.4 Background photons 

4.4.4.1 Photon density 

One of the main parameters affecting ICS above the stellar surface is photon density. 
The initial photon density (at altitude z = 0) highly depends on the temperature of the 
radiating surface. As shown in Chapter [I] (e.g. see Table 1.4), the entire surface has the 
lowest temperature (Tg < 0.8), thus the initial photon density is up to about two orders 
of magnitude lower than warm spot radiation (Tg < 3) and up to about three orders 
magnitude lower than hot spot radiation (Tg < 5). However, the density of the photons 
strongly depends on the distance from the source of radiation (especially for the hot spot). 
Therefore, we used the simplified method presented in Figure |4.19| to calculate photon 
density at a given point L = {r, 6, (p). Then the relative density of photons originating 
from the entire surface can be calculated as 



Ust (e, Tst, L) 



sm 



Ae., 



(4.61) 



no (e, Tst) '" V 2 

where Ustfl (e, Tst) is the density of photons with energy e at the stellar surface with 
temperature Tst, and A^st is the angular diameter of the star at a distance from the star 
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centre r. 

Likewise, we can write a formula for the relative density of photons originating from 
a spot (warm or hot) as 



risp n (e, Tsp) V 2 ' ' 



(4.62) 

""spjU V*-' -^sp^ V ■^ / 

where Ugp^ (e, Tgp) is the density of photons with energy e at the spot surface (either hot 
or warm) with temperature Tsp. The angular diameter of the spot can be calculated as 



A9 = arccos 



rf + rj- 4i?sp 
2rir2 



(4.63) 



here i?sp is the spot radius and ri, r2 are the distances to the outer edges of the spot (see 



Figure 4.19). 
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Figure 4.19: Simplified method used for calculation of a photon density originating from an entire 
stellar surface (blue lines) and from a hot/warm spot (red lines). Here i?sp is a spot radius (either hot 
or warm). Let us note that the simplified method is valid for the entire surface component regardless 
of the component of location L, while for the spot component it can be used only for small values 
of (\). In a more general case the spot should be projected on the surface perpendicular to the radius 
vector r and passing through point L. 



Figure 4.20 presents the dependence of the relative photon density [n (2;) /rio) on 



the distance from the stellar surface. Due to the small size of a polar cap (hot spot, 
-Rhs = 50 m) the density of the photons drops rapidly and already at a distance of about 
z = 150 m it is one order of magnitude lower than at the polar cap surface. On the other 
hand, for a larger size of the warm spot (-Rhs = 1 kni) the photon density is reduced by 
an order of magnitude at a distance of about z = 3 km. From Equation 4.61 it can easily 
be seen that the photon density of radiation from the entire stellar surface decreases by 
an order of magnitude at a distance of about z ~ 3-R ~ 30 km. 
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Figure 4.20: Dependence of the relative photon density on the distance from the stellar surface 
for three different thermal components (the entire stellar surface, the warm spot and the hot spot). 
The following parameters were used for the calculations: star radius R = 10km, warm spot radius 
i?ws = 1km and hot spot radius i^hs = 50 m. 



The very small size of the polar cap also has an additional implication to the back- 
ground photons' density. Namely, the density of the background photons just above the 
polar cap highly depends not only on the distance from the surface, but also on the 
position relative to the cap centre. 



Figure 4.21 presents the dependence of the relative photon density originating from a 



polar cap (the hot spot) on the distance from the stellar surface for three different starting 
points on the polar cap. The distance was calculated for points which follow the magnetic 
field structure of PSR B0656+14. Note that for the extreme magnetic line (which starts at 
the cap edge) already at a distance of about ^2 ~ 5 m the photon density decreases twice, 
while for the central (6*0) and middle line (6'i) the distances are respectively 2:0 ~ 45 m and 
Zi ~ 30 m. This result is important as the background photon density directly translates to 



the particle mean free path in ICS (see Section 4.4.3). This means that for ICS-dominated 
gaps the sparks' height will vary depending on their location. The breakdown of the gap 
(spark) in the central region of a polar cap is easier to develop as the particle mean free 
path is lower, and eventually it will result in lower heights of the central sparks. This 
will influence the properties of plasma produced in the central region of open magnetic 
field lines, and depending on the conditions may result in the formation of plasma either 
suitable to produce radio emission (core emission) or unsuitable to produce radio emission 
(conal emission but with the line of sight crossing the centre of the beam). 

To find the dominant component of thermal radiation at a given altitude we need 
to take into account the initial flux of radiation and how it changes with the distance. 



Below we present the calculations of a radiation flux (Figure 4.22) for PSR B0656+14 



The parameters of an entire surface and warm spot components are in agreement with 
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the observations (see Table 1.4), while the hot spot component was calculated using 
parameters derived from the modelling of a non-dipolar structure of the magnetic field 
(see Chapter [2]). 
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Figure 4.21: Dependence of the relative photon density on the distance from the stellar surface 
for a hot spot component of PSR B0656+14. The relative photon density was calculated for three 
different starting positions: Oq (central), 61 (at the half distance to the edge), and 02 (the cap 
edge). The altitude (z) was calculated for points which follow the magnetic field structure of 
PSR B0656+14. 
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Figure 4.22: Dependence of the radiation flux for three different components (the entire stel- 
lar surface, the warm spot and the hot spot) on the distance from the stellar surface for PSR 
60656+14. The following parameters were used for the calculations: entire stellar surface radiation, 
Tst = 0.7 MK, Est = 20 km; warm spot, Tws = 1.2 MK, R^s = 1.8 km; and hot spot, Ths = 2.9 MK, 
i?hs = 50 m. 



Already at a distance of 240 m the flux of the warm spot radiation becomes higher 
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than the flux of the hot spot radiation. Furthermore, aheady at a height of 750 ni flux 
the radiation originating from the polar cap (hot spot) becomes lower than the flux of 
radiation from the entire stellar surface. With an increasing distance the flux of the warm 
spot decreases faster than the flux of the entire surface radiation and at a distance of 6.3 km 
the thermal radiation from the entire stellar surface becomes the dominant component of 
the background photons. 

The results may suggest that up to a height of about 240 m (for PSR B0656+14) the 
hot spot radiation should be the main source of the background photons involved in ICS. 
However, the actual height is smaller as the results do not include the efficiency of ICS, 
which also depends on the incident angle between the photons and the particles (see the 
next Section). 



4.4.4.2 Photon incident angles 

Another parameter that significantly affects the ICS is the incident angle between the 
background photons and the relativistic particles. Especially for Resonant Inverse Compton 



Scattering is the incident angle of great importance. Figure 4.23 presents the dependence 
of a particle mean free path for ICS on a maximum value of the incident angle ipcrit- If 
incident angles are low, the resonance is outside of the photon spectrum and results in 
very high values of particle mean free paths. The lower the energy of the particle (lower 
Lorentz factor), the incident angles should be larger to ensure that the resonance falls 
into an energy range with high photon density. 
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Figure 4.23: Dependence of the particle mean free path on the maximum value of the incident 
angle V'crit- The particle mean free path Ip was calculated for magnetic field strength B = 10^^ G 
assuming background blackbody radiation with a temperature T = 3MK. Two different particle 
Lorentz factors were used for the calculations: 7 = 10'^ (dashed lines) and 7 = 10'* (solid lines). 
The red lines correspond to Resonant Inverse Compton Scattering, while the blue lines correspond 
to Thermal-peak Inverse Compton Scattering. 
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TICS for a given magnetic field strength and the Lorentz factor of particles is not 
significant (high particle mean free paths) unless the angles of the incident photons are 
high enough. Note the characteristic drop of the particle mean free path for TICS at 
'j/'crit ~ 20° (for 7 = 10"^) and ipcvit ~ 75° (for 7 = 10'^). For such high incident angles the 
resonance takes place at the thermal peak of the background photons. Therefore, TICS 
and RICS are indistinguishable, which results in an almost equal particle mean free path 



(see the text above Figure 4.16 for more details). 

Due to the very small size of the polar cap the influence of the hot spot component 
will by lower not only because of the change of photon density, but also because of the 



rapid change of the incident angle between the photons and particles. Figure 4.24 presents 
the dependence of the maximum incident angle on the altitude above the stellar surface 
for three thermal components (the entire surface, the warm spot and the hot spot). As 
follows from the Figure, already at an altitude oi z ^ 90 m does the maximum value 
of the incident angle between the photons from the hot spot and the particles drop to 
■^crit = 30°, which significantly lowers the efficiency of ICS for this source of background 



photons (see Figure 4.23). Since the size of the warm spot component is larger, the warm 
spot radiation will be significant for up to higher altitudes, but already at a distance of 
z ~ 1.5 km the maximum value of the incident angle also drops to ^/'crit = 30°. 
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Figure 4.24: Dependence of the maximum incident angle on the altitude above the stellar surface 
for three thermal components (the entire surface, the warm spot and the hot spot radiation). 



Note that in the Figure we have calculated the maximum value of the intersection 
angle at altitudes which correspond to radial progression from the stellar surface. In fact, 
the actual maximum value of the incident angle also depends on the structure of the 



magnetic field. Figure 4.25 presents the actual maximum value of the incident angle of 
photons originating from the hot spot for three different magnetic field lines calculated for 
PSR B0656+14. The actual values of the maximum incident angle just above the surface 
exceed 90°, but its rapid decline (especially for extreme lines) causes the radiation of the 
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hot spot component to become insignificant for ICS at relatively low altitudes z ^ 20 m.. 
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Figure 4.25: Dependence of the maximum incident angle on the altitude above the stellar surface 
for the hot spot component of PSR B0656+14. The maximum incident angle was calculated for 
three different starting positions: ^o (central), 0i (at the half distance to the edge), and 02 (the cap 
edge). 



Both the decrease of photon density and the decrease of the maximum inclination angle 
cause the parameters of plasma produced by RICS to highly depend on the properties 
(size and temperature) of the background photons source. The hot spot component will 
be the dominant source of background photons for ICS in the gap region {z < 20 m), 
while the radiation of the warm spot and the entire surface will be the main source of the 
background photons for ICS at higher altitudes. 
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Physics of pulsar radiation 



5.1 Inner Acceleration Region 

5.1.1 Gamma- ray emission 

In our model most of the 7-photons are produced in the Inner Acceleration Region or 
in close vicinity of a neutron star. Due to an ultrastrong surface magnetic field, the 
most energetic 7-photons are produced by Inverse Compton Scattering in the PSG-on 
mode. If a pulsar is in the PSG-off mode, Curvature Radiation produces fewer energetic 
photons than ICS in the PSG-on mode. Photons produced in lAR (both the ICS and 
CR) are absorbed by strong magnetic fields creating positron-electron plasma in the gap 
region, thereby enhancing a cascade, or just above the gap enhancing a secondary plasma 
population. The absorption of 7-photons in close vicinity of NS makes it impossible to 
directly observe the radiation produced in lAR. However, a characteristic of this emission 
defines the parameters of the gap (e.g. multiplicity in the gap region, gap height, etc.), 
and thus the parameters of secondary plasma. 

5.1.1.1 PSG-off mode 

In general, the existence of high potential in lAR (e.g. wide sparks or r^ ^ 1) results in 
solutions for which CR is responsible for the emission of 7-photons. The energy of such 
radiation depends on the Lorentz factor of primary particles and curvature of the magnetic 



field lines. Figures 5.1 and 5.2 present the histogram of photons produced in lAR by CR 



for PSR B0628-28 and Geminga, respectively. The curvature in lAR of Geminga is lower 



(SIq ~ 2.1, see Section 2.4.2), thus the primary particle should be accelerated to higher 
energies in order to produce the required number of photons in the gap region. Eventually 
the higher Lorentz factor of primary particles will result in the emission of 7-photons with 
energy up to 10 GeV for Geminga. On the other hand, the curvature magnetic lines for 



PSR B0628-28 (^Rg = 0.6, see Section 2.4.1) is higher, which reduces the photon mean 



free path and it is possible to produce the required number of photons in the gap region 
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A^^^ for lower the Lorentz factor of primary particles. 

In CR-dominated gaps we can distinguish three types of photons: (I) radiation with 
energy below 1 MeV which is unaffected by the magnetic field (except the splitting) and 
can be detected by a distant observer, (II) soft 7-ray photons which create pairs above 
ZPF, (III) and high energetic 7-photons responsible for pair production below ZPF. In 
an ultrastrong magnetic field the photons from the third group will produce particles just 
after reaching the first threshold. Due to the fact that most CR photons are || -polarised, 
photon splitting is insignificant in cascade pair production in the PSG-off mode. 
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Figure 5.1: Distribution of photons produced in lAR by a single particle for PSR B0628-28. In the 



calculations we used parameters of the gap in the PSG-off mode as presented in Table 3.1 We also 



assumed a linear change in the acceleration electric field (see Equation 3.27) 
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Figure 5.2: Distribution of photons produced in lAR by a single particle for PSR J0633-1746. In 



the calculations we used parameters of the gap in the PSG-off mode as presented in Table 3.1 
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5.1.1.2 PSG-on mode 

When the acceleration potential is low enough (narrow sparks with r] < 1) to satisfy 
the condition for effective ICS {lies ^ Lcc), the gap will operate in the PSG-on mode. 
The energy of ICS radiation in the gap region (RICS) depends on the Lorentz factor of 
primary particles and the strength of magnetic field. In an ultrastrong magnetic field of 
lAR implied by the PSG model, the primary particle loses most of its energy during the 
scattering of background photons. Such extremely energetic photons produce pairs on the 
zero-th Landau level (||-polarised photons) or split to less energetic photons before reaching 



the first threshold (see Section 4.2.6 ). After the photons split the resulting photons are still 
very energetic and create an electron-positron pair enhancing the avalanche production 
of particles. In contrast to the PSG-off, most of the electron-positron pairs in the PSG- 
on mode are created well below ZPF. Furthermore, there is no additional radiation at 



lower energies (e < 1 MeV) which could be detected by a distant observer. Figures 5.3 



and |5.4| present the distribution of photons produced by the first population of newly 
created particles for PSR B0950-I-08 and PSR B1929-I-10, respectively. In both cases the 
energy of the 7-photons ranges from 1 GeV to ~ 20 GeV. The narrow predicted spark 
half-width of PSR B0950-I-08 results in a lower potential in lAR, thus increasing the 
efficiency of ICS (more photons produced by the first population of particles). The particle 



mean free path for ICS is smaller for backstreaming particles (see Section |3.2.1 for more 



details), thus most photons in the PSG-on mode are produced in the direction towards the 
stellar surface. Note that not all photons will produce electron-positron pairs since some 
7-photons are produced so close to the stellar surface that they reach its surface before 
they manage to reach the first threshold for pair production. 
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Figure 5.3: Distribution of photons produced in lAR by the first population of newly created 
particles for PSR B0950+08. In the calculations we used the parameters of the gap in the PSG-on 



mode as presented in Table 3.1 
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Figure 5.4: Distribution of photons produced in lAR by the first population of newly created 
particles for PSR B1929+10. In the calculations we used the parameters of the gap in the PSG-on 



mode as presented in Table 3.1 



5.1.2 X-ray and less energetic emission 

An negligible fraction of energy radiated by a primary particle in the PSG-off mode falls 
in the X-ray band. What is more, in the PSG-on mode all photons produced by ICS have 
energy which exceeds an electron's rest energy by many orders of magnitude. Thus, lAR 
may be responsible only for generating the thermal component of the X-ray spectrum in 
the process of heating the stellar surface. 

5.1.2.1 Thermal emission 



As shown in Section 1.4 , thermal emission is a common feature of neutron stars. Due to the 



large uncertainties in X-ray observations, it is not possible to distinguish all three thermal 
components (entire surface radiation, warm spot component and hot spot radiation) for 
one specific pulsar. Furthermore, only for a few pulsars (e.g. Geminga, PSR B0656+14) 
was it possible to distinguish two thermal components alongside the nonthermal one. In 
this thesis we focus on an analysis of pulsars with a visible hot spot component (6 > 1), 
since only for these pulsars is it possible to estimate the size of the actual polar cap. Most 
of these pulsars are old neutron stars and only for one of them (Geminga) was the whole 



surface radiation found in the X-ray spectrum. Figure |5.5| presents the observed X-ray 
components of the Geminga pulsar: the whole surface radiation, the polar cap (hot spot) 
and the nonthermal component. The maximum of energy for the whole surface radiation 
is in extreme ultraviolet and in soft X-rays for the hot spot component. Taking into 
account the very small area of the polar cap, radiation from the hot spot is unlikely to be 
observed in wavelengths off the maximum. 
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Figure 5.5: Observed flux of radiation for PSR J0633+1746. In the figure we present three 
components of radiation: the nonthermal one (green line), the entire surface radiation (blue line), 
and the hot spot component (red line). The dashed lines correspond to uncertainties in observations 



(see Table 1.4) 



Figure 5.6 presents the X-ray spectrum of PSR B1133+16. The small number of counts 



detected resulted in the fact that only separate fits for the BB and PL components were 
performed. Both the BB and PL fits describe the observed spectrum with similar accuracy. 
In the Figure we present additional thermal components (the entire surface radiation and 
the warm spot) which have not been determined by the observations. The Figure shows 
that the overlapping thermal components can mimic the power-law dependence of the 
spectrum at frequencies below 2 keV. 
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Figure 5.6: X-ray spectrum of PSR 61133+16. In addition to the observed thermal radiation (red 
solid line), two other thermal components are presented: the warm spot radiation (green dashed line) 
and the entire surface radiation (blue dotted line). 
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Although this specific combination of thermal components for PSR B1133+16 would 
result in a photon index greater than the observed one F = 2.51, the spectral fits for all 
pulsars should be extended to include more BB components in order to examine the effect 
of thermal components overlapping at lower frequencies. The results of our calculations 
suggest that the nonthermal X-ray radiation should dominate the spectrum at higher 
frequencies ^ 3 — lOkeV, but the power-law-like behaviour at lower frequencies could be 
the result of the overlapping of thermal components anticipated in the PSG scenario (see 



Section 5.1.2.3). 



5.1.2.2 Nonthermal emission 

The polarisation of ICS radiation in an ultrastrong magnetic field is 50% (one || to every 
±-polarised photon). Synchrotron Radiation of secondary particles created by ±-polarised 
photons would generate hard X-ray photons, however, as was mentioned in Section 4.2.6[ 



these photons will split before they reach the first threshold to produce pairs. Therefore, 
regardless of whether the gap is dominated by CR or by ICS, Synchrotron Radiation in 
lAR is not significant. 

5.1.2.3 Warm spot component 

Apart from the obvious X-ray component corresponding to the whole surface radiation, the 
PSG model can explain both the hot and warm spot radiation. The hot spot radiation is a 
natural consequence of heating the actual polar cap region by the backstreaming particles 



(see Section 1.4.4). As was mentioned in Section 1.4.4, the warm spot component can have 
two different sources: (1) the drastic difference of the crustal transport process due to the 
non-dipolar structure of the surface magnetic field (for young and middle-aged pulsars), 
(11) and a mechanism of heating the surface adjacent to the polar cap. In this section 
we present the second mechanism, i.e. heating of the surface adjacent to the polar cap, 
which can be applied to both young and old pulsars. 



Figure |5.7| presents the mechanism of heating the area adjacent to the polar cap for 
PSR B0950+08. When the gap operates in the PSG-off mode the primary plasma (see 



Section 5.1.3) will lose a significant part of its energy via CR as the particles propagate 
through the region of high curvature. For this particular magnetic line's configuration 
the region of high CR extends up to an altitude about 4 km above the stellar surface. 
The most energetic CR photons emitted in this region have a relatively short mean free 
path and they produce electron-positron pairs in the region of open magnetic field lines. 
However, both the less energetic CR photons and 7-photons produced by SR have a large 
enough photon mean free path to produce pairs in the region of the closed magnetic field 
lines. All newly created pairs move along the closed magnetic field lines and heat the 
surface beyond the polar cap on the opposite side of the star. 
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Figure 5.7: Global structure of magnetic field lines for PSR B0950+08. The structure was obtained 



using two crust-anchored anomalies (see Section 2.4.5). Green lines correspond to the outer open 
magnetic field lines, while the red lines correspond to the closed magnetic field lines at which secondary 
pairs are produced. Blue, yellow and red dots represent the locations of secondary pair production 
for the outer left, the middle and the outer right open field lines, respectively. 



The fraction of energy transferred to the region of the closed field lines highly depends 
on the region of open magnetic field lines considered in CR/SR emission. In the Figure 
we use three different colours (blue, yellow and red) to show the positions of pair creation 
for three characteristic open magnetic field lines (the outer left, the middle and the outer 
right). The simulation results in the following fractions of energy transferred to the 
region of the closed field lines are: 0.02%, 0.1%, 6% for all three lines, respectively. For 
this specific magnetic field configuration the transferred energy fraction increases as we 
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move towards the region with the highest curvature. We can roughly estimate that for 
the proposed magnetic field configuration of PSR B0950+08, about 1% of the outfiowing 
energy is responsible for heating of the surface beyond the polar cap on the opposite side 
of the star. Note that due to strong anisotropy of the outfiowing and backfiowing stream 



of particles (see Section 5.1.3), this fraction could be enough to obtain the warm spot 
component with a luminosity equal or in some cases even higher than the luminosity of 
the hot spot component. 
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Figure 5.8: Global structure of magnetic field lines for PSR B0943+10. The structure was obtained 



using two crust-anchored anomalies (see Section 2.4.4). Green lines correspond to the outer open 
magnetic field lines, while the red lines correspond to the closed magnetic field lines at which secondary 
pairs are produced. Blue, yellow and red dots represent the locations of secondary pair creation for 
the outer left, the middle and the outer right open field lines, respectively. 
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The fraction of energy transferred to the region of closed field lines highly depends 
on the magnetic field configuration. A more complicated structure of the magnetic field 



lines proposed for PSR B0943+10 (see Figure 5.8) results in a much wider area of closed 
field lines at which pairs are created, and hence higher fractions of energy transferred 
to the region of closed field lines. For this specific structure of the magnetic field these 
fractions are: 7%, 1%, 2% for three characteristic lines, respectively. We can roughly 
estimate that about 3 — 5% of the outflowing energy is responsible for the heating. Note 
that the magnetic field structure of PSR B0950+08 results in the heating of only one side 
beyond the polar cap, while in the case of PSR B0943+10 the whole surface around the 
polar cap is heated. The actual size of the warm spot also depends on the magnetic field 
configuration in the heating zone, and can either be decreased or increased. 

5.1.3 Primary plasma 



As we mentioned in Section 3.2.2, PSG-off and PSG-on modes differ essentially by the 
Lorentz factor of primary particles produced in the gap region. Furthermore, different 
scenarios of the gap breakdown (due to surface overheating or due to production of dense 
enough plasma) cause the evolution of primary particles in the two modes to completely 
different. 

We assume that in the PSG-off mode the gap breakdown is due to surface overheating; 
hence the plasma cloud moving away from the stellar surface is a mixture of ions and 
electron-positron plasma. In this scenario the ions are the main source of charge density 



required to screen the gap (see Equation 3.40). As the plasma cloud moves away from 



the stellar surface both the spark height and the spark width increase, which results 
in an increase of the acceleration potential drop. When the particles gain the Lorentz 
factors 7 > 10^, CR begins to produce 7-photons. In the PSG-off mode most of the 



7-photons are created near ZPF (see Figure 3.4). All particles created by 7-photons above 



the ZPF do not contribute to the heating of the surface. Furthermore, the acceleration 
in the upper parts of the gap is relatively weak, and electrons produced in this region 
will also escape from the gap, thus not contributing to the surface heating. Depending 
on the details of the cascade formation, the process described above may result in the 
creation of strong streaming anisotropics, where the flux of the backstreaming particles 
is considerably smaller than the flux of the outstreaming particles. Note that the density 
of the backstreaming particles required to overheat the surface is significantly lower than 



the co-rotational density ncR <^ ^gj (sgc Table 3.1 ). 

In the PSG-on mode the quasi-equilibrium of the flux of backstreaming particles and 
the flux of the polar cap radiation can cause the gap to break only due to the produc- 
tion of dense enough plasma. Thus, the surplus of positrons is the main source of the 
charge in the plasma cloud moving away from the stellar surface. The ICS process re- 
sponsible for the cascade production of particles is effective only in the bottom part of 
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the gap. Hence, the backstreaming electrons will hit the surface with a Lorentz factor 7c 
well below the 7max- As there is no strong pair production near (or above) ZPF, the back- 
streaming/outstreaniing anisotropy arises only due to the difference of the Lorentz factor 
of electrons hitting the stellar surface and the Lorentz factor of positrons accelerated in the 
gap 7max/7c ~ 10. The actual density of newly created plasma to completely screen the 



gap can be calculated only in the full cascade simulation. However, as shown by Timokhin 



(2010), this density should significantly exceed the co- rotational Goldreich- Julian density 
nics ^ ^GJ- We describe the difference between the co-rotational density and the actual 
density of primary plasma required to completely screen the ICS-dominated gap by factor 
A^ics = nics/riGj > 1. 



5.2 Inner magnetosphere of a pulsar 

5.2.1 Gamma- ray emission 

In general there are three processes which can produce 7-ray emission in the inner mag- 
netosphere of a pulsar {Rpc <^ z <^ Rlc)'- CR, ICS and SR. Which of them produces 
the majority of 7-photons depends on the parameters of the primary particles, and thus 
mainly depends on the mode in which the gap operates. Additionally, the efficiency of 
the ICS process strongly depends on the source of the background photons. 

5.2.1.1 Curvature Radiation of primary particles 

When the gap operates in the PSG-off mode, high-energetic particles are produced 
7c > 10^. As they pass the region with high curvature (SfJg ~ 1) they radiate a signi- 



ficant part of their energy through CR (see Section 4.1). 



Figure 5.9 presents the distribution of CR photons produced by a single primary 



particle moving along the open magnetic field line of PSR B1133+16 (see Section 2.4.5 
for the details of the magnetic field configuration). The initial Lorentz factor of the 
particle 7max = 1-7 x 10^ was set according to the value presented in Table 3.1 As the 



particle advanced through the region with high curvature, it lost about 46% of its initial 
energy, which was mainly converted to high-energetic 7-photons with an energy up to 
about 2 GeV. The 7-photons are produced in a region of a strong magnetic field, thus 
after passing a relatively short distance the most energetic photons are absorbed by the 
magnetic field and electron-positron pairs emerge. The red colour in the Figure corres- 
ponds to the final spectrum (after photon splitting, pair production and SR) produced 
by a single primary particle in the PSG-off mode. Most of the energy radiated by the 



primary particle was converted into the secondary plasma (see Section 5.2.3) and only 
about 5% of the particle's initial energy ended in the form of radiation with a cut-off at 
about 30 MeV. 
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Figure 5.9: Final photon distribution produced by a single primary particle for PSR B1133+16. The 
blue line corresponds to the initial CR distribution, while the red line presents the final distribution 
with the inclusion of photon splitting, pair production and SR. 



To increase the amount of photons reaching the observer, the emission zone, i.e. the 
region with the highest curvature, should by located in the area with a weaker magnetic 
field. Such a configuration allows a photon to travel a longer distance before it is absorbed 
by the magnetic field. As a result the electron-positron pairs are created at higher Landau 



levels, which enhances SR. Figure 5.10 presents the distribution of CR photons for PSR 



B0950+08. The calculations were performed for the initial Lorentz factor of the particle 
7max = 2.0 X 10^ (see Table Isll). 
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Figure 5.10: Final photon distribution produced by a single primary particle for PSR B0950+08. 
The blue line corresponds to the initial CR distribution, while the red line presents the final distribution 
with the inclusion of photon splitting, pair production and SR. 



Due to CR the primary particle lost about 40% of its initial energy. In this case 
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about a half of the energy radiated by the primary particle was converted into the second- 
ary plasma and the same amount of energy (about 20% of the particle's initial energy) 
ended in the form of radiation. For both PSR B1133+16 and PSR B0950+08, the max- 
imum of the curvature is of the same order. However, the maximum of curvature for 
PSR B1133+16 is located at an altitude of about 800 m, while for PSR B0950+08 it is 



located at an altitude of about 1.75 km (compare Figures 2.23 and 2.20) 



5.2.1.2 Inverse Compton Scattering of primary particles 



In the PSG-on mode the maximum Lorentz factor of primary particles is in the range of 



10^ — 10^ (see Table 3.1). As it follows from Figures 4.17 and 4.18, the ICS process is 
most effective for particles with a Lorentz factor in the range of 10^ — 10^. Particles with 
high energies (7 > 10^) will upscatter thermal photons only just above the stellar surface. 



where the density of the background photons is very high (see Section 5.1.1.2). Thus, if 
there is no additional source of background photons, the most energetic particles (7 > 10^) 
will escape from the inner magnetosphere without losing their energy by ICS. However, 
the plasma cloud produced by the ICS-dominated gap has a density exceeding the co- 



rotational Goldreich- Julian density even by a few orders of magnitude (see Section 5.1.3). 



Such a high charge density reduces the acceleration (Timokhin, 2010) and, consequently 



the bulk of particles will escape from the lAR with lower Lorentz factors. It is not possible 
to estimate the actual Lorentz factor of particles in the plasma cloud at the moment of gap 
breakdown without performing a full cascade simulation. Thus, in this thesis we assume 
that at the moment of gap breakdown most of the particles will have an energy that is 
about the characteristic value at which the acceleration is stopped by ICS in the bottom 
parts of the lAR 7c. To increase readability for cascade simulations with very low surface 
temperature, in all the Figures of the ICS distribution we present 7-photons produced by 
50 primary particles with Lorentz factors in the range of 0.57c — 27c. 



In Figure [5^11] we present the distribution of ICS photons produced by the upscattering 
of surface thermal radiation with temperature Tg = 0.3 MK for PSR B0834-I-06. Even 
for such a low surface temperature the whole surface radiation is the dominant source 
of background photons for ICS up to altitudes of about one stellar radii. During the 
scattering the primary particles lose about 30% of their initial energy while producing 
7-photons with energy up to 1 GeV. Since the 7-photons are very energetic and are 
produced in a region with a strong magnetic field, they will be absorbed by the magnetic 



field, thus giving rise to the secondary plasma population (see Section 5.2.3). All pairs 



in the inner magnetosphere of a pulsar are created in the nonzero Landau level, thus the 
pair production process is also accompanied by strong SR (see the next section). 
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Figure 5.11: Distribution of ICS photons produced by an upscattering of surface radiation 
(Ts = 0.3 MK) for PSR B0834+06. The plot includes all 7-photons upscatterd by 50 primary 
particles with Lorentz factors in the range of 2.5 x 10^ — 10^. 



A natural way of increasing the number of 7-photons produced by ICS in the inner 



magnetosphere is to increase the number of background photons. Figure 5.12 presents the 
distribution of ICS photons produced by an upscattering of the surface thermal radiation 
with temperature Ts = 0.4 MK for PSR B0834+06. During the ICS process the primary 
particles lose about 65% of their initial energy. For higher surface temperatures the ICS 
produces 7-photons up to higher altitudes (about two stellar radii), thus photons with 
lower energy emerge emin ~ 3 MeV. These less energetic photons will reach the observer, 
but their total energy is significantly lower than the total energy of the secondary plasma 
created by more energetic 7-photons. 
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Figure 5.12: Distribution of ICS photons produced by an upscattering of the surface radiation 
(Ts = 0.4 MK) for PSR B0834+06. The plot includes all 7-photons upscatterd by 50 primary 
particles with Lorentz factors in the range of 2.5 x 10^ — 10^. 
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Note that although for PSR B0834+06 the X-ray spectral fit was performed with 
only one BB component, the surface temperatures used in the calculations (0.3 MK and 
0.4 MK) are in good agreement with the predicted surface temperature of an old neutron 
star. 

Another source of background photons which could be relevant for ICS in the inner 



magnetosphere is the warm spot component. As shown in Section [5.1.2.3[ if the antipodal 
spot operates in the PSG-off mode and if the magnetic field structure is suitable then 
the warm spot is formed in the region adjacent to the polar cap. With a temperature 
lower than the hot spot but a much larger area, the warm spot is the main source of the 
background photons at altitudes up to about half a stellar radius. 



In Figure [5T3] we present the distribution of ICS photons produced by an upscattering 
of warm spot radiation with temperature Tg = 1.0 MK and radius R„s = 1km for PSR 
B0834-I-06. When the warm spot is the main source of background photons, the ICS 
process starts at lower altitudes. As a consequence, the scattering produces photons 
with higher energy and the primary particles lose up to 90% of their initial energy. All 
these high energetic 7-photons are absorbed by the magnetic field producing electron- 
positron pairs. Note that for this specific pulsar the existence of such a strong warm spot 



component is unlikely, but as mentioned in Section |5.1.2.1| the X-ray spectral fits should 
be extended to include more thermal components to put better constraints on the X-ray 
emission of pulsars. 
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Figure 5.13: Distribution of ICS photons produced by an upscattering of warm spot radiation 
(Ts = 1.0 MK, i?ws = 1km) for PSR B0834+06. The plot includes all 7-photons upscatterd by 50 
primary particles with Lorentz factors in the range of 2.5 x 10'^ — 10^. 



5.2.1.3 Synchrotron Radiation 

In both PSG-off and PSG-on modes SR plays a significant role in the generation of soft 



7-ray photons. Figure 5.14 presents the places at which SR-photons are generated (left 
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panel) and the SR spectrum (right panel) in the PSG-off mode for Geminga. The most 
energetic photons are generated close to the stellar surface {z ^ 500 m), while the less 
energetic ones are produced at altitudes z > 2 km, where the magnetic field is weaker. SR 
in the PSG-off mode produces photons with energy in the range of from 30keV to 1 GeV. 
Again, the high energetic 7-photons produce electron-positron pairs in a strong magnetic 
field, thus its observation is not possible. 
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Figure 5.14: Synchrotron Radiation in the PSG-ofF mode for PSR J0633+1746. The left panel 
presents the places at which SR-photons are generated, while the right panel presents the SR-photons 
distribution. Plots were obtained in a cascade simulation calculated for a single primary particle 
moving along the extreme left open magnetic field line. 



In Figure 5.15 we present the places of SR-photon generation (left panel) and the 
energy distribution of photons (right panel) in the PSG-on mode for Geminga. 
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Figure 5.15: Synchrotron Radiation in the PSG-on mode for PSR J0633+1746. The left panel 
presents places at which SR-photons are generated, while the right panel presents the SR-photons 
distribution. Plots were obtained in a cascade simulation calculated for a single primary particle 
moving along the extreme left open magnetic field line. The ICS process was calculated using the 



whole surface radiation with temperature Tg = 0.5 MK (see Table 1.4) 
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The production of SR-photons in the PSG-off mode starts at altitudes about z ^ 1 km 
and ends at altitudes z ~ 4.5 km. Thus, the energy range of SR-photons is narrower, 
with the minimum and maximum photon energy Cmm ~ 40keV and emax ~ 50MeV, 
respectively. Note the significant difference in the number of photons produced by SR in 
the PSG-off and PSG-on modes. The difference is a direct consequence of low secondary 



plasma multiplicity in the PSG-on mode (see Section 5.2.3) 



5.2.2 X-ray emission 

The main source of X-ray photons produced in the inner magnetosphere is SR. As men- 



tioned in Section |5.2.1.1[ to increase the amount of photons reaching the observer the 
emission zone should by located in the area with a weaker magnetic field. In this sec- 
tion we focus on the results of PSR B0943-I-10 and PSR 1929-1-10 for which the proposed 



configuration of a magnetic field satisfies this requirement (see Sections 2.4.4 and 2.4.7) 



In the PSG-off mode most of the X-ray photons are produced by the SR of newly cre- 



ated electron-positron pairs. Figure 5.16 presents the final photon distribution produced 
by a single primary particle of PSR B1929-I-10 in the PSG-off mode. For a single primary 
particle we can estimate that only about 0.7% of the total photon energy is in the range 
of 1 — lOkeV. The bulk of the energy is carried away by newly created particles (73%) 
and high energetic photons (27%). 
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Figure 5.16: Final photon distribution produced by a single primary particle for PSR 1929-1-10 in 
the PSG-off mode. The blue line corresponds to the initial CR photons distribution, while the red 
line presents the final distribution with the inclusion of photon splitting, pair production and SR. 



In Figure 5.17 we present the locations and the photon distribution of SR for PSR 
1929-f 10 in the PSG-off mode. All SR-photons produced closer to the stellar surface 
will contribute to 7-ray emission, while the SR-photons produced at higher altitudes will 
produce photons in the X-ray band. As it results from Figure 2.26[ the curvature at an 
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altitude of 2; ~ 3.5 km is only 50% higher than at z ~ 2 km. Furthermore, before the 
particle reaches the region with a relatively low magnetic field {z ~ 3.5 km), it radiates a 
significant part of its energy at lower heights. 
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Figure 5.17: Synchrotron Radiation in the PSG-off mode for PSR B1929+10. The left panel 
presents the places at which SR-photons are generated, while the right panel presents the SR- 
photons distribution. Plots were obtained in a cascade simulation calculated for a single primary 
particle moving along the extreme left open magnetic field line. 



To increase radiation in the 1 — lOkeV energy band we should apply the magnetic 
field structure with considerably higher curvature at altitudes where X-ray photons are 
generated. Although the curvature will not directly affect the SR, it will enhance CR, 
and thus it will increase the number of pairs produced in the region of a relatively weak 



magnetic field. In Figure 5.18 we present the final photon distribution produced by a single 
primary particle for PSR B0943+10 calculated using the magnetic field configuration as 



presented in Section 2.4.4 




10 10 



10 



10" 10 10' 

e [MeV] 



10 



10 10 



Figure 5.18: Final photon distribution produced by a single primary particle for PSR B0943+10 
in the PSG-off mode. The blue line corresponds to the initial CR photons distribution, while the red 
line presents the final distribution with the inclusion of photon splitting, pair production and SR. 
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For this magnetic field structure about 3% of the total photon energy is in the range 
of 1 — lOkeV. The newly created particles carry away about 63% of the energy radiated 
by the primary particle, while about 37% of the energy remains in the form of photons. 



The structure of the magnetic field of PSR B0943+10 allows enhanced pair production 



in a region of a weaker magnetic field (see Figure 5.19). The SR that accompanies pair 
production at higher altitudes {z > 3 km) essentially increases the amount of energy 
radiated in the 1 — lOkeV energy band. Note, however, that the fraction of energy 
radiated in this band is still relatively low (3%), and in order to be a substantial part of 
the observed X-ray spectrum the strong anisotropy of backstreaming and outstreaming 



plasma is required (see Section 5.1.3). 
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Figure 5.19: Synchrotron Radiation in the PSG-ofF mode for PSR B0943+10. The left panel 
presents the places at which SR-photons are generated, while the right panel presents the SR- 
photons distribution. Plots were obtained in a cascade simulation calculated for a single primary 
particle moving along the extreme left open magnetic field line. 



In the PSG-on mode even for a complicated structure of the magnetic field most of 



the outflowing energy is converted to secondary plasma. Figure |5.20| presents the ICS- 
photons distribution produced in the PSG-on mode for PSR B0628-28. The bulk of 
energy is radiated in the form of high energetic 7-photons which are responsible for pair 
production, and thus the formation of secondary plasma. Taking into account not so 
high backstreaming/outstreaming anisotropy in the PSG-on mode, the ICS process is not 
relevant for the production of X-ray photons. 
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Figure 5.20: Distribution of ICS photons produced by an upscattering of the whole surface radiation 
(Ts = 0.5 MK) for PSR B0628-28. The plot includes all photons upscatterd by 50 primary particles 
with Lorentz factors in the range of 3 x 10^ — 1.2 x 10^. 



The SR which accompanies the pair creation process in the PSG-on mode mostly 



produces soft 7-photons (see the right panel of Figure |5.21 ). Although the secondary 
pairs are produced at similar altitudes in both modes, (compare the left panels of Figures 



5.19 and 5.21), the higher Lorentz factor of secondary plasma produced in the PSG- 



on mode results in higher energy of the SR-photons. The results suggest that when 
the gap operates in the PSG-on mode we should expect lower efficiencies of nonthermal 
X-ray emission than in the PSG-off mode. Note, however, that the final efficiency of 
X-ray radiation in the PSG-off mode highly depends on the backstreaming/outstreaming 
anisotropy and the structure of magnetic field lines. 
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Figure 5.21: Synchrotron Radiation in the PSG-on mode for PSR B0628-28. The left panel 
presents places at which SR-photons are generated, while the right panel presents the SR-photons 
distribution. Plots were obtained in a cascade simulation calculated for 50 primary particles with 
Lorentz factors in the range of 3 x 10^ — 1.2 x 10^ moving along the extreme left open magnetic 
field line. 
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5.2.3 Secondary plasma 

The multipHcity of secondary particles in the PSG-off mode is much higher than in the 
PSG-on mode. However, the primary plasma produced in the JAR of CR-dominated 
gaps has a density considerably lower than the Goldreich- Julian co-rotational density 



see Equation 13.40). Figure 5.22 presents the energy histogram of secondary plasma for 



Geminga (left panel) and PSR B1133+16 (right panel). Despite major differences in the 
magnetic field structure and conditions in the JAR for both pulsars, the secondary plasma 
distribution shows many similarities. The only significant difference is the maximum 
Lorentz factor of secondary plasma, which for Geminga is about 7™^^ ~ 10^, while for 
PSR B1133+16 is is a few times smaller 7^^,^^ ^3x lOl 



By using the overheating parameters presented in Table |3.1| we can roughly estimate 
that the final multiplicity of particles in the plasma cloud in the PSG-off mode ranges 
from M = K ■ Msec - 2 (for Geminga) to M = k ■ Msec ~ 100 (for PSR B1133+16). 
Note, however, that these values do not take into account the anticipated anisotropy of 
backstreaming and outstreaming particles. The existence of such an anisotropy could 
further increase the final multiplicity of particles in the plasma cloud leaving the inner 
magnetosphere. Despite the fact that without a full cascade simulation in the lAR we 
cannot unambiguously determine the final multiplicity in the plasma cloud, it can be 
clearly seen that depending on the details of the gap operating in the PSG-off mode, the 
produced plasma may be suitable (e.g. PSR B1133+16) or unsuitable (e.g. Geminga) to 



generate radio emission (see Section 5.2.3 ). The main factor determining the parameters of 
the CR-dominated gap, and thus determining whether it is possible to effectively produce 



radio emission, is the radius of curvature of the magnetic field lines (see Section 3.3.2) 
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Figure 5.22: Energy histogram of secondary plasma in the PSG-off mode. The left panel was 
obtained in a cascade simulation calculated for a single primary particle moving along the extreme 
left open magnetic field line of PSR J0633-I-1746, while the right panel corresponds to a cascade 
simulation for PSR B1133+16. 



In ICS-dominated gaps, on the other hand, the density of primary plasma produced 
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in lAR exceeds the co-rotational density. Thus, the development of dense enough plasma 



for radio emission is much easier in the PSG-on mode. In Figure 5.23 we present the 
locations of pair production and energy distribution of secondary plasma in the PSG-on 
mode for PSR B0628-28. The final multiplicity of particles in the plasma cloud in the 
PSG-on mode can be calculated as M = Nir,^ x M^pr. As mentioned in Section 15. 1.31 the 



exact value of A^ics can be found only by performing the full cascade simulation in lAR 



but, as shown by Timokhin (2010), we should expect a full screening of the acceleration 



region when Nics reaches a value as high as 20 — 100. Thus we can roughly estimate that 
for the whole surface radiation with temperature Tg = 0.3 MK, the final multiplicity of 
secondary plasma in the PSG-on mode for PSR B0628-28 is of the order of M ;^ 100. 
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Figure 5.23: Secondary plasma produced in the PSG-on mode for PSR B0628-28. The left panel 
presents places at which pairs are produced, while the right panel presents the histogram of particle 
energy. Plots were obtained in a cascade simulation calculated for 50 primary particles with Lorentz 
factors in the range of 3 x 10'^ — 1.2 x 10'^ moving along the extreme left open magnetic field line. 
The ICS process was calculated using the whole surface radiation with temperature Tg = 0.3 MK. 



In the PSG-on mode the main factor which determines the final multiplicity of second- 
ary plasma is the source of the background photons. As shown in Section 4.4.4[ the polar 
cap radiation (the hot spot component) has a negligible impact on the ICS process above 



the lAR. Figure |5.23| presents the location of pair production and energy distribution 
of secondary plasma for PSR B0628-28 calculated assuming the whole surface radiation 
with temperature Tg = 0.5 MK. The increase in the number of background photons res- 
ults in enhancement of the ICS process, and thus an increase of the secondary multiplicity 
Msec ~ 60. For such conditions the final multiplicity of secondary plasma in the PSG-on 
mode is of the order of M ^ 10^ — 10^. 
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Figure 5.24: Secondary plasma produced in the PSG-on mode for PSR B0628-28. The left panel 
presents the places at which pairs are produced, while the right panel presents the histogram of particle 
energy. Plots were obtained in a cascade simulation calculated for 50 primary particles with Lorentz 
factors in the range of 3 x 10^ — 1.2 x 10^ moving along the extreme left open magnetic field line. 
The ICS process was calculated using the whole surface radiation with temperature Tg = 0.5 MK. 



Conclusions 



The hot spot component identified in X-ray observations imphes the non-dipolar structure 
of surface magnetic field. We used the Partially Screened Gap model to explain both the 
X-ray radiation of radio pulsars and production of secondary plasma suitable for genera- 
tion of radio emission. 



A special case (PSR B0943+10) 

Our model predicts two additional sources of X-ray emission: (I) the warm spot compon- 
ent and (II) enhanced SR radiation in the PSG-off mode. The warm spot component is 
associated with particles originating from the antipodal polar cap, while the high lumin- 
osity of X-ray photons produced in the PSG-off mode is a result of strong anisotropy of 
backstreaming and outstreaming particles. 



Very recent results presented by Hermsen et al. (2013) show the anti-correlation of ra- 



dio and X-ray emission of PSR B0943-I-10. The authors suggest an unpulsed, non-thermal 
component in radio-bright mode and a 100%-pulsed thermal component along with a non- 
thermal component in a radio-quiet mode. In our model it is not possible to produce an 
unpulsed, nonthermal X-ray component without the accompanying blackbody radiation 
of the polar cap. Although it is possible to produce nonthermal X-ray radiation which 
obscures the thermal component (strong SR in the PSG-off mode with a high predom- 
inance of outstreaming particles), the resulting radiation should be pulsed. We believe 
that the X-ray radiation of PSR B0943-I-10 in the radio bright mode was misinterpreted 



as the nonthermal one. As shown in Figure 1.7 (panel d), for a derived geometry of PSR 



B0943-(-10 the polar cap produces unpulsed, thermal radiation. Furthermore, as repor- 
ted by the authors, in the radio-bright mode both the absorbed blackbody (BB) and the 



absorbed power-law (PL) models fit the spectrum equally well (see Table S4 in Hermsen 



et al. , 2013). We believe that the observed radiation modes of PSR B0943-I-10 correspond 
to a mode switch between the PSG-on (radio-bright) and the PSG-off mode (radio-quiet). 
When pulsar is in the PSG-on we observe both the radio emission and thermal radiation 
which originates from the polar cap. In the PSG-off mode the secondary plasma is not 
suitable to produce so strong radio emission as in the PSG-on mode, but the polar cap 



radiation is accompanied by pulsed, nonthermal emission produced by SR (see Sec. 5.2.2). 
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Gamma-ray pulsars 



As was shown in Sections 5.1.1 and 5.2.1 7-rays produced in lAR and the inner mag- 
netosphere cannot reach the observer due to efficient pair production in those regions. 
Current models of 7-ray emission propose that the emission comes from outer magneto- 
spheric gaps. The non-dipolar structure of a magnetic field has two key implications 
on 7-ray emission models: (I) the formation of slot gaps is not possible as pairs are pro- 
duced along all open magnetic field lines, (II) the high density of electron-positron plasma 
{up ^ ^Gj) produced in the inner magnetosphere prevents the outer gap formation. The 
high-density plasma which crosses the null line will screen the outer magnetospheric region 
due to plasma separation (acceleration of electrons and deceleration of positrons). Thus, 
the formation of outer gaps is possible only in special cases when the pulsar operates in 
the PSG-off mode and produces secondary plasma with low density np ^ ncj. 



As recently reported by Arka and Dubus (2013): "It is possible for relativistic popu- 



lations of electrons and positrons in the current sheet of a pulsar's wind right outside the 
light cylinder to emit synchrotron radiation that peaks in the sub — GeV to GeV regime, 
with 7-ray efficiencies similar to those observed for the Fermi/LAT pulsars." We believe 
that the observed high-energetic 7-rays are produced in the not yet well explored region 
right outside the light cylinder. 

Radio emission 

Pulsed radio emission remains one of the most intriguing puzzles of astrophysics. It is re- 
markable that despite the large ranges in P, B^, the variations in the pulse profile between 
different classes of neutron stars (young, old, millisecond, magnetars) are similar to those 



within classes (Melrose, 2004). The radio emission of most pulsars can be characterised 
by: a relatively narrow frequency range, ~ 100 MHz to ~ 10 GHz, and a high degree of 
polarisation with a characteristic sweep of the position angle. The extremely high bright- 
ness temperature of pulsar radio emission (typically T^ > 10^^ K) implies that a coherent 
emission mechanism is involved. Many radio emission mechanisms have been proposed, 
but no consensus on a specific emission mechanism has emerged. The radio observations 
alone cannot identify the emission mechanism and, hence, a model of the magnetosphere 
is needed to put constraints on the radio emission model. An acceptable emission mech- 
anism must involve some form of instability to produce coherent radiation. The main 
difficulty in finding a specific emission mechanism is that many of the predicted features 
are common all proposed models. Furthermore, the polarisation can also be regarded as 



generic rather than associated with a specific emission mechanism (Melrose, 2006). 

The X-ray observations have allowed us to put constraints on the polar cap region of 
pulsars. The non-dipolar structure of the surface magnetic field causes plasma to form 
under similar conditions regardless of the global configuration of the magnetic field. We 



5.2. Inner magnetosphere of a pulsar 149 

have showed that depending on the details of lAR, the resulting plasma either meets the 
requirements for efficient radio emission (suitable multiplicity and energy distribution of 
secondary plasma) or is not suitable to produce efficient radio emission (e.g. Geminga). 
Furthermore, the proposed drift model allows to find a connection between radio and 



X-ray emission processes (see Section 3.4). 



The mixed mode 

Although in the thesis we consider the PSG-on and PSG-off mode separately, in a real 
case both of these modes can coexist either on two separate polar caps or on the same 
polar cap occupying its different parts. In the latter case the change of modes is associated 
with varying degrees of intensity of the two modes. Furthermore, if specific conditions 
are met, the ICS process can be a main source of 7-photons in the lower parts of the gap, 
while the CR process can produce 7-photons in the upper parts of the acceleration region. 
In such a case distinguishing between the two modes is even more difficult. 

Summary 

The main propositions associated with this thesis are as follows: 

1. The size of the hot spots implies that the magnetic field configuration just above 
the stellar surface differs significantly from a purely dipole one. 

2. The analysis of X-ray observations shows that the temperature of the actual polar 
cap is equal to the so-called critical value, i.e. the temperature at which the outfiow 
of thermal ions from the surface screens the gap completely. 

3. The non-dipolar structure of a surface magnetic field and the high multiplicity of 
particles produced in lAR prevents the formation of slot and outer gaps. 

4. The PSG model predicts the existence of two scenarios of gap breakdown: the PSG- 
off mode for CR-dominated gaps and the PSG-on mode for ICS-dominated gaps. 

5. The two different scenarios of gap breakdown can in a natural way explain the mode- 
changing phenomenon when both modes produce plasma suitable to generate radio 
emission, and pulse nulling when the radio emission is not generated in one of the 
modes. 

6. The mode changes of the lAR may explain the anti-correlation of radio and X-ray 



emission in very recent observations of PSR B0943-I-10 (Hermsen et al. , 2013). 



7. The regular drift of subpulses can be expected only when the gap operates in the 
PSG-on mode. The proposed model of drift allows to connect the drift information 
obtained by radio observations with the X-ray data of rotation-powered pulsars. 
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